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Abstract 

The scattering problems of a scalar point particle from an assembly of 1 < n < oo non- 
overlapping and disconnected hard disks, fixed in the two-dimensional plane, belong to the 
simplest realizations of classically hyperbolic scattering systems. Their simplicity allows for a 
detailed study of the quantum mechanics, semiclassics and classics of the scattering. Here, we 
investigate the connection between the spectral properties of the quantum-mechanical scattering 
matrix and its semiclassical equivalent based on the semiclassical zeta-function of Gutzwiller and 
Voros. We construct the scattering matrix and its determinant for any non-overlapping n-disk 
system (with n < oo) and rewrite the determinant in such a way that it separates into the 
product over n determinants of 1-disk scattering matrices - representing the incoherent part of 
the scattering from the n-disk system - and the ratio of two mutually complex conjugate deter- 
minants of the genuine multiscattering matrix M which is of Korringa-Kohn-Rostoker-type and 
which represents the coherent multidisk aspect of the n-disk scattering. Our quantum-mechanical 
calculation is well-defined at every step, as the on-shell T -matrix and the multiscattering kernel 
M— 1 are shown to be trace-class. The multiscattering determinant can be organized in terms of 
the cumulant expansion which is the defining prescription for the determinant over an infinite, 
but trace-class matrix. The quantum cumulants are then expanded by traces which, in turn, 
split into quantum itineraries or cycles. These can be organized by a simple symbolic dynamics. 
The semiclassical reduction of the coherent multiscattering part takes place on the level of the 
quantum cycles. We show that the semiclassical analog of the mth quantum cumulant is the 
mth curvature term of the semiclassical zeta function. In this way quantum mechanics naturally 
imposes the curvature regularization structured by the topological (not the geometrical) length 
of the pertinent periodic orbits onto the semiclassical zeta function. However, since the cumulant 
limit m — > oo and the semiclassical limit, Ti — > or (wave number) k — > oo, do not commute 
in general, the semiclassical analog of the quantum multiscattering determinant is a curvature 
expanded semiclassical zeta function which is truncated in the curvature order. We relate the 
order of this truncation to the topological entropy of the corresponding classical system. We 
show this explicitly for the 3-disk scattering system and discuss the consequences of this trun- 
cation for the semiclassical predictions of the scattering resonances. We show that, under the 
above mentioned truncations in the curvature order, unitarity in n-disk scattering problems is 
preserved even at the semiclassical level. Finally, with the help of cluster phase shifts, it is shown 
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that the semiclassical zeta function of Gutzwiller and Voros has the correct stability structure 
and is superior to all the competitor zeta functions studied in the literature. 

PACS numbers: 03.65.Sq, 03.20.+i, 05.45.+b 
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1 Introduction 

The main focus of this manuscript is on the transition from quantum mechanics to semi- 
classics in classically hyperbolic scattering systems, and in particular, on the convergence 
problems of periodic orbit expansions of n-disk repellers. 

1.1 Motivation and historic perspective 

Why more than 70 years after the birth of textbook quantum mechanics and in the 
age of supercomputers is there still interest in semiclassical methods? First of all, there 
remains the intellectual challenge to derive classical mechanics from quantum mechanics, 
especially for classically non-separable chaotic problems. Pure quantum mechanics is linear 
and of power-law complexity, whereas classical mechanics is generically of exponential 
complexity. How does the latter emerge from the former? Secondly, in many fields (atomic 
physics, molecular physics and quantum chemistry, but also optics and acoustics which 
are not quantum systems but are also characterized by the transition from wave dynamics 
to ray dynamics) semiclassical methods have been very powerful in the past and are still 
useful today for practical calculations, from the detection of elementary particles to the 
(radar)-detection of airplanes or submarines. Third, the numerical methods for solving 
multidimensional, non-integrable quantum systems are generically of "black-box" type, 
e.g. the diagonalization of a large, but truncated hamiltonian matrix in a suitably chosen 
basis. They are computationally intense and provide little opportunity for learning how 
the underlying dynamics organizes itself. In contrast, semiclassical methods have a better 
chance to provide an intuitive understanding which may even be utilized as a vehicle for 
the interpretation of numerically calculated quantum-mechanical data. 

In the days of "old" quantum mechanics semiclassical techniques provided of course the 
only quantization techniques. Because of the failure, at that time, to describe more com- 
plicated systems such as the helium atom (see, however, the resolution of Wintgen and 
collaborators |lj|; @ and also H provide for a nice account of the history), they were 
replaced by modern quantum mechanics based on wave mechanics. Here, through WKB 
methods, they reappeared as approximation techniques for 1-dimensional systems and, 
in the generalization to the Einstein-Brillouin-Keller (EBK) quantization, for separable 
problems [[|-[5[] where an n-degree-of-freedom system reduces to n one-degree-of-freedom 
systems. Thus semiclassical methods had been limited to such systems which are classi- 
cally nearly integrable. 

It was Gutzwiller who in the late 60s and early 70s (see e.g. |J and 0) (re-)introduced 
semiclassical methods to deal with multidimensional, non-integrable quantum problems: 
with the help of Feynman path integral techniques the exact time-dependent propagator 
(heat kernel) is approximated, in stationary phase, by the semiclassical Van-Vleck propa- 
gator. After a Laplace transformation and under a further stationary phase approximation 
the energy- dependent semiclassical Green's function emerges. The trace of the latter is 
calculated and reduces under a third stationary phase transformation to a smooth Weyl 
term (which parametrizes the global geometrical features) and an oscillating sum over all 
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periodic orbits of the corresponding classical problem. Since the imaginary part of the 
trace of the exact Green's function is proportional to the spectral density, the Gutzwiller 
trace formula links the spectrum of eigen-energies, or at least the modulations in this 
spectrum, to the Weyl term and the sum over all periodic orbits. Around the same time, 
Balian and Bloch obtained similar results with the help of multiple-expansion techniques 
for Green's functions, especially in billiard cavities, see e.g. 0. 

For more than one degree of freedom, classical systems can exhibit chaos. Generically 
these are, however, non-hyperbolic classically mixed systems with elliptic islands embed- 
ded in chaotic zones and marginally stable orbits for which neither the Gutzwiller trace 
formula nor the EBK-techniques apply, see Berry and Tabor ||. Purely hyperbolic sys- 
tems with only isolated unstable periodic orbits are the exceptions. But in contrast to 
integrable systems, they are generically stable against small perturbations ||. Moreover, 
they allow the semiclassical periodic orbit quantization which can even be exact as for 
the case of the Selberg trace formula which relates the spectrum of the Laplace-Beltrami 
operator to geodesic motion on surfaces of constant negative curvature |J . The Gutzwiller 
trace formula for generic hyperbolic systems is, however, only an approximation, since its 
derivation is based on several semiclassical saddle-point methods as mentioned above. 

In recent years, mostly driven by the uprise of classical chaos, there has been a resurgence 
of semiclassical ideas and concepts. Considerable progress has been made by applying 
semiclassical periodic orbit formulae in the calculation of energy levels for bound-state 
systems or resonance poles for scattering systems, e.g., the anisotropic Kepler problem ||, 
the scattering problem on hard disks []Hj-[I3[], the helium atom etc. (See Ref.[TB] for 



a recent collection about periodic orbit theory.) It is well known that the periodic orbit 
sum for chaotic systems is divergent in the physical region of interest. This is the case on 
the real energy axis for bounded problems and in the region of resonances for scattering 
problems, because of the exponentially proliferating number of periodic orbits, see |]5|Ji7||. 
Hence refinements have been introduced in order to transform the periodic orbit sum in 
the physical domain of interest to a still conditionally convergent sum by using symbolic 
dynamics and the cycle expansion [|l^,|l^,|l^] , Riemann-Siegel lookalike formal and pseudo- 



orbit expansions p0y21| , surface of section techniques ||22|j23[| , inside-outside duality |2 



heat-kernel regularization 0,25] etc. These methods tend to be motivated from other 



areas in physics and mathematics [26] such as topology of flows in classical chaos, the 
theory of the Riemann zeta functions, the boundary integral method for partial differential 
equations, Fredholm theory (see also quantum field theory etc. 



In addition to the convergence problem, there exists the further complication for bounded 
smooth potential and billiard problems that the corresponding periodic orbit sums pre- 
dict in general non-hermitean spectra. This problem is addressed by the Berry-Keating 
resummation techniques pO| , pl| - however, in an ad- hoc fashion. In contrast, scattering 
problems circumvent this difficulty altogether since their corresponding resonances are 
complex to start with. Moreover, the scattering resonances follow directly from the pe- 
riodic orbit sum, as the Weyl term is absent for scattering problems. In fact, it is more 
correct to state that the Weyl term does not interfer with the periodic sum, as a negative 
Weyl term might still be present, see e.g.[l7j]. Furthermore, scattering systems allow for 



a nice interpretation of classical periodic sums in terms of survival probabilities 
In this respect, it is an interesting open problem why these classical calculations do not 
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seem to generate a Weyl term, whether applied for bounded or scattering systems. For 
these reasons, the study of periodic sums for scattering systems should be simpler than 
the corresponding study for bound-state problems, as only the convergence problem is 
the issue. 



1.2 The n-disk repeller: a model for hyperbolic scattering 



Hence, one should look for a simple classically hyperbolic scattering system which can be 
used to address the convergence problem. It should not be too special, as for example the 
motion on a surface of constant curvature, but reasonably realistic and instructive. Eck- 
hardt Jl0| suggested such a system, the "classical pinball machine" . It consists of a point 
particle and a finite number (in his case three) identical non-overlapping disconnected cir- 
cular disks in the plane which are centered at the corners of a regular polygon (in his case 
an equilateral triangle). The point particle scatters elastically from the disks and moves 
freely in between collisions. The classical mechanics, semiclassics and quantum mechanics 
of this so-called three-disk system was investigated in a series of papers by Gaspard and 
Rice, ||iT|-|T3|1, and, independently, by Cvitanovic and Eckhardt |T3|, see also Scherer (jT? 



and Ref . |L5] . It belongs to a class of mechanical systems which are everywhere defocusing, 
hence no stable periodic orbit can exist (see Fig.|l.l|). 



23132321 




2313 



Fig. 1.1. The three-disk repeller with the symbolic dynamics of the full domain. The figure is 
from Ref.@. 

The classical dynamics with one or two disks is simple, there is either no or one trapped 
trajectory. The latter is obviously unstable, since a small displacement leads to a defo- 
cusing after the reflection from the curved surface of disk p| . The two-disk system is 
therefore one of the simplest hyperbolic scattering systems, but it is non-chaotic. How- 
ever, with three or more disks there are infinitely many trapped trajectories forming a 
repeller [[OJ. The periodic orbits corresponding to these trapped trajectories are all iso- 
lated and unstable because of the defocusing nature of the reflections. Note that the one- 
disk and two-disk systems, although classically simple, are nonetheless interesting. The 
quantum-mechanical one-disk scattering system (since it is separable) has been one of the 
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key models for building up the semiclassical theory of diffraction p9|-|3T|. Similarly, the 
two-disk system became the toy ground for the periodic-orbit theory of diffraction |3^j33|| . 
In fact, the two-disk system has infinitely many diffractive creeping periodic orbits which 
can be classified by symbolic dynamics similarly to the infinitely many geometrical orbits 
of the three-disk system. The symbolic dynamics of a general n-disk system is very simple, 
see e.g. [0: periodic orbits can be classified by a series of "house numbers" of the disks 
which are visited by the point particle which follows the corresponding trajectory. Not all 
sequences are allowed: after each reflection from one disk, the point particle has to pro- 
ceed to a different disk, since the evolution between the disk is the free one. Furthermore, 
for general geometries there may exist sequences which correspond to trajectories which 
would directly pass through a disk. The sequences corresponding to these so-called "ghost 
orbits" have to be excluded from the classical consideration. In summary, the geometrical 
periodic orbits (including ghost orbits) are labelled in the full domain of the n-disk re- 
peller by itineraries (= periodic words) with n different symbols (=letters) with the trivial 
"pruning" rule that successive letters in the itinerary must be different. The itineraries 
corresponding to ghost orbits have to be removed or "pruned" with all their sub-branches 
from the symbol tree. Periodic trajectories which have reflections from inside of a disk 
(i.e. the point particle traverses first through a disk and is then reflected from the other 
side of the disk) can be excluded from the very beginning. In fact, in our semiclassical 
reduction of Sec.^ we will show for all repeller geometries with n non-overlapping disks 
that, to each specified itinerary, there belongs uniquely one standard periodic orbit which 
might contain ghost passages but which cannot be reflected from the inside. There is only 
one caveat: our method cannot decide whether grazing trajectories (which are tangential 
to a disk surface) belong to the class of ghost trajectories or to the class of reflected tra- 
jectories. For simplicity, we just exclude all geometries which allow for grazing periodic 
orbits from our proof. Alternatively, one might treat these grazing trajectories separately 



with the help of the diffractional methods of Refs. [pl| , |35 



The symbolic dynamics described above in the full domain applies of course to the equi- 
lateral three-disk system. However, because of the discrete C3 V symmetry of that system, 
the dynamics can be mapped into the fundamental domain (any one of the 1/6-th slices 
of the full domain which are centered at the symmetry-point of the system and which 



exactly cut through one half of each disk, see Fig.lO). 




Fig. 1.2. Equilateral 3-disk system and its fundamental domain. 
In this fundamental domain the three-letter symbolic dynamics of the full domain reduces 
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a two-letter symbolic dynamics. The symbol '0', say, labels all encounters of a periodic 
orbit with a disk in the fundamental domain where the point particle in the corresponding 
full domain is reflected to the disk where it was just coming from, whereas the symbol 
'1', say, is reserved for encounters where the point particle is reflected to the other disk. 

Whether the full or the fundamental domain is used, the 3-disk system allows for a unique 
symbolic labelling (if the disk separation is large enough even without non-trivial pruning). 
If a symbolic dynamics exist, the periodic orbits can be classified by their topological 
length which is defined as the length of the corresponding symbol sequence. In this case 
the various classical and semiclassical zeta functions are resummable in terms of the cycle 
expansion [0,|1^] which can be cast to a sum over a few fundamental cycles (or primary 
periodic orbits) tf and higher curvature corrections C m of increasing topological order m: 

S / m 

The curvature C m in fll.l]) contains all allowed periodic orbits of topological length m for a 
specified symbolic dynamics and suitable "shadow-corrections" of combinations (pseudo- 
orbits) of shorter periodic orbits with a combined topological length m. 

Common to most studies of the semiclassics of the n-disk repellers is that they are 
"bottom- up" approaches. Whether they use the Gutzwiller trace formula 0, the Ru- 
elle or dynamical zeta function |36| , the Gutzwiller- Voros zeta function |37| , their starting 



point is the cycle expansion [T^,|TP|J2^] . The periodic orbits are motivated from a semiclas- 
sical saddle-point approximation. The rest is classical in the sense that all quantities which 
enter the periodic orbit calculation as e.g. actions, stabilities and geometrical phases are 



determined at the classical level (see however Refs. |38|-^0[ where leading ^-corrections 
to the dynamical zeta function as well as the Gutzwiller- Voros zeta function have been 
calculated). For instance, the dynamical zeta function has typically the form 

( \E) = l[(l-t p ), tp = ~j= e'W*^- 1 * * ■ (1.2) 



The product is over all prime cycles (prime periodic orbits) p. The quantity A p is the 
stability factor of the p-th cycle, i.e., the expanding eigenvalue of the p-cycle Jacobian, S p 
is the action and v v is (the sum of) the Maslov index (and the group theoretical weight 
for a given representation) of the p-th cycle. For n-disk repellers, the action is simply 
S p = kL p , the product of the geometrical length L p of the periodic orbit and the wave 
number k = \JlmE in terms of the energy E and mass m of the point particle. The 
semiclassical predictions for the scattering resonances are then extracted from the zeros 
of the cycle-expanded semiclassical zeta-function. In this way one derives predictions of 
the dynamical zeta function for the leading resonances (which are the resonances closest to 
real /c-axis). In the case of the Gutzwiller- Voros zeta function also subleading resonances 
result, however, only if the resonances lie above a line defined by the poles of the dynamical 
zeta function |]4l|- ^3| , p!5| . The quasiclassical zeta function of Vattay and Cvitanovic is entire 
and gives predictions for subleading resonances for the entire lower half of the complex 
plane |43[ . 
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1.3 Objective 



As the n-disk scattering systems are generically hyperbolic, but still simple enough to 
allow for a closed-form quantum-mechanical setup [JD3] and detailed quantum-mechanical 
investigations P^JlSj], we want to study the structure of the semiclassical periodic sum 



for a hyperbolic scattering system in a "top-down" approach, i.e. in a direct derivation 
from the exact quantum mechanics of the n-disk repeller. This is in contrast to the usual 
semiclassical "bottom-up" studies of the n-disk repellers which can be affected by un- 
controlled operations during the long and mostly formal derivation from the Gutzwiller 
trace formula. Especially regularization prescriptions, like the cycle expansion, have to be 
added from the outside in order to get converging semiclassical predictions. 

Hence, for any n-disk scattering problem with a finite number of non-overlapping discon- 
nected disks we want to construct a direct link from the defining exact S-matrix to the 
pertinent semiclassics (in terms of a suitable periodic-orbit expansion) with the following 
qualifications: 

(i) The derivation should lead to a unique specification of the periodic orbits for a given 
n-disk geometry. The method should be able to handle n-disk geometries which allow 
for ghost orbits, i.e., periodic orbits existing in any of the pertinent "parent" disk- 
systems (defined by the removal of one or more disks) which are blocked by the 
return of at least one of the removed disks. 

(ii) Since disk-systems are known where the semiclassics is strongly governed by diffrac- 
tive orbits (see [Q and especially |k| for the 2-dimensional scattering analog of the 



two- well-potential problem), diffractive periodic orbits should emerge together with 
their standard partners, 
(iii) The subleading stability structure of the standard periodic orbits should follow from 
this derivation in order to discriminate between the Gutzwiller- Voros zeta function 



and other competitors, e.g., the dynamical zeta function of Ruelle or the qua- 
siclassical zeta function of Vattay and Cvitanovic [I3|; in other words, we want to 
derive the semiclassical spectral function. 

(iv) The setup of the starting-point, the quantum-mechanical side, should not be plagued 
by formal or uncontrolled manipulations or assumptions. Especially, if the quantum- 
mechanical side does not exist without a suitable regularization prescription, the 
latter should be provided before the semiclassical reduction is performed. This should 
exclude that the semiclassical sums encounter hidden problems which are already 
present at the quantum-mechanical level. 

(v) The link between the exact quantum mechanics and semiclassics should not only 
allow for the computation of scattering resonance, but should be valid for all values 
of the wave number, also away from the resonances and from the real axis, modulo the 
boundary of semiclassical convergence, as this issue can only be addressed during the 
link-procedure. Branch cuts and singularities on the quantum-mechanical side have 
to be taken into account of course. 

(vi) The spectral function should not only result in a formal sense, but, if necessary, with 
a pertinent regularization and summation prescription that should not be imposed 
from the outside. 

(vii) Most importantly, the derivation should be well-defined and allow for a test of the 
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summation prescription of the period-orbit expansion. If potential problems occur, 
they should be pinpointed in the derivation. 



1.4 Outline 



The manuscript is organized as follows. We begin in Sec.|2] with the standard approach 
relating quantum-mechanical and semiclassical resonances for n-disk repellers. 



Generalizing the work of Gaspard and Rice |13| to non-overlapping n-disk problems of 
arbitrary geometry and disk sizes we construct in Sec.^| the S-matrix from stationary scat- 
tering theory. Details of this calculation are relegated to App.[B|. Utilizing the machinery 
of trace-class operators which are summarized in App.|A] we construct the determinant of 
the n-disk S-matrix as the product of n one-disk determinants and the ratio of the deter- 
minant and its complex conjugate of the genuine multiscattering matrix. It is shown how 
the latter determinants split under symmetry operations. The proofs for the existence of 
the determinants are relegated to App.|C| and the comparison to alternative formulations 
of the multiscattering kernel can be found in App.|D]. 

In Sec.|] we state the link between the exact determinant of the n-disk S-matrix and the 
Gutzwiller-Voros curvature expansion. We discuss the semiclassical limit of the incoherent 
part, whereas the actual calculation is reported in App.||. It is shown that, under the 
semiclassical reduction of the quantum traces, the Plemelj-Smithies recursion relation 
for the quantum cumulants transforms into the recursion relation for the semiclassical 
curvatures which are known from the cycle expansion. 

The actual semiclassical reduction is worked out in Sec|| We start with the construction 
of the quantum cycles or itineraries built from the convolution of a finite number of multi- 
scattering kernels and show that they have the same symbolic dynamics in the full domain 
of an arbitrary n-disk system as their semiclassical counterparts, the geometrical periodic 
orbits. We discuss the case that the quantum-mechanical cumulant sum incorporates ge- 
ometries which classically allow for non-trivial pruning and hence for periodic orbits which 
pass undisturbed straight through a disk, see Refs.|S^,5S]. We show how these ghost or- 



bits cancel against their "parent" periodic orbits resulting from itineraries without the 
disk which is affected by the ghost passage. For the general case of an arbitrary quantum 



cycle Sec. 5.3 mirrors the semiclassical reduction of the convolution of two multiscattering 



kernels studied in detail in App.[F| with the help of the Watson contour integration and 



suitable deformations of the paths in the complex angular-momentum plane. In Sec.|5.4j the 



geometrical limit of a quantum cycle is studied, which is generalized to the case of r times 



repeated cycles in the following section. In Sec. |5.6| the ghost cancellation rule for arbitrary 
cycles is derived. The semiclassical diffractive creeping contributions are constructed and 
studied in Secs. |5.7| and |5.8| . Sec.|5] ends with the proof that an arbitrary quantum itinerary 
reduces semiclassically to a periodic orbit of Gutzwiller-Voros stability, such that the link 
between the exact determinant of the n-disk S-matrix and the Gutzwiller-Voros curvature 
sum is established. 

Numerical tests of the semiclassical curvature expansion can be found in SecM for the 
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example of the three-disk system in the Ai symmetry-class representation. First, the exact 
quantum-mechanical data are compared to the semiclassical predictions of the Gutzwiller- 
Voros zeta function, the dynamical zeta function ]36| and the quasiclassical zeta function 



suggested in f43fl , where all three semiclassical zeta function are expanded in curvatures 
which are truncated at finite order. Secondly, the exact cluster phase shift (defined by the 
phase of the determinant of the multiscattering matrix) is compared with the semiclassical 
predictions of the three zeta functions. Although all three zeta functions seem at first sight 
empirically equivalent, as they all predict the same leading resonances closest to the real 
fc-axis, this comparison shows clearly which of the three is superior and is hence the 
candidate for - at least - the FAPP ("for all practical purposes") zeta function. Sec.[] 
ends with an order-by-order comparison of the exact cumulants with their semiclassical 
counterpart, the curvatures for the Gutzwiller-Voros zeta function. From these numerical 
data we extract an empirical truncation rule for the curvature expansion as a function 
of the wave number. We relate this rule to the uncertainty bound resulting from finite 
quantum-mechanical resolution of the exponentially proliferating details of the classically 
repelling set. 

Sec.|7| concludes with a summary. Here we emphasis the preservation of unitarity under 
the semiclassical reduction, the decoupling of the incoherent one-disk from the coherent 
ra-disk determinants, and the particularities, when bounded domains are formed in the 
case of (nearly) touching disks. Furthermore, the resonance data are correlated with the 
truncation from the uncertainty bound. We discuss the relevance of those periodic orbits 
whose topological order exceeds the uncertainty bound. Arguments are presented that 
the Gutzwiller-Voros zeta function ought to be interpreted in the asymptotic sense as an 
truncated sum, whether it converges or not. The conclusions end with a discussion on h 
corrections. 



Note that the contents of Apps.jB] and |C| are based on M. Henseler's diploma thesis |46| 



while Sec. [3] as well as Secs|||3] and AppsJA].l-2 have partial overlap with Refs. [pTFH and 
48] , respectively. 
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2 Semiclassical resonances of the n-disk system 



The connection between exact quantum mechanics, on the one side, and semiclassics, 
on the other, for the n-disk repellers in the standard "bottom-up" approach, is rather 
indirect. It has been mainly based on a comparison of the exact and semiclassical pre- 
dictions for resonance data. In the exact quantum-mechanical calculations the resonance 
poles are extracted from the zeros of a characteristic scattering determinant (see ref. [0 
and below), whereas the semiclassical predictions follow from the zeros (poles) of one of 
the semiclassical zeta functions. These semiclassical quantities have either formally been 
taken over from bounded problems (where the semiclassical reduction is performed via 
the spectral density) [0,0] or they have just been extrapolated from the corresponding 
classical scattering determinant [f42| , |43| . Our aim is to construct a direct link between the 
quantum-mechanical and semiclassical treatment of hyperbolic scattering in a concrete 
context, the n-disk repellers. 



Following the methods of Gaspard and Rice [13] we will construct in Sec.w and App. 



the pertinent on-shell T-matrix which splits into the product of three matrices, namely 
C(fc)M _1 (/c)D(A;). The matrices C(k) and T)(k) couple the incoming and outgoing scat- 
tering wave (of wave number k), respectively, to one of the disks, whereas the matrix 
M(fc) parametrizes the scattering interior, i.e., the multiscattering evolution in the mul- 
tidisk geometry. The understanding is that the resonance poles of the n > 1 disk prob- 
lem can only result from the zeros of the characteristic determinant detM(&); see the 
quantum-mechanical construction of Gaspard and Rice jl3| for the three-disk scattering 
system fl0||IT] , |T4|| . Their work refers back to Berry's application [0,0] of the Korringa- 
Kohn-Rostoker (KKR) method ]51] to the (infinite) two-dimensional Sinai-billiard prob- 
lem which, in turn, is based on Lloyd's multiple scattering method p2] , [5^[ 
cluster of non-overlapping muffin-tin potentials in three dimensions. 



for a finite 



The resonance poles are calculated numerically by solving detM(fc) = in a finite, but 
large basis, such that the result is insensitive to an enlargement of the basis (see, e.g., 
44|). On the semiclassical side, the geometrical primitive periodic orbits (labelled by p) 
have been summed up - including repeats (labelled by r) - in the Gutzwiller-Voros zeta 
function |5|,|37 



Z GY (z;k) = expl 



i{z n n p (k)y 



V r=l ' L A; 




z n n p (ky 



A 3 

1 i.p 



(2.1) 



the dynamical zeta function of Ruelle |36 



1 „ 



C \z; k) = exp - 53 E ~J' nv W = II (1 - zHpt 



(2.2) 



L p r=l ) p 

(which is the j = part of the Gutzwiller-Voros zeta function) or the quasiclassical zeta 
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function of Vattay and Cvitanovic 



EST 



Zqcl(z; k) = exp i - H - 



p j=0 z=o 



(2.3) 



which is an entire function. In all cases t p (k) 



iv v TT/2 



/J\A P \ is the so-called p 



th 



cycle, rip is its topological length and z is a book-keeping variable for keeping track of 
the topological order. The input is purely geometrical, i.e., the lengths L p , the Maslov 
indices u p , and the stabilities (the leading eigenvalues of the stability matrices) A p of the 
pth p r j m itive periodic orbits. Note that both expressions for the three zeta functions, 
either the exponential one or the reformulation in terms of infinite product (s), are purely 
formal. In the physical region of interest, they may not even exist without regularization. 
(An exception is the non-chaotic 2-disk system, as it has only one periodic orbit, t Q (k).) 
Therefore, the semiclassical resonance poles are computed from these zeta functions in 



the curvature expansion |4^JT^JT^] up to a given topological length m. This procedure 
corresponds to a Taylor expansion of, e.g., Zqy(z; k) in z around z = up to order z m 
(with z set to unity in the end), e.g., 



z;k) = f-zT l — V (2.4) 




1 1 1 1 

n p =l n„/=l A„ A„/ 



+ 



The hope is that the limit m — > oo exists — at least in the semiclassical regime Re k ^> 1/a 
where a is the characteristic length of the scattering potential. We will show below that 
in the quantum-mechanical analog — the cumulant expansion - this limit can be taken, 
but that there are further complications in in the semiclassical case. 

The cycle expansion is one way of regularizing the formal expression of the Gutzwiller- 



Voros zeta function (|2.1). Another way would be the multiplication with a smooth cutoff 
function, as it is customary in quantum field theories, see e.g. 0. This is, in principle, 
allowed. In order to be able to compare quantum mechanics with semiclassics, however, the 
very same cutoff function has to be introduced already on the quantum level. Candidates 
for such cutoff functions which work on the quantum side and on the semiclassical side 
are not so obvious, see e.g. App. |A.4j . They have to be formulated in terms of the T-matrix 
or the multiscattering kernel and would introduce further complications. Fortunately, the 
quantum-mechanical side of the present problem exists without further regularization. 
Thus there is no need for an extra cutoff function. 

As mentioned above, the connection between quantum mechanics and semiclassics for 
these scattering problems has been the comparison of the corresponding resonance poles, 
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the zeros of the characteristic determinant on the one side and the zeros of the Gutzwiller- 
Voros zeta function or its competitors - in general in the curvature expansion - on the 
other side. In the literature (see, e.g., Refs. [fT2|JT7| , [T5|1 based on Ref. f54|| or f5"5| ) the link is 
motivated by the semiclassical limit of the left hand sides of the Krein-Friedel-Lloyd sum 
for the (integrated) spectral density |56|j5?H and |52|,|53[| 



lim lim (iV (n) (k + ie; b) - N {0) (k + ie; b)) = — ImTr In S(Jfc) , (2.5) 
lim lim (p (n) (A; + ie;6) - p^(k + ie;b)) = — ImTr^- hi S(A;) . (2.6) 



See also Ref. [|58[] for a modern discussion of the Krein-Friedel-Lloyd formula and Refs. [|55| , |5"9f 
for the connection of ( [2.6|) to the the Wigner time delay. In this way the scattering prob- 
lem is replaced by the difference of two bounded reference billiards (e.g. large circular 
domains) of the same radius b which finally will be taken to infinity, where the first con- 
tains the scattering region or potentials, whereas the other does not (see Fig. f2.lp . Here 
p( n \k; b) (N^ n \k; b)) and p^(k; b) (N^(k; b)) are the spectral densities (integrated spec- 
tral densities) in the presence or absence of the scatterers, respectively. In the semiclassical 
approximation, they will be replaced by a Weyl term and an oscillating sum over periodic 
orbits ||. 



• b \ b 



Fig. 2.1. The "difference" of two bounded reference systems, where one includes the scattering 
system. 

Note that this expression makes only sense for wave numbers above the real A;-axis. Es- 
pecially, if k is chosen to be real, e must be greater than zero. Otherwise, the exact left 
hand sides ( |2.5|) and ( |2.6| ) would give discontinuous staircase or even delta function sums, 
respectively, whereas the right hand sides are continuous to start with, since they can 
be expressed by continuous phase shifts. Thus, the order of the two limits in (p.5[) and 



J|) is essential, see e.g. Balian and Bloch [54| who stress that smoothed level densities 



should be inserted into the Friedel sums. In Ref.[T2[], chapter IV, Eqs. (4.1-4), the order 
is, however, erroneously inverted. 

Our point is that this link between semiclassics and quantum mechanics is still of very 
indirect nature, as the procedure seems to use the Gutzwiller-Voros zeta function for 
bounded systems and not for scattering systems and as it does not specify whether and 
which regularization has to be used for the semiclassical Gutzwiller trace formula. Neither 
the curvature regularization scheme nor other constraints on the periodic orbit sum follow 
naturally in this way. For instance, as the link is made with the help of bounded systems, 
the question might arise whether even in scattering systems the Gutzwiller-Voros zeta 
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function should be resummed a la Berry and Keating |2l| or not. This question is answered 
by the presence of the ie term and the second limit. The wave number is shifted from the 
real axis into the positive imaginary fc-plane. This corresponds to a "de-hermitezation" 
of the underlying exact hamiltonian - the Berry-Keating resummation should therefore 
not apply, as it is concerned with hermitean problems. The necessity of the +ie in the 
semiclassical calculation can be understood by purely phenomenological considerations: 
Without the ie term there is no reason why one should be able to neglect spurious periodic 
orbits which solely are there because of the introduction of the confining boundary. The 
subtraction of the second (empty) reference system helps just in the removal of those 
spurious periodic orbits which never encounter the scattering region. The ones that do 
would still survive the first limit b — > oo, if they were not damped out by the +ie term. 

Below, we will construct explicitly a direct link between the full quantum-mechanical 
S-matrix and the Gutzwiller-Voros zeta function. It will be shown that all steps in the 
quantum-mechanical description are well defined, as the T-matrix and the matrix A 
\I 1 are trace-class matrices (i.e., the sum of the diagonal matrix elements is absolutely 
converging in any orthonormal basis). Thus the corresponding determinants of the S- 
matrix and the characteristic matrix M are guaranteed to exist, although they are infinite 
matrices. 
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3 The n-disk S-matrix and its determinant 



Following the methods of Berry jKJ and Gaspard and Rice |0| we here describe the 
elastic scattering of a point particle from n hard disks in terms of stationary scattering 
theory. Because of the hard-core potential on the disk surfaces it turns into a boundary 
value problem. Let tp(r) be a solution of the pertinent stationary Schrodinger equation 
at spatial position r: 



(V 2 + k 2 ^j ip(r) = , r outside the n disks, 

ip{r) = , fon the surfaces of the disks, 

where E = h 2 k 2 /2m is the energy of the point-particle written in terms of its mass m 
and the wave vector k of the incident wave. After the wave function ip{r) is expanded in 
a basis of angular momentum eigenfunctions in two dimensions, it reads 



oo 
m=— oo 



where k and are the length and angle of the wave vector, respectively. The scattering 
problem in this basis reduces to 



(V 2 . + fc 2 ) ip m {r) = , f outside the disks; 

Vm(0 = , fon the disk surfaces. 

For large distances from the scatterers (kr — > oo) the spherical components i\) m can be 
written as a superposition of in-coming and out-going spherical waves, 



-| oo 

^(f ) ~ —L= Yl [<W e-^-i-'-f ) +S W eH^-f-'-D] e^'*- , (3.1) 
\ 2,71 kr m ' = _oo 

where r and $ r denote the distance and angle of the spatial vector f as measured in the 
global two-dimensional coordinate system. Eq.( |3.1|) defines the scattering matrix S which 
is unitary because of probability conservation. In the angular-momentum basis its matrix 
elements S mm ' describe the scattering of an in-coming wave with angular momentum m 
in an out-going wave with angular momentum m! . If there are no scatterers, then S = 1 
and the asymptotic expression of the plane wave e lfc r in two dimensions is recovered from 
i/j(r). All scattering effects are incorporated in the deviation of S from the unit matrix, 
i.e., in the T-matrix defined as S(k) = 1 — iT(k). In general, S is non-diagonal and 
therefore non-separable. An exception is the one-disk problem (see below). 

For any non-overlapping system of n disks (of, in general, different disk-radii aj, j = 
1, . . . , n) the S-matrix can be further split up. Using the methods and notation of Gaspard 
and Rice |TB[ this is achieved in the following way (see also Ref. |k| and App. [B] for a 
derivation of this result): 
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S mm ' (k) — 5 mm i iT mm , (k) 

= 5 mm , - iCJ(k)(M-\k))fnj; mf (k) . (3.2) 

Here the upper indices j, j' — 1, . . . , n < oo label the n different disks, whereas the lower 
indices are the angular momentum quantum numbers. Repeated indices are summed over. 
The matrices C- 7 and D- 7 depend on the origin and orientation of the global coordinate 
system of the two-dimensional plane and are separable in the disk index j: 



3 1 



iy{ m ? = -TrafJ m >- l '(kR jl )J v (ka j >)e im#i V , (3.4) 

where Rj and are the distance and angle, respectively, of the ray from the origin 
in the 2-dimensional plane to the center of disk j as measured in the global coordinate 
system (see Fig. |3.1| ). H\ l \kr) is the ordinary Hankel function of first kind and Ji(kr) 
the corresponding ordinary Bessel function. The matrices C- 7 and D J parameterize the 
coupling of the incoming and outgoing scattering waves, respectively, to the scattering 
interior at the j th disk. Thus they describe only the single-disk aspects of the scattering 
of a point particle from the n disks. The matrix M- 7 - 7 has the structure of a Kohn-Korringa- 
Rostoker (KKR)-matrix, see Refs. |4^j50|j5^| , 



M# = 5 jf 5 u > + (1 - 5 jf ) ^ -4^r H^(kR jf ) T jf (l, I') . (3.5) 

a j' Hi, '(kdj 



without Ewald resummation [p^fl , as the number of disks is finite. Here Rjy is the separa- 
tion between the centers of the jth and j'th disk and Rjj> = Rj'j, of course. The auxiliary 
matrix Tjji(l,V) = e^ la i'j~ l ( a jj' _7r )) contains - aside from a phase factor - the angle ctj/j 
of the ray from the center of disk j to the center of disk j' as measured in the local 
(body-fixed) coordinate system of disk j (see Fig.|3.1|). 



Note that T jy (/,/') = Ty^l', I) )*. The "Gaspard and Rice prefactors" of M, 

i.e., (ira/2i) in ||13|| , are rescaled into C and D. The matrix A = M — 1 contains the 
genuine multidisk "scattering" aspects of the the n-disk problem, e.g., in the pure 1-disk 
scattering case, A vanishes. When (M -1 )-" is expanded as a geometrical series about the 
unit matrix {5 JJ }, a multiscattering series in "powers" of the matrix A is created. 

The product CM _1 D is the on-shell T-matrix of the n-disk system. It it the two- 
dimensional analog of the three-dimensional result of Lloyd and Smith for a finite cluster 
of non-overlapping muffin-tin potentials. At first sight the expressions of Lloyd and Smith 
(see Eq.(98) of [|53[] and also Berry's form [49| for the infinite Sinai cluster) seem to look 



simpler than ours and the original ones of Ref. ||13||, as, e.g., in M the asymmetric term 
ajJi(kaj)/ajiH\}\kaji) is replaced by a symmetric combination, Ji(kaj) / H^lkaj) . Under 
a formal manipulation of our matrices we can derive the same result (see App. 0). In fact, 
it can be checked that the (formal) cumulant expansion of Lloyd's and our M-matrix are 
identical and that also numerically the determinants give the same result. Note, however, 
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Fig. 3.1. Global and local coordinates for a general 3-disk problem. 



that in Lloyd's case the trace-class property of M is lost, such that the infinite deter- 
minant and the corresponding cumulant expansion converge only conditionally, and not 
absolutely as in our case. The latter fact is based on the trace-class properties of the 
underlying matrices and is an essential precondition for all further simplifications, as e.g. 
unitary transformations, diagonalization of the matrices, etc. 

A matrix is called "trace-class" , if and only if, for any choice of the orthonormal basis, the 
sum of the diagonal matrix elements converges absolutely; it is called "Hilbert-Schmidt", 
if the sum of the absolute squared diagonal matrix elements converges (see M. Reed and 
B. Simon, Vol.1 and 4 |60| , |6l| and App. [A] for the definitions and properties of trace-class 
and Hilbert-Schmidt matrices). Here, we will only list the most important properties: 

(i) Any trace-class matrix can be represented as the product of two Hilbert-Schmidt 
matrices and any such product is again trace-class. 

(ii) A matrix B is already Hilbert-Schmidt, if the trace of B^B is absolutely convergent 
in just one orthonormal basis. 

(iii) The linear combination of a finite number of trace-class matrices is again trace-class. 

(iv) The hermitean-conjugate of a trace-class matrix is again trace-class. 

(v) The product of two Hilbert-Schmidt matrices or of a trace-class and a bounded 
matrix is trace-class and commutes under the trace. 

(vi) If a matrix B is trace-class, the trace tr(B) is finite and independent of the basis. 

(vii) If B is trace-class, the determinant det(l + zB) exists and is an entire function of z. 

(viii) If B is trace-class, the determinant det(l + zH) is invariant under any unitary trans- 
formation. 

In App. |C| we show explicitly that the /-labelled matrices S^ n \k) — 1, C^(k) and D J (A;) 
as well as the {/, j}-labelled matrix A(k) = lSA(k) — 1 are of "trace-class", except at the 
countable isolated zeros of H^'(kaj) and of DetM(/c) and at k < 0, the branch cut of the 
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Hankel functions. The ordinary Hankel functions have a branch cut at negative real k, such 
that even the /c-plane is two-sheeted. The last property is special for even dimensions and 
does not hold in the 3-dimensional n-ball system |46|J62]] . Therefore for almost all values of 
the wave number k (with the above mentioned exceptions) the determinant of the n-disk 
S-matrix exist and the operations of (3J3) are mathematically well defined. We concentrate 
on the determinant, detS, of the scattering matrix, since we are only interested in spectral 
properties of the n-disk scattering problem, i.e. resonances and phase shifts, and not in 
wave functions. Furthermore, the determinant is invariant under any change of a complete 
basis expanding the S-matrix and therefore also independent of the coordinate system. 



det*S (n) = deti (l - iCM _1 D) = exp tr { In (l - iCM _1 D 



: exp 



oo -JV 

£ N tTl 



oo • jV 



\ N=l 

(oo 
~ ^ N Tl ' 
N=l iV 

:expTr L ln(l-iM -1 DC 

Det L (M _1 (M-iDC 
Det L (M-iDC) 
Det L (M) ' 



(CM _1 D 
(M _1 DC 



N 



Detr 1- 



iM _1 DC) 



(3.6) 



We use here exp tr In notation as a compact abbreviation for the defining cumulant expan- 
sion ([A. 71) , since det(l + /iA) = exp(— J2n=i ^~n~ t r (A N )), is only valid for |/i|max|Aj| < 1 



where Aj is the i-th eigenvalue of A. The determinant is directly defined by its cumulant 
expansion (see Eq.(188) of Ref . [plj] and Eq.( |A.7| ) of App. 0.2) which is therefore the 
analytical continuation of the e trlog -representation. 

The capital index L is a multi- or "super"-index L = On the l.h.s. of Eq. ( |3.6| ) the 

determinant and traces are only taken over small I, on the r.h.s. they are taken over the 
super-indices L = (/, j). In order to signal this difference we will use the following notation: 
det . . . and tr . . . refer to the \m) space, Det . . . and Tr . . . refer to the super-space. The 
matrices in the super-space are expanded in the complete basis {\L)} = {\m;j)} which 
refers for fixed index j to the origin of the jth disk and not any longer to the origin of 
the 2-dimensional plane. In deriving (|3.6|) the following facts were used: 



(a) D J , G 7 are of trace-class in the {|/)} space (see App.|C|). 

(b) As long as the number of disks is finite, the product DC - now evaluated in the 
super-space {\L)} - is of trace-class as well (see property (iii)). 

(c) M — 1 is of trace-class (see App. y). Thus the determinant DetM(fc) exists. 

(d) Furthermore, M is bounded (since it is the sum of a bounded and a trace-class 
matrix) . 

(e) M is invertible everywhere where Det^M(/c) is defined (which excludes a countable 
number of zeros of the Hankel functions H^\kaj) and the negative real k there 
is a branch cut) and nonzero (which excludes a countable number of isolated points in 
the lower /c-plane) - see property (e) of App. [A].2. Therefore and because of (d) the 
matrix M _1 is bounded. 
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(f) The matrices CM _1 D, M _1 DC, are all of trace-class as they are the product of 
bounded times trace-class matrices and tr m [(CM _1 D) 7V ] = TrM[(M _1 DC) Ar ], because 
such products have the cyclic permutation property under the trace (see properties (hi) 
and (v)). 

(g) M — iDC — 1 is of trace-class because of the rule that the sum of two trace-class 
matrices is again trace-class (see property (hi)). 

Thus all traces and determinants appearing in Eq. (|3.6| ) are well-defined, except at the 
above mentioned isolated k singularities and branch cuts. In the {|m; j)} basis the trace 
of M — 1 vanishes trivially because of the 5jj> terms in (|3.5| ). One should not infer from 
this that the trace-class property of M — 1 is established by this fact, since the finiteness 
(here vanishing) of Tr(M — 1) has to be shown for every complete orthonormal basis. After 
symmetry reduction (see below) Tr(M — 1), calculated for each irreducible representation 
separately, does not vanish any longer. However, the sum of the traces of all irreducible 
representations weighted with their pertinent degeneracies still vanishes of course. Semi- 
classically, this corresponds to the fact that only in the fundamental domain there can 
exist one-letter "symbolic words" . 

After these manipulations, the computation of the determinant of the S-matrix is very 
much simplified in comparison to the original formulation, since the last term of Eq. ( |3.6| ) 
is completely written in terms of closed form expressions and since the matrix M does 
not have to be inverted any longer. Furthermore, as shown in App. [FJ.3, one can easily 
construct 



HP(ka f 



M» -iDL,C m ,£ =5^ i — — (iy 



Hf } (kaj> 

M - iff J ° - H ^ kR ^ r »'<!> > ( 3 - 7 ) 



where H^(kr) is the Hankel function of second kind. The first term on the r.h.s is just 
the S-matrix for the separable scattering problem from a single disk, if the origin of the 
coordinate system is at the center of the disk (see App. 0.2): 



H^\ka 
Hp\ka 



After fl3.7|) is inserted into (|3.6|) and Q3.8|) is factorized out, the r.h.s. of ( |3.6|) can be 
rewritten as 

= De tL[ M W -.DWC(k)l = , A Det,M(r)> 

v ' Det L M(k) K JJ J j Det L M(A;) v ; 

where {H%\z)}* = H${z*) has been used in the end. All these operations are allowed, 
since M(/c) — 1, M(/c) — iD(A;)C(A;) — 1 and S^(k) — 1 are trace-class for almost every 
k with the above mentioned exceptions. In addition, the zeros of the Hankel functions 
Hffi(kaj) now have to be excluded as well. In general, the single disks have different sizes 
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and the corresponding 1-disk S- matrices should be distinguished by the index j. At the 
level of the determinants this labelling is taken care of by the choice of the argument 
Note that the analogous formula for the three-dimensional scattering of a point particle 
from n non-overlapping spheres (of in general different sizes) is structurally completely 
the same [i6,62], except that there is no need to exclude the negative fc-axis any longer, 



since the spherical Hankel functions do not posses a branch cut. In the above calculation 
it was used that T*-,{l,V) = Tjji(—l, —I') in general |46j and that for symmetric sys- 



tems (equilateral 3-disk-system with identical disks, 2-disk system with identical disks): 
r*^/, /') = I') (see [|13j). Eq.( |3.9| ) is compatible with Lloyd's formal separation of the 

single scattering properties from the multiple-scattering effects in the Krein-Friedel-Lloyd 
sum, see, e.g., p. 102 of Ref.[|53| (modulo the above-mentioned conditional convergence 
problems of the Lloyd formulation). Eq. (|3.9|) has the following properties: 

(i) Under the determinant of the n-disk -matrix, the 1-disk aspects separate from the 
multiscattering aspects, since the determinants of the 1-disk S*- 1 - 1 matrices factorize from 
the determinants of the multiscattering matrices. Thus the product over the n 1-disk 
determinants in ( |3.9| ) parametrizes the incoherent part of the scattering, as if the n-disk 
problem just consisted of n separate single-disk problems. 

[ii) The whole expression (|3.9|) respects unitarity as is unitary by itself, because of 



(Hffl(z))* = H$(z*) and as the quotient of the determinants of the multiscattering ma- 
trices on the r.h.s. of ( |3.9| ) is manifestly unitary. 

(iii) The determinants over the multiscattering matrices run over the super-index L of 
the super-space. This is the proper form for the symmetry reduction (in the super-space), 
e.g., for the equilateral 3-disk system (with disks of the same size) we have 



Det L M 3 - dis k = det iAi M Al det^M^ (det /E M E ) 2 , (3.10) 
and for the 2-disk system (with disks of the same size) 

Det L M 2 - disk = det, Ai M Al det^M^ det* Bi M Bl deti B2 M B2 , (3.11) 

etc. In general, if the disk configuration is characterized by a finite point-symmetry group 
Q, we have 



Det L M n _ disk = [] (det /c M Dc (A;)) dc , (3.12) 

c 

where the index c runs over all conjugate classes of the symmetry group Q and D c is the 
c th representation of dimension d c For the symmetric 2-disk system, these represen- 
tations are the totally symmetric A 1; the totally anti-symmetric A 2 , and the two mixed 
representations Bx and B 2 which are all one-dimensional. For the symmetric equitriangular 
3-disk system, there exist two one-dimensional representations (the totally symmetric Ai 
and the totally anti-symmetric A 2 ) and one two-dimensional representation labelled by 
E. A simple check that DetM(A;) has been split up correctly is the following: the power 
of H^(kaj) Hankel functions (for fixed m with — oo < m < +oo) in the denomina- 
tor of He [det/ c MD c (fc)] dc has to agree with the power of the same functions in DetM(A;) 
which in turn has to be the same as in Ilj=i (detS^/ca,-)). Note that on the l.h.s. the 
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determinants are calculated in the super-space {L}, whereas on the r.h.s. the reduced 
determinants are calculated, if none of the disks are special in size and position, in the 
normal (desymmetrized) space {/} (however, now with respect to the origin of the disk 
in the fundamental domain and with ranges given by the corresponding irreducible repre- 
sentations). If the n-disk system has a point-symmetry where still some disks are special 
in size or position (e.g., three equal disks in a row the determinants on the r.h.s. 



refer to a correspondingly symmetry-reduced super-space. This summarizes the symmetry 
reduction on the exact quantum-mechanical level. It can be derived from 



Det L M = exp 



: exp 



■ exp 



: exp 




(3.13) 



where U is unitary transformation which makes A block-diagonal in a suitable trans- 
formed basis of the original complete set {|m; j)}. These operations are allowed because 
of the trace-class-property of A and the boundedness of the unitary matrix U (see also 
property (d) of App. 0.2). 
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4 The link between the determinant of the S-matrix and the semiclassical 
zeta function 



In this chapter we will specify the semiclassical equivalent of the determinant of the 
n-disk S-matrix. As detS^ in (|3.9| ) factorizes into a product of the 1-disk determinants 
and the ratio of the determinants of the multiscattering matrix, DetM(fc*)'/DetM(A;), the 
semiclassical reduction will factorize as well into incoherent one-disk parts and an coherent 
multiscattering part. Note, however, that there is an implicit connection between these 
parts via the removable one-disk poles and zeros. This will be discussed in the conclusion 
section |7]. 

In App. [E|, the semiclassical expression for the determinant of the 1-disk S-matrix is 



constructed in analogous fashion to the semiclassical constructions of Ref.[44| which in 



turn is based on the work of Ref.[29l: 



detjS (1) (A;a) « (e" 1 



(4.1) 



{UZi t 1 -e+iWM])" 

with the creeping exponential (for more details, see App. [E] and the definitions of App. |F.4|) 



v t {ka) = ka + e +i ^ 3 (ka/6) 1/3 q e + • • • = ka + ia e {ka) + ■■■ , (4.2) 
v t {ka) = ka + e~ i7r/3 (ka/6) 1/3 q e + ■ ■ ■ = ka - i{a e {k*a))* + ■■■ 

= (Mk*a))* , (4.3) 

and N{ka) = (7ra 2 k 2 )/47r+- ■ ■ the leading term in the Weyl approximation for the staircase 
function of the wave-number eigenvalues in the disk interior. From the point of view of 
the scattering particle the interior domains of the disks are excluded relatively to the free 
evolution without scattering obstacles (see, e.g., Therefore the negative sign in front 
of the Weyl term. For the same reason, the subleading boundary term has here a Neumann 
structure, although the disks have Dirichlet boundary conditions. Lets us abbreviate the 
r.h.s. of ( |4.1| ) for a specified disk j as 



det,S«(* ai ) ^ (e-^)) 2 Z i-**0)(fc'%)' Z^^a£ ^ ^ 

where ^i-di s k(i) (kcij) and Zi-disk(r) (^ a j) are the diffractional zeta functions (here and in the 
following we will label semiclassical zeta-functions with diffractive corrections by a tilde) 
for creeping orbits around the j'th disk in the left-handed sense and the right-handed sense, 
respectively (see Fig. |4.1|) . The two orientations of the creeping orbits are the reason for 
the exponents two in ( |4.1| ). Eq.( f4.1| ) describes the semiclassical approximation to the 
incoherent part (= the curly bracket on the r.h.s.) of the exact expression ( p.9| ). 

We now turn to the semiclassical approximation of the coherent part of ( |3.9| ), namely 
the ratio of the determinants of the multiscattering matrix M. Because of the trace-class 



24 



Fig. 4.1. Left- and right-handed diffractive creeping paths of increasing mode number £ for a 
single disk. 

property of A = M — 1, the determinants in the numerator and denominator of this 
ratio exist individually and their semiclassical approximations can be studied separately. 
In fact, because of DetM(A;*) t = (DetM(/c*))*, the semiclassical reduction of DetM(fc*) t 
follows directly from the corresponding result of DetM(ft) under complex conjugation. 
The semiclassical reduction of DetM(/c) will be done in the cumulant expansion, since 
the latter is the defining prescription for the computation of an infinite matrix that is of 
the form 1 + A where A is trace-class: 



z 2 

Det [1 + zA(k)} = 1 - (-z)Tr[A(k)} - — {Tr [A 2 (£;)] - [TrA(/c)] 2 | + • • • 

oo 

= E^JA) withQ„(A) = l (4.5) 

n c =0 

where we have introduced here a book-keeping variable z which we will finally set to one. 
This allows us to express the determinant of the multiscattering matrix solely by the 
traces of the matrix A, Tr[A m (/c)] with m = 1, 2, 3, • • •. The cumulants and traces satisfy 
the (Plemelj-Smithies) recursion relations (|A.16|) 



Q nc (A) = -J2 (-l) m+1 Qn c - m (A)Tr [A m ] for n c > 1 (4.6) 



n c m=l 



in terms of the traces. In the next section we will utilize Watson resummation tech- 
niques p4] , p9| which help to replace the angular momentum sums of the traces by contin- 



uous integrals which, in turn, allow for semiclassical saddle-point approximations. With 
these techniques and under complete induction we will show that for any geometry of n 
disks, as long as the number of disks is finite, the disks do not overlap and grazing or 
penumbra situations |3~1"|,|3"5| are excluded (in order to guarantee unique isolated saddles), 



the semiclassical reduction reads as follows: 



Tr [A m (A0] ^ (~l) m E E S m,m p n p ^\ y + creeping p.o.'s (4.7) 

1 ~ W 



P r>0 



with inputs as defined below (|2.3|) . The reduction is of course only valid, if Rek is suffi- 
ciently large compared to the inverse of the smallest length scale of the problem. The right 
hand side of Eq. ( |4.7| ) can be inserted into the recursion relation ( |4.6| ) which then reduces 
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to a recursion relation for the semiclassical approximations of the quantum cumulants 

C nc (s.C.) = J2 Cn c -m(s.C.) S m,rn p n p P . - y for Tl c > 1 (4.8) 

U c m =l p r>0 1 - (x;J 

where we have neglected the creeping orbits for the time being. Under the assumption 
that the semiclassical limit Re k — > oo and the cumulant limit n c — > oo commute (which 
might be problematic as we will discuss later), the approximate cumulants C nc (s.c.) can 
be summed to infinity, Y^ c =o znc C nc (s.c), in analogy to the exact cumulant sum. The 
latter exists since A is trace-class. The infinite "approximate cumulant sum" , however, is 
nothing but the curvature expansion of the Gutzwiller-Voros zeta function, i.e., 

oo 

^Gv(-2;^)|curv.rcg. = (s.C. ) , (4.9) 

n c =0 

since Eq. Q4.8|) is exactly the recursion relation of the semiclassical curvature terms || . 

If, in addition, the creeping periodic orbits are summed as well, the standard Gutzwiller- 
Voros zeta function generalizes to the diffractive one discussed in Refs. p2H54]j which we 
will denote here by a tilde. In summary, we have 



DetM(fc) ^> Z GV (£0Uv.re g . (4.10) 
for a general geometry and 

detM Dc (A;) ^ Z Dc (A;)| curv . reg . (4.11) 

for the case that there is a finite point-symmetry and the determinant of the multiscatter- 
ing matrix splits into the product of determinants of matrices belonging to the pertinent 
representations D c , see Eq.( |3.12| ). Thus the semiclassical limit of the r.h.s. of Eq. (|3^) is 



3=1 

■mN(k aj )\ 2 Zi- disk (i)(k*aj) Zi- Aisk{v )(k*aj) [ Z GY (k*) 



n(«- 



Z\-d\sk(\){kaj) Zi- disk ( r )(kaj) J Z G y(k) 

(4.12) 



where, from now on, we will suppress the qualifier • • • | C urv.rc g .- For systems which allow 
for complete symmetry reductions (i.e., equivalent disks under a finite point-symmetry 
with ctj = aWj) the link reads 



det,S (n) (£;) = (det ; S (1) (A;a 



n c (det ic M Dc (/c)) f 
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(< 



-i7r7V(fca)^| 2n [ ^l-disk(l)(fc*g) ^l-disk(r) (jgg) | x 
^l-disk(l)(^«) ^l-disk(r) (kd) 



n 



x3^t (4,3) 

in obvious correspondence. Note that the symmetry reduction from the right hand side 
of Q4.12] ) to the right hand side of ( |4.13|) is compatible with the semiclassical results of 
Refs. [l65|J& 



In the next section we will prove the semiclassical reduction step (|4.7|) for any n-disk 
scattering system under the conditions that the number of disks is finite, the disks do not 
overlap, and geometries with grazing periodic orbits are excluded. We will also derive the 
general expression for creeping periodic orbits for n-disk repellers from exact quantum 
mechanics and show that ghost orbits drop out of the expansion of Trln(l + A) and 
therefore out of the cumulant expansion. 
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5 Semiclassical approximation and periodic orbits 



In this section we will work out the semiclassical reduction of Tr [ A m (k) ] for non-overlapping, 
finite ra-disk systems where 

Mi 1 ' = (1 " M a M^\ ^' H i-i>(kRjj') ■ (5.1) 

As usual, cij, ay are the radii of disk j and j', 1 < j, f < n, Rjj> is the distance between the 
centers of these disks, and aj>j is the angle of the ray from the origin of disk j to the one 
of disk j' as measured in the local coordinate system of disk j. The angular momentum 
quantum numbers I and I' can be interpreted geometrically in terms of the positive- or 
negative- valued distances (impact parameters) l/k and I'/k from the center of disk j and 
disk j', respectively, see f49fl. 



Because of the finite set of n disk-labels and because of the cyclic nature of the trace, 
the object Tr[ A m (k) } contains all periodic itineraries of total symbol length m with an 
alphabet of n symbols, i.e. A jlj2 A j2js ■ ■ ■ A jm - ljm A jmjl with ji G {1, 2, ... , n}. Here the 
disk indices are not summed over and the angular momentum quantum numbers are 
suppressed for simplicity. The delta-function part (1 — Sjf) generates the trivial pruning 
rule (valid for the full n-disk domain) that successive symbols have to be different. We 
will show that these periodic itineraries correspond in the semiclassical limit, kay ^> 1, 
to geometrical periodic orbits with the same symbolic dynamics. For periodic orbits with 



creeping sections P^ , |45| , |32T -p4| the symbolic alphabet has to be extended. Furthermore, 
depending on the geometry, there might be non-trivial pruning rules based on the so-called 
ghost orbits, see Refs.R,E9|. We will discuss such cases in Sec.|5\2 . 



5.1 Quantum itineraries 



As mentioned, the quantum-mechanical trace can be structured by a simple symbolic 
dynamics, where the sole (trivial) pruning rule is automatically taken care of by the 1 — Sjj> 
factor appearing in A.jf t . Thus we only have to consider the semiclassical approximation 
of a quantum-mechanical itinerary of length m: 



+00 +00 +OG +00 +00 

■■= Y Y Y ■■■ Y Y AfZ-'A/r ••• Af- 7" a/7 : (5.2) 

— OO ^2 — — ^3 — — OO l 7n — \= — OO / m = — OO 

with ji G {1, 2, ... , n}. This is still a trace in the angular momentum space, but not any 
longer with respect to the superspace. Since the trace, TrA m , itself is simply the sum of 
all itineraries of length m, i.e. 



n n n 71 n 

Tr A m = E E E • • ■ E E A A^- ; • ■ • A 3m - ljm A jmn , (5.3) 

.71 = 1.72=1 .73 = 1 im-l=ljm = l 
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its semiclassical approximation follows directly from the semiclassical approximation of 
its itineraries. Note that we here distinguish between a given itinerary and its cyclic 
permutation. All of them give the same result, such that their contributions can finally 
be summed up by an integer- valued factor n p := m/r, where the integer r counts the 
number of repeated periodic subitineraries. Because of the pruning rule 1 — 8jj/, we only 
have to consider traces and itineraries with n > 2 as A.\f, = implies that Tr A = in the 
full domain. 



We will show in this section that, with the help of the Watson method |6^J29|1 (studied 



for the convolution of two A matrices in App.[F] which should be consulted for details), 
the semiclassical approximation of the periodic itinerary 



A' 



JU2 A 



A 



becomes a standard periodic orbit labelled by the symbol sequence jij'2 • • ■ jm- Depending 
on the geometry, the individual legs — > ji — > ji + \ result either from a standard 
specular reflection at disk ji or from a ghost path passing straight through disk ji. If 
furthermore creeping contributions are taken into account, the symbolic dynamics has to 
be generalized from single-letter symbols {ji} to triple-letter symbols {ji,Si x with 
li > 1 integer- valued and Sj = 0, ±1 ^3 By definition, the value = represents the 
non-creeping case, such that {ji,0 x i { } = {ji,0} = {ji} reduces to the old single-letter 
symbol. The magnitude of a non-zero li corresponds to creeping sections of mode number 
\£i\, whereas the sign Sj = ±1 signals whether the creeping path turns around the disk ji 
in the positive or negative sense. Additional full creeping turns around a disk j' can be 
summed up as a geometrical series; therefore they do not lead to the introduction of a 
further symbol. 



5.2 Ghost contributions 




Fig. 5.1. The ghost itinerary (1,2,3,4). 

An itinerary with a semiclassical ghost section at, say, disk ji will be shown to have the 
same weight as the corresponding itinerary without the ji th symbol. Thus, semiclassically, 
they cancel each other in the Trln(l + A) expansion, where they are multiplied by the 
permutation factor m/r with the integer r counting the repeats. E.g. let (1,2,3,4) be a 

# 1 Actually, these are double-letter symbols as Sj and k are only counted as a product. 
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non-repeated periodic itinerary with a ghost section at disk 2 steming from the 4th-order 
trace TrA 4 , where the convention is introduced that an underlined disk index signals a 
ghost passage (see Fig. [5.1| ). Then its semiclassical, geometrical contribution to Tr ln(l+A) 
cancels exactly against the one of its "parent" itinerary (1,3,4) (see Fig. |5.2j ) resulting from 
the 3rd-order trace: 




Fig. 5.2. The parent itinerary (1,3,4). 



_i (a a 1 '- A 1,3 A 3 - 4 A 4 ' 1 U-fu 1 ' 3 A 3 ' 4 A 4 ' 1 

^ \^ ghost ghost geom geom J g \ geom geom geom j 
1,3 A 3,4 A 4,l n 

(5.4) 



.1+1) AM A 3,4 A 4 

1 / geom geom geom 



The prefactors +1/3 and —1/4 are due to the expansion of the logarithm, the factors 3 
and 4 inside the brackets result from the cyclic permutation of the periodic itineraries, 
and the cancellation stems from the rule 



. . . aM±1 A i±l,i+2 = Ki,i+2 . . _ /k c\ 

ghost ghost geom " \ U ' U J 

We have checked this rule in App. [F.6| for the convolution of two A-matrices, but in Sec. |5.6| 
we will prove it to hold also inside an arbitrary (periodic) itinerary. Of course the same 
cancellation holds in case that there are two and more ghost segments. For instance, 
consider the itinerary (1, 2, 3, 4, 5, 6) with ghost sections at disk 2 and 5 resulting from the 
sixth order trace. Its geometrical contribution cancels in the trace-log expansion against 
the geometrical reduction of the itineraries (1,2, 3, 4, , 6), (1, 3, 4, 5, 6) from the 5th-order 
trace with ghost sections at disk 2 or 5, respectively, and against the geometrical reduction 
of the itinerary (1, 3, 4, 6) of the 4th-order trace with no ghost contribution: 



--(6 A 1 ' 1 A-' 3 A 3 ' 4 A 4 '- A-' 6 A 6 - 1 

g ^ ghost ghost geom ghost ghost geom 

-lI^a 1 '- a-' 3 a 3 ' 4 a 4 ' 6 a 6 ' 1 -^a 1 ' 3 a 3 ' 4 a 4 '- a-' 6 a 6 ' 1 "l r^fii 

' g ^ ghost ghost geom geom geom ' " geom geom ghost ghost geomy V"" u / 

-T4A 1 ' 3 A 3 ' 4 A 4,6 A 6 - 1 

a V geom geom geom geoi 



^ V geom — geom — geom - geom 

4-9-lU 

geom geom geom geom 



1+2 -11 A 1,3 A 3 ' 4 A 4,6 A 6,1 =D (<, 7) 

J geom geom geom geom w • V ' / 



30 



Again, the prefactors —1/4, +1/5, —1/6 result from the trace-log expansion, the factors 
4, 5, 6 inside the brackets are due to the cyclic permutations, and the rule (|5.5| ) was used. 
If there are two or more ghost segments adjacent to each other, the ghost rule ( |5.5| ) has 
to be generalized to 



» i, i+l » i+l , i+2 . i+k , i+k+l . i+n— l ,i+n 

ghost ghost ghosts ghosts 



» i.i+2 
• ■ A — 

ghost 
. i, i+3 
^ ghost 

A M+ n 

geom 



••A' 
•■A' 



i +k , i+k+l 

hosts 
i+k , i+k+l 
ghosts 



■ A- 
A 



i+n— l,i+n 
ghosts 
i+n— l,i+n 



ghosts 



(5.8) 
(5.9) 
(5.10) 



Finally, let us discuss one case with a repeat, e.g. the itinerary (1,2,3,4,1,2,3,4) with 
repeated ghost sections at disk 2 in the semiclassical limit. The cancellations proceed in 
the trace-log expansion as follows: 



-s( 3A i 



A 2 -' 3 A 3 ' 4 A 4 ' 1 A 1 '- A-' 3 A 3 ' 4 A 

ghost ghost geom geom ghost ghost geom geom 



. 4,1 



2 

ghost 



a2,3 A 3,4 a4,1 a 1,3 *3,4 a4,1 



A 3,4 a 1,1 \L,a a 

geom geom geom geom geom geom 



geom 
4,1 



geom geom geom gc 



1.3 



.3,4 



4,1 



-i + l-I 

2 2 



A 1 ' 3 A 3 ' 4 A' 

geom geom geom 



.4,1 



(5.11) 



Note that the cyclic permutation factors of the 8th and 6th order trace are halved because 
of the repeat. The occurrence of the ghost segment in the second part of the 7th order 
itinerary is taken care of by the weight factor 7. 



The reader might study more complicated examples and convince him- or herself that 
the rule ( [5 . 1 Q| ) is sufficient to cancel any primary or repeated periodic orbit with one or 
more ghost sections completely out of the expansion of Tr ln(l + A) and therefore also 
out of the cumulant expansion in the semiclassical limit: Any periodic orbit of length 
m with n(< m) ghost sections is cancelled by the sum of all 'parent' periodic orbits of 
length m — i (with 1 < i < n and i ghost sections removed) weighted by their cyclic 
permutation factor and by the prefactor resulting from the trace-log expansion. This is 
the way in which the non-trivial pruning for the ra-disk billiards can be derived from the 
exact quantum-mechanical expressions in the semiclassical limit. Note that there must 
exist at least one index i in any given periodic itinerary which corresponds to a non-ghost 
section, since otherwise the itinerary in the semiclassical limit could only be straight and 
therefore non-periodic. Furthermore, the series in the ghost cancelation has to stop at 
the 2nd-order trace, TrA 2 , as TrA itself vanishes identically in the full domain which is 
considered here. 
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5. 3 Semiclassical approximation of a periodic itinerary 



The procedure for the semiclassical approximation of a general periodic itinerary, Eq. fl5.2p , 
of length m follows exactly the calculation of App.[F] for the convolution of two A-matrices. 
The reader interested in the details of the semiclassical reduction is advised to consult 
this appendix before proceeding with the remainder of the section. First, for any index 
i, 1 < i < m, the sum over the integer angular momenta, will be symmetrized as in 
Eq.flPjjjD with the help of the weight function d(k) [d(k^0) = 1, d(k=0) = 1/2]. 



A nj ' 2 ■ ■ ■ A Ji ~ lM A M:,i+1 • • • A- 7 " 1 - 71 

+oo +oo +oo +oo +oo +oo 

= y y ... y y y ... y Af/ 2 ---Af-f A/f +1 •••Af-/ 1 

l\=— 00/2— ~oo li—\=— OO li=— OO li+l=— 00 l 7n =—oo 
+00 +00 +00 +00 +00 +00 

= E E E E ••• E E ^-^-iKy-O 

Zi=0 l 2 =0 k-i=0 k=0 h+i=0 l m =0 Bi,...,a m =-l,l 

x A 3132 . . . A ii-lii A . . . A imjl 

(5.12) 



Furthermore, the angles A«j 4 = aj i+1 ^ — otj^^ [the analogs of a^y — a^y in Eq. ([F.3|) 1 
will be replaced by Aay jC ^ = Aa Jt — <7j27r where Oi = 0, 2, 1. This will be balanced by 
multiplying Eq. (|5.2|) with (— where a[ = o~i for o~i = 1 and zero otherwise. The three 
choices for Oi are, at this stage, equivalent, but correspond in the semiclassical reduction 
to the three geometrical alternatives: specular reflection at disk ji to the right, to the 
left or ghost tunneling. In order not to be bothered by borderline cases between specular 
reflections and ghost tunneling, we exclude disk configurations which allow classically 
grazing or penumbra periodic orbits [0,0 • 

Then, the sum over the integer angular momentum Zj will be replaced by a Watson contour 
integration over the complex angular momentum v j 

+00 1 1 

y (-iy^d(k)X k = -Uv i —— e~^< X Vi , (5.13) 

f n 2l J SUl ZA-7T 

as in Eq. (F.4|) . The quantity abbreviates here 



x JjJMiiL H {i) , kR (/,•/,' ; . )<^ /A "^ 

H U ( ka 3i) «i=-l,l 



-Y h (5.14) 



where the expression has simplified because of ( fca^ ) / -f^ 1 / (^a^ ) = J/ i ( fca^ ) / Hjy (kaj. ) , 
since is an integer. The quantity Y\ i abbreviates the sum in (|5.14|) . The next steps are 
completely the same as in App. |F.l| - |F.2| . The paths below the real axis will be trans- 
formed above the axis. The expressions split into a sin(i/j7r)-dependent contour integral 
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in the upper complex plane and into a sin(z/j7r)-independent straight-line integral from 
ioo(l + iSi) to — ioo(l + iSi). Depending on the choice of <jj, the sum ( |5.13|) becomes ex- 
actly one of the three expression (|F.15| ), ( |F.16| ) or ( |F.17| ), where the prefactor Wjtfi in 



App. |F.2| should be, of course, replaced by all the ^-independent terms of Eq. (|5.2|) and 



where are substituted by by ji-i, ji, ji+i- The angular momenta I and I" are here 

identified with s^i x and x respectively. After the Watson resummation of 
the other sums, e.g., of the U-i sum etc., / has to be replaced by and I" by Ui+i. 



If the penumbra scattering case |31| , |35| is excluded, the choice of o~i is, in fact, uniquely 
determined from the empirical constraint that the creeping amplitude has to decrease 
during the creeping process, as tangential rays are constantly emitted. In mathematical 
terms, it means that the creeping angle has to be positive. As discussed at the beginning 
of App. [F.2| , the positivity of the two creeping angles for the left and right turn uniquely 
specifies which of the three alternatives <7j is realized. In other words, the geometry is 
encoded via the positivity of the two creeping paths into a unique choice of the <7j. Hence, 



the existence of the saddle-point Q5.15Q is guaranteed 



The final step is the semiclassical approximation of the analog expressions to Eqs. (|F.15|) - 
( |K.17| ) as in App. [F.3HF.5| . Whereas the results for the creeping contributions can be di- 
rectly taken over from Eqs. ([F.34|) - (|F.43|) , there is a subtle change in the semiclassical 



evaluation of the straight-line sections. In the convolution problem of App. |F.3j and [b',5 



we have only picked up second-order fluctuating terms with respect to the saddle solution 
v' s from the v' integration. Here, we will pick up quadratic terms {}>i— vf) 2 from the ^ 
integration and mixed terms {v~ vt){vi±x— from the neighboring z/j_i and Ui+i inte- 
grations as well. Thus instead of having m one-dimensional decoupled Gauss integrations, 
we have one coupled m-dimensional one. Of course, also the saddle-point equations [the 
analog to Eq.(|F.35|) or (|F.44|)1 are now coupled: 



(1 - fy,i)2 arccos^/A;^] = arccos[(V* - z^_ x ) / 'kR^^] + arccos[(z/* - v? +1 )/kR jdi+1 ] 

~ i a n+ik - a n-in ~ ( 5 - 15 ) 

where the saddle v\ of the ith integration depends on the values of the saddles of the 
(z — l)th and (i + l)th integration and so on. Indeed, all m saddle-point equations are 
coupled. This corresponds to the fact that the starting- and end-point of a period orbit 
is not fixed from the outside, but has to be determined self-consistently, namely on the 
same footing as all the intermediate points. 

In order to keep the resulting expressions simple we will discuss in the following subsection 
just the geometrical contributions, and leave the discussion of the ghost and creeping 
contributions for later sections. 



5.4 Itineraries in the geometrical limit 



We will prove that the itinerary (— l) m A- 7U2 A- ?2J3 • • • A jm ~ ljm A jmjl leads, in the semiclas- 
sical reduction, to the following geometrical contribution: 
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oifc-Li-»m-i2m(7r/2) 

(-l) n A 3lj2 A» • • • A jm - ljm A jmjl | geom = — f ^ , (5.16) 

where the factor (— l) m results from the trace-log expansion Tdn(l — [—A]), as the peri- 
odic orbit expansion corresponds to this choice of sign. The quantity Li^ m is the length 
of the periodic orbit with this itinerary. Ai^ m > 1 is the expanding eigenvalue of the 
corresponding monodromy matrix and /ii_> m = 2m is the corresponding Maslov index 
indicating that the orbit is reflected from m disks (all with Dirichlet boundary condi- 
tions). Thus, for n-disk Dirichlet problems, the Maslov indices come out automatically. 
[Under Neumann boundary conditions, there arises an additional minus sign per disk 
label jj, since { {ka^)} / {±H$ {ka k )} ~ -H^(ka H )/ H^(ka H ) in the Debye ap- 
proximation. The minus sign on the right-hand side cancels the original minus sign from 
the trace-log expansion such that the total Maslov index becomes trivial. Otherwise, the 
Neumann case is exactly the same.] If the itinerary is the rth repeat of a primary itinerary 
of topological length p, the length, Maslov index and stability eigenvalue will be shown 
to satisfy the relations: L\^ m = rLi_ p , pi-, m = rp±^ p and Ai_ m = (Ai^ p ) r . 

Let us define the abbreviations 



di-i,* = jRl_ di - [(^ti " vl)/k} 2 = d^ x (5.17) 



Pi = ^al-(ut/ky (5.18) 

Li-\ t i = di-i 7 i — pi-i — pi = (5.19) 
5vi = Vi — v- 
~ 8ui 

QV. 



a i-l,i 

with i evaluated modulo m, especially i = is identified with % = m and i — m + 1 with 
i — 1. The quantity f^-i,, is the geometrical length of the straight line between the impact 
parameter v\_ x jk at disk and the impact parameter vf/k at disk ji in terms of the 
saddle points vf_ t and vf. The latter are determined by the saddle-point condition (|5.15|) 
which can be re- written for non-ghost scattering (o"j 7^ 1) as a condition on the reflection 
angle at disk jf 



0j. = arcsin [u? / ka 



arcsm 



'HI 



kR 



H-lH 



+ arcsin 



kR 



HH+l 



+ { a H+lH- a H-xk- <T i 2 ' K ) 



(5.20) 



Thus, pi is the radius of the disk ji times the cosine of the reflection angle and Lj-i^ 
is the geometrical length of the straight-line segment between the (i— l)th and ith point 
of reflection. Under the condition that the disks do not overlap, the inequalities Li_i ti < 
di-\ t i < Rj^ji hold and exclude the possibility that the reflection points are in the mutual 
shadow region of disks. For each itinerary there is at most one reflection per disk-label ji 
modulo repeats, of course. 



Then in analogy to App. F.5| the geometric limit of the itinerary ( |5.2|) becomes 
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•l) m A J1 -' 2 A-' 2 -' 3 • • • A Jm ~ i: ' m A Jmn l 
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3 ifc(di_l i— 2/^)— i7r 



if/4. 



i»/4 



X 



ddvi 



- e~ i7r / 4 oo 
j k J27=l L i-l,i~' lm7r 



-iTT/4 . 



\Di- 



11/2 



(5.21) 



where we have used that £™ x 2p^ = £™(Pi-i + Pi) since p = p m . Li^ m = YT=i is 
the total geometrical length of the geometrical path around the itinerary, see App. G of 
Ref . ]49 1 . Note that we used the saddle-point condition ( 5.15|) in order to remove not only 
the linear fluctuations, but all terms of linear order in the z/|'s from the exponents. Only 
the zeroth-order terms and the quadratic fluctuations remain. Di^ m is the determinant 
of the m x m matrix Fi_> m (=F) with 



Fi,i 
F 

1 m,m 

Fj.i+i 

F m ,i 
F- • 



2^m,l 

Pi 

2dj-i,i _ 
Pi 

2dm — 1 ,m 
pm 

d 1/2 . 

m— 1 ,m 
wV2 1 



1 



= ai 



i— l,i 



di,i+l 

dm — l,m 



— r 

- = Fi 
L l,m 

otherwise 



cij for 2 < z < m — 1 

"m 

for 1 < % < m — 1 

b m ,i 



(5.22) 



for m > 3. [For m = 2 the off-diagonal matrix elements read instead 

jl/2 



Fi, 



jl/2 

F _ "2,1 "1,2 _ 9 
* 2,1 — — — H — — — * ■ 



Jl/2 Ji/« 
"1,2 a 2,l 



1/2 



(5.23) 



The corresponding diagonal matrix elements are given as above, but simplify because of 
^2,1 = di t 2-] Thus in general, the determinant reads 



det 



( ai 6i >2 

&1.2 «2 &2,3 

6 2 ,3 "3 &3,4 

63 4 04 




\b m ,i 








&m,l \ 




' • m _2, m -l 

ci m —\ b m — i, m 

^m— l.m "m 



(5.24) 
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Note that determinants of this structure can also be found in Balian and Bloch and 
Berry . Our task, however, is to simplify this expression, such that the stability struc- 
ture of an isolated unstable periodic orbit emerges in the end. In order to derive a simpler 
expression for Di_». m , let us consider the determinant D^Xm °f the auxiliary mxm matrix 
F^ m (=F(°)) which has the same matrix elements as F with the exception that h m .\ = 0. 
The original determinant D%^ m can now be expressed as 

A-m = Dfl m - d -^Dfl m _ l + 2(-l)-+ 1 (5.25) 

where the last term follows from n^Li &i,i+i — 1- Here and in the following Dj® k is defined 
as the determinant of the auxiliary (k—l+1) X (k—l+1) matrix F^ fe with matrix elements 
F^ fc |ij = Fij for I <i,j < k(< m). Furthermore we define Z^+i-^ = F^^ = D^\ 
I . The D^ k determinants fulfill the following recursion relations 



(0) _ 


' 2dk-i,k 


j dk-l,k 


/—♦A; — 


\ Pk 


dk,k+l 


(0) _ 


f2dt- ltl 




Z-+fc — 


\ Pi 


^2,2+1 / 



\ Pi "2,2+1 / "2,2+1 

such that -D^°] fc can be constructed from all the lower determinants and Df2,j with 

I < i < j < k. For example, 

n (o) _ n (o) , ^fcA 2d m ,i n ( 0) , 2di, 2 (o) 

^l-fjfc — — 1 -^l-^fc-l ~' - V / 1 -^l-^O "I ^l-fl — • • • 

"fc,fc+l \ Pi P2 

Pi Pfc / 

(5.28) 

and 



(0) = + ( _ 1)fe -2 +1 ( j _ + ^2+1 (0) 



j% = --p^n.-i + (-1)™ i - -^A^-i + 

"fe,fc+l \ Pi Pz+i 



2-*/ 



+ --- + (-l) fe - Z+1 ^A ( 2 fc -il (5-29) 
Pk J 

as can be shown by complete induction. Note that the product dk y k+iDi^k is a multinomial 
in djj + i/pi where, for each index j, the djj+i/ Pi factors appear at most once. 

Replacing the u\\ m term in ( |5.25| ) by the r.h.s. of Eq.( |5.28| ) and using the relation ( |5.27| ) 
in order to simplify the expression 



-l,m / n (o) n (o) v 

m,l 
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recursively, we finally find after some algebra that 



n _ 2d w _i ;m n (p) 2(j m _2, m -i n (p) . yi _ l 2d Tn -i-i tm -j n (p) 

Pm Pm—1 Pm—i 

+ ... + (_i r -l^i jD (°) 
Pi 

2rf m — l,m r-)(0) , '^ J d m _i rn j-j(O) / -, \j— 1 2c? m — i )J7l t-j(O) 

-^2->m-l "i -^3->m-l ~ ' ' ' + l — -U U i^m-\ 

Pi P2 Pi 

+ --- + (-ir- 1 ^^^°l m -i. (5.30) 

Pm-1 

By complete induction it can be shown that Di^ m is a multinomial in Id^x^j pj of order 
m where the single factors appear at most once and the highest term has the structure 
n£Li 2di-i,i/ Pi- Thus, all the dj-i/s are in the numerators, whereas all the p^'s appear 
the denominators of this multinomial. We will show in Sec .6. 91 that 



-lf Al-»m + 



Ai. 



- 2 



(5.31) 



where Ai_» m is the expanding eigenvalue of the monodromy matrix which belongs to that 
period orbit which is given by the geometric path of the periodic itinerary. If the result 
of Eq. (|5.31|) is inserted into Eq. (|5.21|) the semiclassical reduction (|5.16|) is proven. 



5.5 Itineraries with repeats 



In the following we will discuss modifications, if the periodic itinerary is repeated r times, 
i.e., let m = p x r still be the total topological length of the itinerary, whereas p is the 
length of the prime periodic unit which is repeated r times: 



A' 



3U2 



P^H3i+\ . . . A 



Jmjl 



A 



nn 



••A 



(5.32) 



The length and Maslov index of the itinerary are of course r times the length and Maslov 
index of the primary itinerary A 7 ' 15 ' 2 • • • A- 7 ^ 1 , e.g., L\^ m = rL\^ p . The non-trivial point 
is the structure of the stability determinant Di_, m . Here we can use that the matrix 
Fi^ m has exactly the structure of the matrices considered by Balian and Bloch in Ref . J7|] , 
Sec. 6 D. Let Fi_> p be the corresponding matrix of the primary itinerary with matrix ele- 
ments as in Eqs.( |5.22"D [where m is replaced by p of course]. Following ref. J7| we furthermore 
define a new matrix F p (x) with matrix elements 



F p (x)|i,p = e" lx Fi_ +p |i iP 
F p (x)| P ,i = eix Fi^ P | P ,i 
F p (x)\i,j = F i^ P \i,j otherwise . 
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The determinant of the total itinerary is then |7j 



D^ m = J] £i^(e i2 ™ /r ) (5.33) 

n=l 



in terms of the rth roots of unity, since after r repeats, the prefactor in front of the (1, m) 
and (m, 1) matrix elements must be unity in order to agree with the original expression 
( |5.22| ). Let us furthermore define 



/3 P = \ (^V^i-p(O) - VIWtOJ (5-34) 
then, according to Balian and Bloch J7|, Sec.6D, 



D^ m =(a r p -(-(3 p y) 2 

= a 2 p r + ff-2(-a p P P Y ■ 

In our case we have the further simplification (in analogy to Eq.( |5.25[) ) 



-a p /3 p = ~(D 1 ^ p (0)-D 1 

=-i( 2 (-ir i -{-2(-ir 1 } 



as the corresponding two matrices differ only in the sign of the their (l,p) and (p, 1) 
elements. Especially we now have (3 P = (— l) p+1 /a p , such that 



= (ajr + -4^ - 2(-ir 



which corresponds to the usual form 



D^ m = (-ir (A^ p + J- - 2j (5.35) 

if Ai^ p is identified with (— l) p ctp. Note that from this the structure of Eq.( |5.31| ) follows 
for the special case r = 1. Thus we have achieved so far two things: we have proven that 
the determinant Di__ m organizes itself in the same way as a monodromy matrix does and, 
in fact, that it can be written in terms of a monodromy matrix M^p with eigenvalues 
Ai_ p , 1/Ai _+ p as follows 



D^ m = (-ir +1 det(MU - 1) . (5.36) 



38 



What is left to show is that Mi^ p is the very monodromy matrix belonging to the periodic 
orbit with the itinerary as in Eq. (|5.32|) . This will be done in Sec. [5.9| . But first, we will 
complete the study of the geometrical sector by deriving the ghost subtraction rules, and 
furthermore discuss periodic orbits with creeping contributions. 



5.6 Ghost rule 



Let us now imagine that the itinerary (|5.2|) has, at the disk position ji, an angular domain 
that corresponds to a ghost section, i.e. <7j = 1, i fixed: 



Because of the cyclic nature of the itinerary we can always choose the label ji away from 
the first and end position [remember that at least two disk positions of any periodic 
orbit must be of non-ghost nature]. In this case there are four changes relative to the 
calculation in Eq. (|5.21| ), see also App. [F.6| : first, the path of the Vi integration is changed, 



second, there is a minus sign, third, the saddle-point condition at disk ji is given by 
Eq.( |5.15|) with a' = 1 and not by (|5.20|) , fourth, the pi terms are absent. As in App. |F.6| , 



the saddle condition (|5.15|) at the jjth disk implies that dj-i^ + d iii+ i = c/j_ l i+1 . We can 



use this in order to express the length of the ghost segment Lj_ 1>i+1 between the reflection 
point at disk and the next reflection point at disk ji + \ in terms of the quantities 
defined in Eq. ( |5.19|) : 



— di-i.i+i — Pi-i — Pi+i — Li-i t i + 2pi + Li t i + i — pi-i — pi + \ . (5.38) 
Thus, by adding and subtracting the pi contributions we get 

A J1J2 . . . A ji—lji A 3iH+ 1 . . . \jmjl I 

1 * 1 * — 1 * — 1 * Igeom 



J] j d5vj \ J d5z/ l (-l)e- i -/ 2 e ifc2ft+i ^ ( ^ ) 



2_2_ 

Pi 



xTII e 1 \\-- 1/2— e ^ dl ~ h 

a — i7r/4 

( m p i7r/4 } m e ' 00 

UiV2^ ii=i 



. (5.39) 



\ D 9i 11/2 11/2 
I 1— >ml I I— mi I 

[Note that the exponent of the ghost itinerary is exactly the same as of the one of its parent, 
the same itinerary without the disk ji, whose geometrical path has the length Y4Z1 Li-i,i + 
£i-i,i+i+X}™i+2 ^1-1,1] I n writing down the last-but-one line we have cancelled the overall 
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minus sign by exchanging the upper and lower limit of the 5p{ integration. In addition, 
the following substitutions were applied: 

bvi =: d][\ j bvi for I ^ i, i + 1 , 

Sui =: e m ' 2 ( di - 1 > idi ' i +A 1/2 Sui , (5.40) 



In this way, the integration path and phase of the ith term agree with the ones of the 
other terms. Dfi^ m is the determinant of the mxm matrix Ffl^ m (=F 9i ) which is affected 
by the substitutions in the following way: 



F9i 

l,k 

F ft 

i,i 

F9i 

F 9 '- -, 



F; jfc for Z, A; ^ z,z + 1 
1 



\4-l,i4— l,i+l 
/ dj-iA 1 / 2 = 

V 4,2 + 1 / 



1/2 



*M-1 



r t+l,i 



4- 



iFj,i+i — iF 



(5.41) 



J?9i 



+l,i+2 



Pi+l di,i+l 4+l,i+2 

/ 1/2 
4-l,i+l \ ' 



/ 4-i,i+i \ 

V 4+l,i+2 / 



i+2,i+l 



where F/^ are the matrix elements as defined Eqs.( |5.22"|) ; i.e., 



D{ 



^ m = det (F* 



det 



i-l,t-l 



TT>9i 





1 F^ +1 

n f 34 F 9i F ft 

u r i+l,i r i+l,i+l r i+l,i+2 

Fj-Wi+l F i+ 2,j + 2 



V : 



(5.42) 



We now subtract the ith row times Ff i _ l from the (i— l)th and the ith row times Ff l i+1 
form the (i+l)th, as both operations leave the determinant Df i ^ m unaffected. Using that 
the ghost segments add, i.e., tZj-i^+i = dj-i,, + the numerators of the terms in the 

(i, i) and (z+l,i+l) matrix elements can be simplified. The determinant Dfi^ m , expressed 
via the transformed mxm matrix F 9i , reads 



D{ 



m m = det [F 9, 



det 



Fj-i^-i 


T?9i 
t i-l,i 


T?9i 

r i— 1,2+1 








1 










T?9i 


■pSi 


r i+l,i+2 








F 9i 


Fi+2,i+2 



V : 



7 



(5.43) 
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where 



TT>9i 

r i— 1 

1—1,4 

r t+l,i 

r 1 + 1,4 + 1 



2d, 



/ dj-2,i-i \ 
\ rfj— l,i+l / 



1 - 



1.2 



(5.44) 



2d, 



Pi+1 



- 1 - 



*j+l,i+2 



Note that we do not have to specify the elements on the ith row explicitly, as the ones 
on the ith line satisfy = 5{j. For the same reason we can remove the ith line and 
row altogether without affecting the result for the determinant. In doing so, we exactly 
recover the determinant -Di^i-i,j+ 2 -^m and matrix F 1 ^j_ li j +2 -^ m of the parent itinerary of 
the considered "ghost" . [The parent itinerary has the same sequel of disk indices except 
that the disk ji is missing.] 



D't 



det 



■ ' ' Fj+i^-i Fj +1>i+1 Fj +1> .j +2 

Fj + 2,i+l Fj + 2,i+2 



V 

detFx. 



7 



>i— l.i+2— »m 



£>i- 



►t— l,i+2— >m 



(5.45) 



The contribution of the ghost segment itself to the total "stability" of the itinerary in 
the geometric limit, i.e. to the stability factor of the corresponding periodic orbit, is just 
trivially one. As also the geometrical lengths and signs of both itineraries are the same, 
we have finally found that 



A 3132 . . . p^li-inpjiji+i . . . J^jmjl 



gcom 



A-? 1 "? 2 . . . pji-lji + l . . . pjmjl 



geom j 



(5.46) 



i.e., the ghost cancellation rule ( |5.5| ). Of course, the calculation of this section can trivially 
be extended to itineraries with more than one ghost (with and without repeats) as the 
operations in Eqs.( 5.39f) , ( |5.40 ) and ( |5.43 ) are local operations involving just the segments 
with disk labels ji, and ji+i- Thus they can be performed successively without any 
interference. Furthermore, as the transformations of the pairs (Si/^, 5vk+i) in (|5.40|) can be 
done iteratively (and in any order) for k — i, i + 1, • • -, the generalization to the extended 
ghost cancellation rule ( |5.10| ) is trivial as well: 



A 3132 . . . A ji-lji A 3iji+l A H+lji+2 . . . A jmjl I 

1 * 1 * — - * ^i- 1 * I geom 

A 3132 . . . A ji-lji \jiji+2 . . . A jmjl I 

- * -it- — ' * — : x I geom 

■ A jlj2 . . . A ji-lji + l A ji+lji + 2 A jmjl I 

Jr\. 1 \ * Jr\- | geom 

: A jlh ■ ■ ■ A j *- lj >+ 2 ■ ■ ■ A jmjl | geom (5.47) 
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etc. 



5. 7 Itineraries with creeping terms 



Let us now study an itinerary of topological length m which has, in the semiclassical limit, 
m — 1 specular reflections and a left-handed or right-handed creeping contact with one 
disk (which we can put without lost of generality at the end position), i.e. 



A 3132 A 3m-23m-l A 3m-l3m A 3mjl I 

" " |sc+creep • 



(5.48 



We mark those disk positions with creeping contributions by a tilde. Using the results 
and methods of App. |F.5| and Sec. [5.3| , we find the following result for the itinerary ( |5.48| ) 



sc+creep 



oo 



,itt 



I 12 Co 



a 



,]rn 



m-l,s m im 1-^1— »m— 1 



(0) 



i«m^ m Aaj m -7rcr m +arccos 



x e 



1/2 ~ 

1 — e i27r ^™ 



(5.49) 



Here, z/£ = ^e{kaj m ) is the £ m th zero of the Hankel function Hjp(kaj m ) in the upper 
complex v m plane [and Ci m = Ce m (kaj m ) is the creeping coefficient as given by Eq.( |F.33| ), 
see App. [F.4|| , A<x, m , a m and a' m are defined in Sec.|5.3|, pi and Li_i,j are defined in 

Eqs. Q5.171 ) - ( |5.19D and v\ is the solution of the saddle-point equation ( [5.15D [where in 
the cases i=l and i=m—l, the respective saddles v?_ x and uf +1 have to be replaced by 
s mV£m\- Furthermore, the following additional definitions have been introduced 



d m -i, Sm e m = s R 



Jm — ljm 



v" m - x - s m v tn 



2\ 1/2 



= L m _ lsm g m + p. 



m— 1 



d. 



R 



3™31 



1\ 1/2 



(5.50) 
(5.51) 



which correspond to the geometrical lengths to the surface of disk j m if ug is approximated 



by kdj m (see below). Finally, D^ m _ l is the determinant of the matrix F[' 



(o) 



*m— 1 
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with the matrix elements 



Fi i = - 1 - =ai for 2 < i < m - 2 

P _ 2d a m t m ,l _ 1 _ rf s m £ m ,l _ ~ 

1 ' 1 Pi di,a ~~ 

"p 2dm — 2 } m — 1 -I d m — 2,m — 1 ~ 

" m— l.m— 1 n i <7 „ — ™"m— 1 

rm—l u m-l,sm*m _ 

F v+ i = (f^)=F i+M =6 M+ i for2<z<m-2' ( ' j 

/ j' \l/2 

Fi, 2 = %7i^ = F 2,1 = 61,2 

\ m,l / 

Fjj = otherwise 

Note that D^ m _ l and Fj^. m-1 have exactly the form of the determinant D^\, m _ 1 and the 



matrix F^Xm-i defined in Sec.|T4] if the tilded lengths d Sm ^i and <i m _i jSm £ m are replaced 



by the "normal" geometrical lengths d SmTn i and d m -\ , 



d m ~i, Sm £ m ~ { R2 j m -dm ~ ( s ™a jm - v s m -xl k Y} = GU_i )Smm (5.53) 
<Lw ~ {^Li - ^ a im - "Ilk) 2 } = d„i (5.54) 
As discussed in App.|F.4j, this approximation is justified in the leading Airy approxima- 

1 /3 

tion, where terms of order h ' or higher are anyhow neglected. To this order we can 
approximate vi m everywhere by ka>j m , except in the "creeping" exponential, since there 
the error would be of order Ti~ 1 ^. Note that, in order to be consistent, we have to ap- 
proximate T>t m ~ ka,j m in the saddle-point conditions for v s m _ x and v\ as well. Thus, in 
this approximation, the saddles are manifestly real. Hence only in the overall factors in 
the exponents we keep the 0{h~ 1 ^) term of 



■ TV ( hjQ/ 



;j m 



6 J 1 m 

= ka jm + Su im . (5.55) 

For all the other terms the errors from neglecting 5vz m are, at least, of order 0({5Ve m } 2 / k) = 
0(H ' ) or even of 0(5V e Jk) = 0{h 2/3 ). The expansion of the products kd m -i >Sm t m and 
kd Sm e m ,i in the exponents leads to potentially dangerous linear terms of order 5ue m - How- 
ever, they cancel exactly against the terms in the expansion of the arccosines combined 
with those contributions which result if i>i m « kdj m is inserted into the saddle-point 
relations for v s m __ x and v\ . 

We will show below that in the leading Airy approximation, d m -\ Bm i m .jj corre- 

sponds exactly to the effective radius of the creeping periodic orbit -R^U m defined in 
Ref . p2| . The latter quantity is constructed, as in Eq. ([F.41|) , in terms of the length seg- 



ments U-x,i = Li-i^ between the [i— l)th and zth point of reflection (if 2 < i < m— 1), 
the length segments i TO) i = d mt \—pi and l m -i,m = d m -i, m ~ Pm-i between the (creeping) 
impact parameter at disk j m and the first or last point of reflection, respectively, and pi 
[= the radius of the disk ji times the cosine of the reflection angle]: 
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m— 1 

->eff 



R m ^ m ~ hnA IT + K ih,i+l) 
i=l 

= ^m^ m _i(l + «m-lVl,ro) > (5.56) 

where the curvature Kj is given by the recursion relation 

iU = —r^n + - (5-57) 

K i-1 + H— l,i Pi 

with 1/ko = 0. The proof of the equivalence of cf m _i )Srr j m |-Di°l m -il and i?f_^ m uses the 
following relations, which can be derived from Eqs. (|5.56|) and ( 5.57|) by complete induc- 
tion: 

(9 Rcff O Dcff 9 E>eff \ 

1 + + + • • • + ^Wl) /..^ (5.58) 

Pi P2 Pj-1 J 

K m^j-l \ Pi P2 Pj-1 / 

where R^^o = an d Rm^i = lm,i- For the right hand side of the proof, the recursion 
relations ( |5.28| ) are applied to the combinations \DfXj-i\dj-ij ~ D^Xj-idj-ij = R m <->j'- 



Rm^i — —Rm^i-i + (— l) 3 ' 1 ( 1 — — - H — — - 



Pi P2 



; 1!' ' 2/ '"' v : U; ., (5.60) 
Pi-i / 



with R m ^o = 0. In the induction assumption one can use that, for 1 < j < m. the 



quantities R m ^j and -R^Li are related as 



Rm^ = R^ + Pi ■ (5.61) 

This follows from the difference between and Ij-ij- Note that R m ^i = rf m ,i = -^m,i+Pi 
and R^^i = ^m,i satisfy trivially this induction ansatz. By complete induction it can 
now be shown that the recursion relation ( |5.60|) , applied to R m ^ m = -Di-l m -i^m-i,m, can 
be rewritten as 



f 9 f? e ^ 9 f? e ff 9 F? e ^ 

5 _ peflf i j I i | ZJX m^l , zrL TO^i , , ^^m^m-l 

-H-m^m -'i' m ^ m _i ~r 'm— l,m ) "T ' ' ' "T " " ' ~r 



4.9/ [ (-l)-- 1 -! (-l)- 2 ^!-^, 



Pi Pi Pm-1 

eff 



2 pi 
^ i-) J^m^i J^ m ^i _|_ _ _ _ _|_ ^m^m-l -fl- m ^ m _i 



Pi Pm-1 
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where the last bracket vanishes identically because of 



E> f? e ff f? f? e ff _|_ 9 Rcff 



Pi Pi-1 

for 1 < j < m. This equality can be derived from the induction ansatz Eq.( [5.61 ), if 
Eq.( |5.58 ) is inserted for the remaining R^^j on the left hand side and for the remaining 

i$*Li-i on the right hand side. Thus, the equivalence R m ^ m = Dfl m _ x d w _ 1>m = R^ m 
is established in the leading Airy approximation. 

We therefore get 



a 



co „i7r/12 /"( 

y^ e W m 

l m =ls m =±l (kdjj) 1 ' 6 V-^m' 



sc+creep 
\ 1/2 



,7 m 



JfcL 



m [ — >m 



(Sm) 



x e 

x 



1 



a,- — arccos 



3 i27r(fca 3m +5^ m ) 



_i/fc 



^| /fc 



JraJl 



(5.62) 



where L m ^ m = L m l + Y^=2 1 + ^m-\,m is the total length of the straight geometrical 
sections. The impact parameter vljk is given by the solution of the saddle-point equation 
( |5.15| ) where, in the cases i=l and i—m—1, the respective saddles vl_ x /k and z// +1 /fc, have 
to be replaced by s m aj m . 

In summary, all the quantities entering the semiclassical-creeping limit of the itinerary 
( |5.62j ) with one creeping section have geometrical interpretations: 

The integer index t > 1 enumerates the creeping modes around the boundary of disk 
j m . With increasing £, the impact parameter (or distance of the creeping path from 
the surface of the disks) and the "tunneling" suppression factor increases. 
The index s m = ±1 distinguishes between creeping paths of positive sense or negative 
sense around a surface section of disk j m . 

The coefficient e 17r / 12 Ci m /(kaj rn ) 1 ^ 6 is proportional to the product of the two creeping 
diffraction constants at the beginning and end of the creeping segment along the 
boundary of disk j m which parameterize the transition from a straight section to a 



ii 



m 



,iv 



(vi 



creeping section and vice versa, see [29,30,32 



The second prefactor is the inverse square root of the effective radius R^^ m , in 
units of disk radius a Jm . It is the geometrical amplitude, i.e., the geometrical stability 
factor. 

The i>i independent terms in the exponents are just \k times the sum of all lengths 
of the straight geometrical segments of the periodic itinerary, i.e., L m ^ m (s m ) = 

The geometrical length along the creeping section times ik is given by the sum of all 
exponential terms that are proportional to a Jm . 



45 



(vii) The creeping "tunneling" suppression factor is given by the imaginary part of U£ m or 

(viii) The denominator 1 — e l27TUl ™ results from the summation of all further complete 
creeping turns around the disk j m , in terms of a geometrical series |32fl . Note that 
the apparent poles at 1 — e l2nl/lm = cancel against the corresponding semiclassical 
poles of one-disk S-matrix, S 1 (A;aj m ). In fact, the zeros of 1 — e l27T1/im are given by 
vi m = /(integer) and are nothing but the zeros of the Hankel function H\ l \kdj m ) in 
the Airy approximation. 



5.8 More than one creeping section 



The rth repeat of the itinerary (|5.48|) follows simply as the sum, XifEs ra ! over the rth 
power of the summands on the r.h.s. of Eq. ( |5.49| ). As in the case of geometrical itineraries, 
this rule is trivial for the occurring prefactors, signs, phases and exponential terms. The 
non-trivial point is the behavior of the determinant D^Xm-i under the rth repeat. How- 
ever, as the corresponding matrix F^ m _ 1 has zero (l,m— 1) and (m—1,1) matrix ele- 
ments, such that repeats cannot couple here, the determinant of the rth repeat corresponds 
exactly to the rth power of the determinant of the primary itinerary. For the same reason, 
also the determinants and corresponding effective radii of itineraries, with more than one 
creeping contact (i.e., with at least two disks j\ and j m with creeping contacts), decou- 
ple from each other. The corresponding semiclassical result for such an itinerary is thus 
the multiple sum, J2e t Y^ Si J2e m J2 Sm over the products of the corresponding itinerary from 
disk ji to disk j m and the itinerary from disk j m to disk ji, each individually given by the 
suitably adjusted summand on the r.h.s. of Eq.( |5.49| ), e.g.: 



£3132 . . . J^3i-l3i J^3i3i+l . . . p^3m-l]m J^JmJl 

oo oo 

~E E E E 

4=1 s i =±ie m =i s m =±i 

x (-1 



1/2 



sc+creep 



1 



(fci*) 1 /" I J2* 



I _ g2n(kaj.+6i>£.) 



x e 



i(kaj i +5f£. ) I (2—rj / i — Si(Ti)iT+SiAaj i —arccos 



a A . — s 



1—1' 



a, a. —s 



V s ,, Ik 



R 



x(-l) 



1/2 



x e 



/ fa,,— sm^^i/fc 

i{ka jm +5vt m )l (2-a' m -s m cr m )n+s m Aaj - •■ ^ ■ ^- ' 



"f/fc 



etc. [If there are two (or more) creeping contacts next to each other, e.g., ji = j m -i, then, 
in the above formula, the corresponding impact parameters v? +1 /k and v s m _ x jk have to 
be replaced by a Jm and aj m _ 15 respectively] 

The physical reason for the simple rule of piecing together creeping paths, is the point-like 
contact at e.g. disk j% between the creeping sections on the one hand and the geometrical 
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sections on the other hand which is mediated by the diffraction constants Ci . . [Mathemat- 
ically, this corresponds to the fact that v^. is uniquely determined as the ijth zero of the 
Airy integral and not by a semiclassical saddle-point equation that would couple with the 
saddle-point equations at the disks and ji+x] Because of this point-like contact [the 
independent determination of V(\ the semiclassical itineraries multiply for fixed value of 
the mode numbers £j and creeping orientation Sj. Especially, if we limit the mode number 
to i = 1, periodic orbits with common creeping sections can exactly be split up into their 



primary periodic orbits, see Ref. [[32 



Finally, the ghost cancellation works for itineraries with creeping sections in the same way 
as for itineraries which, semiclassically, are purely geometrical. The reason is two-fold: 
First, by construction (see App. [F.5| ), ghost segments can only occur in the geometrical 



part of the creeping itinerary. Second, the ghost cancellation rules of Sec.[5]6] are based 
on the local properties of the segments i— 1 — > i and i — > Let us now assume, 

for simplicity, that the disk ji is cut by the ghost section. If there is no creeping at the 
neighboring disks and ji+i, the reduction of the stability matrix F^ ^ and of the phases 
and lengths of the segments is precisely the same as in the purely geometrical case (see 



the substitutions ( [5.40D and the analogous steps of Sec.|5.6|). If there is a creeping contact 



at disk or/and disk the substitutions ( |5.40| ) simplify, as bv^\ or/and <5z4+i do 
not exist. Thus, the elements Ffl l i or/and Ff+ lji+1 of Eq.( |5.41| ) are zero and the ith row 



of the determinant ( |5.42| ) has only to be subtracted from the (i + l)th or the (i — l)th row 



or from none, in order that Eq. (|5.42|) becomes Eq. fl5.43p . The reduction in the lengths 
and phases hold in these cases as before. 

In summary, the ghost cancellation works for geometrical orbits with creeping sections 
as well as for purely geometrical orbits, studied in Sec.[5.6[ Semiclassically, neither the 
ghost itineraries nor their parent itineraries [which have the same symbol sequence except 
that the ghost labels are removed subsequently] contribute to the semiclassical trace-log 
expansion and to the cumulant/curvature expansion. Thus, one can omit these "ghost- 
affected" periodic orbits altogether from the curvature expansion. 

Deep inside the negative complex fc-plane the limitations of the first Airy correction 
introduce rather big errors, see Ref . |67| . In this case it is advisable to use to the original 
expression (|5.49|) for semiclassical "creeping" itineraries with Vi m and Cg m , as given in 
Eqs.flra|) and ( jFgg) , instead of Eq.(pbp. 



To summarize, for the special case of n-disk repellers, the creeping periodic orbits of 



Ref. [32 1 have been recovered directly from quantum mechanics, whereas the construction 



of Ref.p2| has relied on Keller's semiclassical theory of diffraction |3"0fl. Furthermore, the 



symbol dynamics has to be generalized from the single- letter labelling {ji} to the two- 
letter labelling {ji, Si x £i} with Sj = 0, ±1 and ^ = 1, 2, 3, • • -. 



5.9 Geometrical stabilities 



In this subsection we will return to purely geometrical periodic orbits and show that 
Eqs. (|5.35p and (|5.36|) are correct, i.e. that the determinant Di^ m satisfies in fact 
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-l)^[A 1 _ m + -^—-2) (5.63) 



^ m = (-l) m+1 det(M 1 _ m -l 

Ai 

irrespective, whether there are repeats or not. Here M^ m is the 2x2 dimensional real 
monodromy matrix of the purely geometrical periodic orbit of total topological length m 
(m = pr if there are r repeats of a primary orbit of topological length p), that is, the 
semiclassical limit of the itinerary A^ 1 ■ ■ ■ A 7 ™- 71 . Because of phase-space conservation, 
the determinant of Mi^ m is unity. For this reason and as the matrix elements of 
are real (see below), the two eigenvalues of the matrix are related as Ai^ m and 1/Ai_> m . 
We do not have to treat repeated orbits explicitly here, as this case was already studied 



in Sec. 5.5 



In Ref.||15|] it was shown that, for any two-dimensional scalar billiard problem (whether a 
bound state problem or a scattering problem), the monodromy matrix M^ m of a periodic 
orbit with m collisions with the billiard walls is given by the 2x2 dimensional Jacobian 
belonging to the infinitesimal evolution of the vector (Sp±,5x±) T perpendicular to this 
classical trajectory in phase space, i.e., by the product 



Mi_ wl = nT<_ M R i . (5.64) 

8=1 

Here the matrix 

! 1 \ (5.65) 




parametrizes the translational (straight ray) evolution of the vector (Sp±, Sx±) T [or rather 
(S6 p , Sx±) T with 9p being the angle of the momentum p, since the modulus of p is conserved 
anyhow] between the (i— l)th and ith collision where Lj_i 5 j is the corresponding length 
segment. As usual i=0 should be identified with i—m. The matrix 




Ri = n ' r ' (5.66) 



parametrizes the evolution of the vector (56 p ,5x±) T from immediately before to imme- 
diately after the ith collision. The quantity pi=aiCos9i is, in general, the product of the 
local radius of curvature a« and the cosine of the reflection angle 9i at the ith collision with 
the billiard walls. Especially for our n-disk scattering problems, a, is of course nothing but 
a,j v the radius of the disk jj, whereas $i should be identified with the scattering angle 9j i 
of Eq.( |5.20| ), the solution of the saddle-point equation. Since the determinants detTj.^ 
and detRj are trivially unity, the determinant of the product Mi^ m is unity as well, as it 
should because of Liouville's theorem. Furthermore, the matrix elements of Mi^ m have 
to be real, since the matrices detT^i^ and detRj are real, by definition. Thus the two 
eigenvalues of have the structure A!^ m and 1/Ai_,. m . 
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5.9.1 Monodromy matrix in closed form 

In the following we will construct a closed-form expression for the matrix elements of the 
matrix Mi_, n , 1 < n, by complete induction. Let us denote these matrix elements as 



Mi_>„ 



A n B n 
C n D n 



(5.67) 



By inserting Eqs. (|5.65 ) and ( |5.66 ) into Eq. (|5.64 ), one can show that 



1 



2_ 

Pi 

D x = 1 



Pi 



and that the matrix elements of M^ n+1 and are related as follows: 



A n +l 


= A n 


+ B n L Utn+ i 


B n +\ 


= B n 


I 1 I ^^ n ' n+1 

V Pn+l 


Cn+l 


= C n 


+ D n L nj1l +i 


Dn+1 


= D n 


1 2L„ ))1+1 

\ Pn+l 



+ A r 



+ C n - 



Pn+l 

2 

Pn+l 



(5.68) 



(5.69) 
(5.70) 
(5.71) 
(5.72) 



In order to be able to perform the induction step, we do not make use of the cyclic 
permutation, i.e., in the following we do not replace L 0) i with L n \ or L n +i,i, respectively, 
and L n n+ i or L n+ x,n+2 with L n l , but keep the original labelling. From Eqs. ( |5.68|) and 
( p. 72] ) it follows, by complete induction, that 



71—1 9 peff 

i=t pi 

n I no peff 

B.. = E 1+ £ ^f' 

i=l \ j=i+l Pi 
ri peff 



Dr 



n 9 peff 

h P* 



9 n peff Reff 9 



j=l 



Pi 



(5.73) 

(5.74) 

(5.75) 
(5.76) 



with -R^„ given as in Eq.( |5.56j ) where we identify k-i^ with L^+j-ij+j. In analogy to 
Eq.( |5.58 ), we can derive the recursion relation 



R 



.eff 



P eff 

' j— *n— 1 



n-l-j 2 

1+ J2 R f-i+i— 

i=l Pj+i 



J n—l,n 1 



(5.77) 
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where Rf^ n \j> n = 0. Thus R^ n should not be mixed up with the quantity R c ^-, n °f 
Sec.|5]7] that rather corresponds here to R^ n with Z ,i = £0,1, of course. Note that the 
first iteration of Eq.( |5.77| ) leads to R^j +X = Ljj +X . For later purposes we also define 
here the effective radius Rf!_ n which is, of course, equal to Rfi> n and which satisfies the 
recursion relation 



n-X-j 2 

RC j^n = R< j+l^n + \ 1 + E ~ R n-i+-n | > (5.78) 

i=l Pn-i 

where again R c ^_ n \j> n = 0, such that R^_ x< _ n = L n _ Xn . 



The second equation of Eq.( 5.74 ) follows trivially from ( 5.77|) . By inserting the Ansatze 



( |5.75| ) and ( |5.76| ) into the induction step ( |5.71|) , one can easily show, with the help of the 
recursion relation ( |5.77|) (for the case j =0), that the result for C n+X is given by Eq. ( |5.75|) , 
with n replaced by n + 1. Similarly, by inserting the Ansatze (|5.73|) and the last identity 
of ( |5 . 74| ) into the induction step ( |5.69| ), one finds that A n+X is compatible with ( |5.73| ). 
Here we used the identity R^ n = 0. Applying the recursion relation Q5.771 ) to R^ n+1 , 
it easy to show that D n+X is compatible with Eq.( |5.75| ) as well. Finally, for proving that 
B n+ i is compatible with Eq. ( |5.74| ), one inserts the first equation of ( |5.74j ) and ( |5.73|) into 



( |5.7(J| ), uses Eq. ( [5. 77] ) for re-expressing R^. n+X and the fact that L n n+1 = R^ nJrX . 

Having a closed form expression for the matrix elements of M^,,, we could now construct 
the corresponding eigenvalues A.f^ m . But, in fact, we only need the linear combination 
(— l) m (Ai^ m + 1/Ai^ m ) which is equal to the sum A m + D m . 

In summary, we have now a closed form expression for the right hand sides of Eq. ( |5.63| ) 



-l) m+1 det(M^ m - 1) = (-l) m ( A^ m + ' 



2 ^2 



= 2(1 + (-ir +1 ) + Y, R&i- + E -Rtm , (5-79) 

i=X Pi i=l Pi 

where we used the identity R^ m = RfL m in writing down the last relation. 
5.9.2 Stability determinant in closed form 

In analogy to the definitions of Sec.^T^ [see Eq. ( |5.60| )] we define here Dy^ x _^ k _ x dk-x,k = 
Ri^k = Ri^-k- These quantities satisfy, according to Eq. (|5.26|) , the recursion relations 



- - (- 2 r- - 1 1 ^ 

Ri^k + Ri-^k-x = Ri->k-x Ri-*k-x + Ri^k-2 3 dk-i± (5.80) 



Pk-l 

and, according to Eq.( |5.27|) , the recursion relations 



'ft-2,fc-l . 
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Ri^k + Ri+i^k — di,i+\ — Ri+i<-k — j 

\Pl+l «/+l,i+2 



Ri+i<-k + R, 



By complete induction, these recursion relations can be summed up to give 



Ri- 



Ru 



V i=l Pl+i J 



-R 



M+l 



/ fc-J-1 

1+ E(-i 

\ 1=1 



Pk- 



Witll -R;^fc|/>A: — -R;^fcL>fc — 0. 



(5.81) 



(5.82) 
(5.83) 



According to Eq. ( |5.30|) of Sec.|5.4|, the stability determinant -Di_> m can be rewritten in 
terms of the i^_>fc's as follows 



m r\ m— 1 o 

Di^m = EMr - *^- + E • 

8=1 i=l Pi 



Adding and subtracting Eq. ([5.79|) to ( |5.84j) we get 



(5.84) 



m— 1 



^ m = 2(i+(-ir +1 ) + E^-+ E 



8=1 
m 



Pi 



8=1 



Pi 



m—1 



+ 2((-i) m -i) e ((-i) m ~ 4 R ^-RoQ - + E - ((-i)' 

8=1 A 8=1 A 



eff 

i<— m 



(5.85) 



The equality fl5.63| ) is established, if we can show that the second line of ( |5.85| ) is identically 



zero. 



Note that the effective radius R m ^j, in the creeping case, fulfills the recursion relations 
( |5.82| ) and ( |5.83|) as well, see e.g., Eq. (|5.60|) . However, as here do,i(= d>m,i) — Po + ^o.i + Pi 
and d m _i, m = p m _i + L m -i,m + d m , whereas in Sec|T7| d ,i = A),i + p x and cf m _i, m = 
p m _i + L m _i jm , the relation between the Ri-^/s and the -Rff^-'s have to be modified in 
comparison to the relation ( |5.61| ) between the R m ^/s and the .R^L^'s. In fact, instead of 
Eq.( |5.61|) , we get 



i?o^i = 1 1 + Po 

-Rj^m = S 1 + Prr 



5 



<9L ,i J 
9 



ncff , "^O^i -"fl-H-l 

i<o-H H 7 

> -kj-l.i 



Pi 



pefF I - n '8<— m -"'i+l-i— m 



Pi 



(5.86) 
(5.87) 



where the differentiations with respect to Lo,i and L m -i j7n produce the additional po 
and p m pieces in <i ,i and d m _i, m , respectively, relative to Eq. ( |5.61| ). As in Sec.|5.7|, these 
relations can be proven by complete induction. 



51 



Now, by solely inserting Eqs.( |5.86| ) and ( |5.87| ) into the second line of Eq.( |5.85| ) and col- 
lecting terms, we get for this expression 



2 ^2 



2((-i) m -i) + E ((-i) m_i Ro^i-Rf^) - + E- Rn-m-R. 

Pi i=1 Pi 



■eft 
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m— 1 r> peff v 



1=1 



Pi 



(5.88) 



With the help of the recursion relations ( |5.77|) and ( |5.78| ) this expression can be rewritten 
as follows 
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(5.89) 



Under p Q =p m , Ro^ m =R ^ m and Rfi^ m —R^_ m , the expression simplifies [furthermore, 
note that R^_ i4 _ m and -R^i are independent of L ,i and L m _i )m , respectively, if i < m — 1]: 
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(5.90) 



The identity ( |5.63|) is therefore established. In summary, we have proven that the geomet- 
rical semiclassical limit of a quantum itinerary for any non-overlapping n-disk system [see 
Eq.( |5.21[) 1 is exactly the corresponding periodic orbit with the Gutzwiller weight. Hence, 
the validity of Eq.( |5.16|) for any non-overlapping finite n-disk system (with the exclusion 
of the grazing geometries) is shown in the semiclassical limit. Note, however, that this is 
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no general proof of the convergence of that the curvature series, since two limits are in- 
volved: the semiclassical limit p/Ti = k — > oo (or h — > 0) and the cumulant limit m — > oo. 
In general, these two limits do not commute. For purely chaotic classical n-disk systems 
with a positive value for the topological entropy, the exponential proliferating number of 
orbits and, therefore of classical input, is not compatible to the just algebraically rasing 
number of operations, needed to solve for the zeros of the quantum determinant of the 
multi-scattering kernel. In these cases, the curvature sum of the periodic orbits has to 
deviate from the cumulant sum involving the quantum itineraries. The semiclassical limit 
and the cumulant limit should better not commute. We will study this numerically in the 
next section. 
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6 Numerical tests of semiclassical curvature expansions against exact data 



In this section which overlaps partly with Ref. |7[ we test the predictions of the curva- 
ture expanded Gutzwiller-Voros zeta function, the dynamical zeta function |^6| and the 



quasiclassical zeta function of Refs. [12.53] against the exact quantum-mechanical data for 
the 3-disk-system in the A i -represent at ion. 



As mentioned in the introduction, the 3-disk repeller [ p!0| , p!l| ,p^- p!5|j42]] is one of the sim- 
plest, classically completely chaotic, scattering systems and provides a convenient numer- 
ical laboratory for computing exact quantum-mechanical spectra as well as for testing the 
semiclassical ideas. It consists of a free point particle which moves in the two-dimensional 
plane and which scatters off three identical hard disks of radius a centered at the cor- 



ners of an equilateral triangle of side length R, see Fig. |6.1| . The discrete C3 V symmetry 
reduces the dynamics to motion in the fundamental domain (which is a l/6th slice of 
the full domain and which exactly contains one half of one disk), and the spectroscopy 
to irreducible subspaces A 1; A 2 and E. All our calculations are performed for the fully 
symmetric subspace Ai [13|J34 . 




R 



R 





R 



Fig. 6.1. The three-disk system with center-to-center separation R = 6a. 



The genuine multiscattering data in the Ai subspace are computed from the determi- 
nant det[l + Aa^)] where the multiscattering kernel Aa^), expressed in the angular 



momentum basis relative to the half-disk in the fundamental regime, reads |13 



A Al (fc)m,m' =d(m)d{m) {cos ^(5m - m')j H { ^_ ml {kR) 



H^(ka) 



(6.1) 



with < m, m' < oo and 



d{m) :- 



a/2 for m > 
1 for m = . 



As A is trace-class for any n-disk geometry, the determinant exists and can numerically 
be calculated in a truncated Hilbert space. The Hilbert space is here the space of angular 
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momentum eigenfunctions on the surface of the half-disk in the fundamental domain 

which can be truncated by an upper angular momentum m max . From the study, in App.0, 
of the asymptotic behaviour of A m m / with respect to the angular momentum one can 
derive the following inequality for the truncation point m max : 

g 

»max ~ -jka « l.hka . (6.2) 

This agrees, of course, with the numerical findings. The truncated matrix Ma 1 (A;) is then 
numerically transformed to an upper triangular matrix and the determinant is calculated 
from the product of the diagonal elements. This procedure is faster than the computation 
of the determinant from the product of the eigenvalues of Ma^) (see ( |A.10| )). The 



numerical results for both ways agree, of course, up to computer accuracy. The zeros 
of the determinant, detMA 1 (fc), in the lower complex wave-number plane determine the 
scattering resonances, whereas the phase of the determinant evaluated on the real fc-axis 
gives the cluster phase shift. The cumulants can be constructed either from the Plemelj- 
Smithies recursion formula ( |A.14j ) or from the multinomials of the eigenvalues ( |A.11|) . The 
latter procedure is numerically more stable, especially deep inside the negative complex 
wave-number plane. This concludes the numerical setup for the exact calculation. 

As shown in Secs.[| and [5] the classical analog of the characteristic determinant (actu- 
ally of det{l + zA(k)} to be precise) is the semiclassical zeta function of Gutzwiller 0] 
and Voros |37j which, prior to a regularization, is given by Zcv( z ',k) (see ( |2.1|) ). How- 
ever, in the literature there exist other competitors for a semiclassical zeta function, e.g., 
the dynamical zeta function (q 1 ^; k) of Ruelle |]36]| (see ( |2.2|) ) which is the j=0 part of 



the Gutzwiller- Voros zeta function as well as the quasiclassical zeta-function Z qc \(z;k) 
of ref.||| (see (|01)). As usual, for all three choices, t p (k) = e lkL p-' lv v^l^ /\^ p \\ i s the p th 



primary cycle, n p its topological length, L p is its geometrical length, v v its Maslov in- 
dex together with the group theoretical weight of the studied C 3v -representation (in the 
present case the Ai-representation), and A p its stability (the expanding eigenvalues of the 
stability matrix) — see refs. [p^| , |43| for further details. The variable z is a book-keeping de- 
vice for keeping track of the topological order in the cycle- or curvature expansion [ |18| , |19H 



(see Eqs.(p.4|) and (fO|)). In the following, the various curvature-expanded zeta functions 
are truncated at a given curvature (i.e., total topological) order n c . The semiclassical pre- 
dictions for the scattering resonances are determined from the zeros of these truncated 



zeta functions, the predictions for the cluster phase shifts discussed in Sec. |6.2| from the 
phases on the real fc-axis and the curvatures from the terms of order z m in the curva- 
ture expansion. Input data for the lengths L p , stabilities A p and Maslov indices v v of 
the periodic orbits of the 3-disk system in the Ai-representation have been taken from 
Rosenqvist p8lp9j, Scherer JTIJ and Eckhardt |70 . 



6.1 Exact versus semiclassical resonances 



In this chapter we compare the numerically computed exact quantum-mechanical reso- 
nances of the 3-disk repeller with the corresponding semiclassical predictions of the three 
semiclassical zeta functions: the Gutzwiller- Voros zeta function ( |2.1|) , the dynamical zeta 
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function ( 2.2 



) and the quasiclassical zeta function (|2.3| ). 



For the 3-disk-repeller with center-to-center separation R=6a, we have computed all exact 
quantum-mechanical Ai resonances (numerically determined from the zeros of defM^ {k)) 
as well as all the corresponding approximate ones (from the zeros of the at finite curvature 
order n c truncated zeta functions) in the wave- number window: < Refc < 250/a and 
> Irak > — 1.6 /a. This window contains several hundreds of leading and sub leading 
resonances, from the lowest ones onwards. In Figs. |G.l|-|G.9l, for increasing curvature order, 



the resonances are plotted as the real part of the wave number (resonance "energy" ) versus 
the imaginary part of the wave number (resonance "width" ) . 



Some features of the resonance spectra allow for an immediate interpretation p~3|JTT|JT5|1 : 



The mean spacing of the resonances is approximately 2n/L, where L is the average of 
the geometrical lengths of the shortest periodic orbits, namely the lengths L and L\ of 
the two periodic orbits of topological length one. The data also exhibit various beating 
patterns resulting from the interference of the periodic orbits of nearly equal length; e.g., 
the leading beating pattern is of order 2ir/AL, where AL is the difference of the lengths 
Li and L . 



In Figs.|G.l] - |G.4| a comparison is made from the first- up to the fourth order in the 



curvature expansion. Already at fourth order the four leading resonance bands are well 
approximated by the Gutzwiller-Voros zeta-function (in fact, for Refc ^75/ a already the 
second curvature order is enough to describe the first two leading resonance bands). This 
is in agreement with the rule of thumb that any new resonance band is linked with a new 
curvature or cumulant order. Neither the dynamical zeta-function nor the quasiclassical 
one perform as well to fourth order. The reason is that the quasiclassical as well as the 
dynamical zeta-function predict extra resonances which are absent in the exact quantum- 
mechanical calculation. Thus the third and fourth curvature order of these zeta-functions 
are distributed over the average of the third and fourth resonance bands and the spurious 
extra resonances. In the window plotted one can classify the exact data into four leading 
resonance bands closest to the real wave-number axis and 2 subleading ones shielded by 
the leading resonances. Thus, just periodic orbits of topological length up to four are 
needed in order to reproduce the qualitative trend of the exact data closest to the real 
axis. The 3-disk-system has 8 periodic orbits up to this topological length. Actually, the 
3-disk-system with center-to-center separation R = 6a is not very chaotic at these k 
values. All experimentally accessible spectral data in this regime (which can be extended 
up to Re k ps 950/a as only about there the subleading resonance bands mix with the four 
leading ones) can be parameterized by 16 real numbers, i.e., 8 periodic orbit lengths, 8 
stabilities, and 8 Maslov indices. Experimentalists can stop here. The subleading bands 
are completely shielded (up to Refc ~ 950/a) by the above mentioned four bands. The 
subleading bands (below Refc ~ 950/a) are only of theoretical interest, as they can be 
used to test the semiclassical zeta functions. 



In Figs. G.5 a comparison is made up to fifth curvature order. The Gutzwiller-Voros 



zeta-function does at least as well as in Fig.|G.4|a for the leading four resonance bands 



but now it also describes the peak position of the fifth resonance band for large enough 
values of Kek. Note the diffractive band of exact resonances from k ps (0. — i0.5)/a to 
k ps (100. — il.6)/a which our semiclassical zeta functions fail to describe. As shown 
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in Refs. |p2|-434f the diffractive band of resonances can be accounted for by inclusion of 
creeping periodic orbits which have been omitted from our semiclassical calculations. The 
dynamical and quasiclassical zeta functions show a slight improvement with respect to 
the four leading resonance bands; however, no agreement with the fifth one. 



In Figs. |G.6| a comparison is made up to sixth curvature order. The Gutzwiller-Voros 



zeta-function fails for the third resonance band below Re A; ~ 20/ a, for the fourth below 
Re/c m 45/ a, for the fifth and sixth below Kek m 70 /a and 80/ a, respectively. Below 
these values, the last two curvature orders try to build up an accumulation line. Above 
these values, the qualitative agreement with the data is rather good. The dynamical zeta- 
functions at this order just improves the description of the four leading resonance bands. 
Furthermore, it builds up a sharp line of accumulation for the subleading resonances, the 
border of convergence controlled by the location of the nearest poles of the dynamical zeta 
function, see Pl"P2"| . The quasiclassical zeta function also improves the description of the 



four resonance bands, although it is still not of the same quality as the Gutzwiller-Voros 
one even at two curvature orders lower. Note that the quasiclassical zeta function is trying 
to build up two bands of spurious resonances in agreement with our rule of thumb. 



In Figs. |G.7| a comparison is made up to seventh curvature order. The first four resonance 



bands of the Gutzwiller-Voros zeta function have converged and the accumulation line 
has moved up. Only above Re/c ~ 140/a the fifth and sixth resonance band emerge, now 
with improved accuracy, however. Also the seventh resonance band is approximated. The 
dynamical zeta-function now clearly produces its line of convergence (the accumulation 
line of resonances). Above this line, the resonances (except the ones very close to the accu- 
mulation) are approximated as well as in the Gutzwiller-Voros case; below, no agreement 
is found. At this order the quasiclassical zeta-function is doing as well as the Gutzwiller- 
Voros zeta-function did already at curvature order four. None of the subleading bands are 
described by the quasiclassical zeta-function. Instead another band of spurious resonances 
emerges. 

In Figs. |G.8] and |G.9| the comparison is made up to the eighth and twelfth curvature or- 
der, respectively. The border of convergence of the Gutzwiller-Voros zeta-function has 
now moved (in the plotted region) above the fifth and sixth band of the exact resonances. 
It has moved also closer to the very sharp accumulation line of resonances of the dynam- 
ical zeta function. However, these lines are still not identical even at twelfth curvature 
order. The subleading quasiclassical resonances have stabilized onto the spurious bands. 
Furthermore, some subleading resonances move further down into the lower complex k- 
plane E^. Eventually (see also [0), starting with curvature order 10 and 12 the fifth and 
sixth resonance bands are approximated — in addition, to four or six spurious resonance 
bands, respectively. Thus the quasiclassical zeta-function seems to find the subleading 
resonance bands, but at the cost of many extra spurious resonances. Note that at these 
high curvature orders the quasiclassical zeta-function has numerical convergence prob- 
lems for large negative imaginary k values (especially for low values of Re A;). This is 



# 2 Note that the quasiclassical results of this figure are directly comparable with the results 
of the so-called Quantum- Fredholm determinant of Ref.[^] (see Fig. 4b in Ref . ) as both 
calculations involve periodic orbits of topological length up to eight. As we now know, all the 
subleading resonances of that figure have nothing to do with quantum mechanics. 
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in agreement with the expected large cancellations in the curvature expansion at these 
high curvature orders. Furthermore, periodic orbits of larger topological order than twelve 
would be needed to falsify the success of the quasiclassical zeta function, since it barely 
manages to approximate the two bands of subleading resonances at this curvature order. 

Qualitatively, the results can be summarized as follows. The Gutzwiller-Voros zeta-function 
does well above its line of convergence, defined by the dynamical zeta-function, already at 
very low curvature orders where the dynamical zeta-functions still has problems. Below 
this line we observe that the Gutzwiller-Voros zeta-function works only as an asymp- 
totic expansion. However, when it works, it works very well and very efficiently This 
implies that the additional (A p -dependent) terms of the Gutzwiller-Voros zeta function, 
relative to the simpler dynamical zeta function, are the correct ones. This is of course 
in agreement with the findings of our semiclassical reduction in Sec.[5| Eventually, the 
dynamical zeta-function does as well for the leading resonances as the Gutzwiller-Voros 
one. As experimentally these are the only resonances accessible, one might - for prac- 
tical purposes - limit the calculation just to this zeta function, see, however, Sec.|6.2|. 
The quasiclassical zeta-function seems to find all subleading geometrical resonances. Un- 
fortunately, the highest periodic orbits at our disposal are of topological length 12; the 
very length where the sixth resonance band seems to emerge. Thus higher orbits would 
be needed to confirm this behavior. But all this comes at a very high price: The rate of 
convergence is slowed down tremendously (in comparison with the asymptotically work- 
ing Gutzwiller-Voros zeta-function), as this zeta-function is producing additional spurious 
resonance bands which do not have quantum-mechanical counter parts, but only classical 
ones f47j. Without a quantum calculation, one could therefore not tell the spurious from 



the real resonances. 

As a by-product we have a confirmation of our empirical rule of thumb that 'each new 
cumulant or curvature order is connected with a new line of subleading resonances'. This 
rule therefore relates the curvature truncation limit, m —* oo, either to the limit Irak — » 
— oo, if there is no accumulation of subleading resonances, i.e., if the zeta function is 
entire [f|^,(i3|], or to the formation of an accumulation band of resonances. Both facts 



support our claim that, in general, the curvature limit m — > oo and the semiclassical limit 
Re k — > oo cannot and should not commute deep inside the lower complex fc-plane, as the 
subleading resonances of increasing cumulant order are approximated worse and worse. 
Only an asymptotic expansion should be possible, in agreement with our findings for the 
Gutzwiller-Voros zeta function. 



6.2 Exact versus semiclassical cluster phase shifts 



In the last chapter the semiclassical zeta functions were judged by the comparison of 
their resonances predictions with the exact resonances poles (especially the subleading 
ones), as was done in the past, see e.g. 



Refs. 13,15,42lB4lHg 



4fl. Since the deviations 



between the zeta functions themselves and from the exact data are most pronounced for 
the subleading resonances (which are shielded by the leading ones), one could argue that 
empirically it does not matter which of the three zeta functions are used to describe the 
measured data, since all three give the same predictions for the leading resonances |42]J43[ 
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Below, however, we will show that even experimentally one can tell the three semiclassical 
zeta functions apart and that, in fact, the Gutzwiller-Voros one is by far the best. 



6.2.1 Cluster phase shifts 

In Sec. | the exact and semiclassical expressions for the determinant of S-matrix for non- 
overlapping n-disk systems have been constructed. For the case of the 3-disk system they 
read 



det,S<3> W = (det l SM(to))' d f' M * ( *f» t (jW^ 
V ' \ 1 y ') detM Al {k) det,M A2 (A;) (detjM E (fc)) 



/ e -i7rAr(fc)\ 2x3 f -^l-disk(l)(fc*) -Z"l-disk(r)(fc*) 
\ ^l-disk(l) (k) ^l-disk(r) 



3 

X 



z Al (k*)* z M {k*y z E (k*f 

X ■ 



Z Al (k) Z A2 (k) z E (k) 



2 



(6.3) 



where the tilde indicates that diffractive corrections have to be included, in general. Espe- 
cially for the Ax-representation of the 3-disk system we therefore have the relation between 
the quantum-mechanical kernels and the Gutzwiller-Voros zeta functions 



det ; M Al (^) f sx . Z Al (k*)* 

det*M Al (A;) Z Al (k) ' 1 ' j 

where we have now neglected diffractive corrections. As argued in the conclusion section |f] 
both sides of Eq. and Eq. fl£T4|) respect unitarity; the quantum- mechanical side exactly 
and for the semiclassically side under the condition that the curvature expansion converges 
or that it is truncated. As all the n-disk resonances for non-overlapping n-disk repellers are 
below the real fc-axis, the border of absolute convergence, defined by the closest resonances 
to the real axis [T5|,^2j is inside the lower complex wave-number plane and unitarity on 
the real axis is guaranteed. Thus, if the wave number k is real, the left hand sides and also 
the right hand sides of eqs.( |6.3| ) and ( |6.4| ) can be written as exp{i2i](k)} with a real phase 
shift i](k). In fact, we can define a total phase shift for the coherent part of the 3-disk 
scattering problem (here always understood in the Ax-representation) for exact quantum 
mechanics as well as for the three semiclassical candidates by: 



- det,M(fc*)t 



(6.5) 



&etM(k) 

(6.6) 



,2ir, dyn (fc) ._ Co (k*) , _s 
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This phase shift definition should be compared with the cluster phase shift given in Sec. 4 
of Lloyd and Smith |53|]. For a separable system, as e.g. the 1-disk system (in the angular 
momentum representation), the cluster phase shift just corresponds to the sum 



V (k) = m{k) , (6.9) 

!=— oo 

as the S-matrix of the one-disk system (evaluated with respect to the center of the disk) 
reads 



-H?\ka ) 
H\ 1] {ka) 

e 2i ^5 u , , (6.10) 



such that 



+oo 

detS(fc) = n e2Mk) ■ ( 6 - n ) 

l=— oo 

Let us once more stress: the coherent or cluster phase shift is an experimentally accessible 
quantity: from the measured differential cross sections the elastic scattering amplitudes 
have to be constructed. This leads to the full phase shift of the 3-disk system including 
the contribution from the single disks. However, the incoherent part can be subtracted 
by either making reference experiments with just single disks at the same position where 
they used to be in the 3-disk problem or by numerical subtractions as the one-disk phase 
shifts are known analytically, since the system is separable, see ( |6.10|) and ( |6.11| ). In this 



way one can separate the incoherent phase shifts from the coherent ones. 

Thus i] qm (k) is "measurable" in principle. We next use these cluster phase shifts in order to 
discriminate between the various zeta functions. Below, we compare the exact quantum- 
mechanical cluster phase shift r] qm with 

(i) the semiclassical cluster phase shift t/gv(&) of the Gutzwiller-Voros zeta function 



(ii) with the semiclassical cluster phase shift ^dyn(^) of the dynamical zeta function (2.2), 

(iii) and with the semiclassical cluster phase shift rj qc i(k) of the quasiclassical zeta function 
f|27 



The zeta functions in the numerator as well as in the denominator of Z(k)*/Z(k) have 
been expanded to curvature order (=topological length) 12. For the Gutzwiller-Voros zeta 
function this is an overkill as already curvature order 4 should describe the data below 
Kek = 950/a. In fact, we have not seen any difference in the Gutzwiller-Voros calculation 
between the curvature order 3 and 12 results for k < 120/a and up to figure accuracy. 
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Curvature order 2, however, gives in the regime 100/a < k < 120/a noticeable deviations. 
On the other hand, as mentioned in Sec. |6.1| , the quasiclassical zeta function has problems 
for lower curvature orders with predicting the (sub-)leading resonances; therefore, these 
high curvature orders are used in order to give the quasiclassical zeta function as fair 
a chance as possible. The coherent phase shifts are compared in the window 104/ a < 
k < 109/a, which is a typical window narrow enough to resolve the rapid oscillations 
with k sufficiently large such that diffractive effects can be safely neglected. Furthermore, 
although we have no physical interpretation in terms of the S-matrix, we also compare 
in the same window the exact quantum-mechanical product detM(A;)detM(A;*)"'" with the 
squared modulus of the Gutzwiller-Voros zeta function Zgv(^)^gv(^*)*, the dynamical 
zeta function CcT^Ko" 1 ^*)*' an d the quasiclassical zeta function Z qc i(k)Z qc \(k*)* . Here 
k is taken to be real and the case of the 3-disk system in the Ai-representation with 
center-to-center separation R = 6a is studied. 



Consider finally the general quasiclassical zeta functions of Ref. 
ratio 



and especially the 



Z(k) :-. 



F+a,k)F-a,k) 



F_(l,k)F + (l,k) 
with F + (j3, k; z) and F_(j3, k; z) being defined as follows 



(6.12) 



OO 1 

F + (P,k;z)=exp{-Y,T,- 



>>H p {k)] 



[p] t p (k) 



,ri-/3+ 7 



00 1 A r \z 
F_((3,k-z)=exp{ - V E-rrr^ " 



Here the subleading factor (1 + | | 2/3 4 ) of Eq.(ll) in ref. Jy| has been removed as in 
Eq.(12) of ref. |43|]. When Eq.( |6.12| ) is used, the corresponding coherent phase shift 



z(k*y 



(6.13) 



works on the real wave-number axis and in the limit n —>■ oo (where n is the curvature 
order) as well as the original Gutzwiller-Voros zeta function. Hence, it does not matter 
here whether the Gutzwiller-Voros zeta function is directly expanded in the curvature ex- 
pansion or whether the individual determinants F + (|, fc), k), F_(|, fc) and F + (|, /c) 
are each expanded in separate curvature expansions up to the same curvature order and 
then inserted in the ratio ( |6.12| ). Note that the presence or absence of the subleading 
factor (1 + |Ap| 2/3 ~ 4 ) in the definitions of F + ((3, k; z) and F_(/3, k; z) does not change the 
results up to figure accuracy. 
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Fig. 6.2. (a) The coherent cluster phase shifts of the 3-disk scattering system in the 
Ai-representation with center-to-center separation R = 6a. The exact q uan tum-mechanical data 
are compared to the predictions of the Gutzwiller-Voros zeta function ( p.ip , the dynamical zeta 
function (2.2) and the quasiclassical zeta function ( |2.3| ) calculated up to 12th order in the cur- 
vature expansion, (b) The same for the squared moduli of the exact spectral determinant and 
the semiclassical zeta functions. The predictions of the Gutzwiller-Voros zeta function and exact 
quantum mechanics coincide within the resolutions of the plots. 

Let us stress that phase shifts are not only of theoretical interest, as are the sub leading 
resonances (which are completely shielded by the leading resonances), but hard data which 
can be extracted, in principle, from measured differential cross sections. 



In summary, even empirically, one can tell the three semiclassical zeta functions apart 
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and see which is the best. Again the Gutzwiller-Voros one — whether used directly or 
whether defined as the ratio ( |6.12| ) of four quasiclassical determinants as in ref.[|43| - is 
by far the best. 



6.3 The quantum-mechanical cumulant expansion versus the semiclassical curvature ex- 
pansion 



In this subsection it will be shown that the Gutzwiller-Voros zeta function approximates 
its quantum-mechanical counterpart, the characteristic KKR-type determinant |5l|j53|j49[| , 



only in an asymptotic sense, such that it always should be understood as a truncated series. 

As shown in Sec.|], the characteristic determinant and the Gutzwiller-Voros zeta-function 
are related as 



detM(fc) ^ Z GY {k) . (6.14) 

Let Qm{k) denote the m th cumulant of detM(fc) - i.e. the term proportional to z m in 
the Taylor expansion of det{l + zA(k)} - which satisfies the Plemelj-Smithies recursion 
relation ( ^4.6[ ) (see also App.[A|). Since the Plemelj-Smithies recursion formula is plagued 
by cancellations of very large numbers, we have not used the Plemelj-Smithies recursion 
relations for our numerical calculation of Q m (k), but instead we construct this quantity 
directly from the eigenvalues {Xj(k)} of the trace-class matrix A(k), i.e. 



Qm(k) 



(6.15) 



1<71<— <Jm<00 



(see again App.|A]for more details). Unfortunately, a semiclassical analog to this exact for- 
mula has not been found so far. Thus C m (k), the corresponding semiclassical mth-order 
curvature term, of ZQy(k), can only be constructed from the semiclassical equivalent of 
the Plemelj-Smithies recursion relation ( |4.8D which exactly corresponds to the standard 
curvature expansion of refs. pT9|,|42|JT5[1 ) and is therefore inherently plagued by large can- 



cellations. The cumulant and curvature expansions, truncated at nth order, read: 



detM(k)\ n =J2Qm(k) 

n 



(6.16) 
(6.17) 



m=0 



Let us recapitulate what we already know about these series. From Sec.f| together with 
the appendices [A| and |C| we deduce that the cumulant sum 

n 

lim detM(A;)| n = lim ^ Q m (k) = detM(fc) (6.18) 

m=0 
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is absolutely convergent, i.e. 



^ \Qm(k)\ < oo , (6.19) 

m=0 

because of the trace-class property of A(k) = M(/c) — 1 for non-overlapping, disconnected 
n-disk systems. On the other hand, as discussed in Refs.[|4T]-|4"5|], the Gutzwiller-Voros cur- 



vature sum converges only above an accumulation line (running below and approximately 



parallel to the real wave-number axis, see Sec. |6.1|) which is given by the first poles of 
the dynamical zeta function, CcT 1 ^)' or the leading zeros of the subleading zeta function. 
However, as shown in Sec. |6.1| , even below this boundary of convergence the truncated 
Gutzwiller-Voros curvature sum, Zgv(^)U approximates the quantum- mechanical data 
as an asymptotic series. 

In addition, a very important property for the discussion of the cumulant and curva- 
ture terms is the existence of the scaling formulas (established by us numerically) which 
relate the mth cumulants or curvatures inside the complex wave-number plane to the 
corresponding quantities on the real fc-axis: 



Q m (Re k + ilm k) ~ Q m (Re k) e - mLlmk (6.20) 
C m (Refc + iImfc)~C m (Refc)e- m£lmfc . (6.21) 

(For this to hold, diffractive effects have to be negligible, i.e. — Imk <C Refc.) Here L m 
R — 2a is the average of the geometrical lengths of the shortest periodic orbits, the two 
orbits of topological length one. The scaling can be motivated by the approximate relation 
Tr[A m (A;)] « {TrA(/c)} m which, of course, cannot be exact, as otherwise the cumulants 
would be identically zero. Nevertheless, the overall behaviour follows from this, since 



Tr [A (Re A; + ilm A;)] ~ Tr [A(Refc)] e" 



-Lira k 



From Fig. |6.3| one can deduce that the deviations between quantum-mechanical cumu- 
lants and semiclassical curvatures (as evaluated on the real fc-axis) decrease with increas- 
ing Re k, but increase with increasing curvature order m. The value of Re k where the 
quantum-mechanical and semiclassical curves join is approximately given by Re ka ~ 
2 m+1 . Approximately the same transition points can be generated from a comparison of 
the phases of the cumulant and curvatures. 



By varying the center-to-center distance we have numerically verified that the above 
limits generalize to the following relations valid on the real wave-number axis (k real and 
positive) : 



C m {k) « Q m (k) with 1> \C m {k)\n \Q m {k)\ if ka>2 m ~ 1 - (6.22) 

a 



and 
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Fig. 6.3. Comparison of the squared moduli (on a logarithmic scale) of the first seven quan- 
tum-mechanical cumulant terms, \Q n (k)\ 2 , with the corresponding semiclassical curvature terms, 
\C n (k)\ 2 , of the Gutzwiller-Voros zeta function (2.1) evaluated on the real wave-number axis k. 
The system is the Ai three-disk repeller with center-to-center separation R = 6a. 



l»|C m (A;)|»|Q m (A;)| if ka< 2 m ' 



(6.23) 



What is the interpretation of fl6.22|) and ( |6.23|) ? For fixed k, even in the regime, where 
ZGv{k)\ n converges, e.g., on the real fc-axis, the Gutzwiller-Voros zeta function is only an 
asymptotic approximation to the true quantum-mechanical cumulant sum, since for m > 
m cr it, defined by ka ~ 2 mcrit_1 -, the exact quantum- mechanical cumulants Q m (k) and the 
semiclassical curvatures C m (k) are grossly different. These deviations can be enhanced 
by the m th derivative, m > m cvlt) with respect to the book-keeping variable z, since this 
operation eliminates all approximately equal terms, such that the corresponding cumulant 
and curvature series are transformed to completely different expressions. The fact that 
ZQy(k)\ n - even in its convergence regime - is only an asymptotic expansion to the exact 
quantum mechanics is normally not visible, as the terms in ( |6.22|) are exponentially small 
on or close to the real axis and therefore sum to a tiny quantity. In other words, close to 
the real axis the absolute error \C n (k) — Q n {k) \ for m > m cr i t is still small. The relative 
error \C n (k) / Q n (k)\ on the other hand is tremendous (see Fig. |6.3|) . Deeper inside the 
negative complex wave-number plane, however, under the scaling rules fl6.20| ) and ( |6.21| ), 
the deviations Q6.23Q are blown up, such that the relative errors \C n {k)/Q n {k) \ eventually 
become visible as absolute errors \C n (k) — Q n (k)\ in the resonance calculation of Sec.|6.1|. 
If Imk is above the boundary line of convergence, these errors still sum up to a finite 
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quantity which might, however, not be negligible any longer, as was the case on or close 
to the real fc-axis. Below the convergence line these errors sum up to infinity. Thus the 
Gutzwiller-Voros curvature expansion Z G y(k)\ n does not suddenly become an asymptotic 
approximation to detM(/c)| n , it always is an asymptotic approximation (as shown by 
the relative error \C n (k)/Q n (k)\), even in its convergence regime above the accumulation 
line and even on the real axis, where the zeta function itself is in its domain of absolute 
convergence 
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Thus, the value of Imk where - - for a given m - - the Z G y(k)\ n sum deviates from 
detM(fc)| n is governed by the real part of k and the scaling rules ( |6.20|) and ( |6.21 



It has nothing to do with the boundary line between the convergence region and the 
asymptotic region of Z GV (fc), as the asymptotic expansion is given by a finite sum of 
all terms satisfying (|6.22|) . Therefore, the truncated Gutzwiller-Voros expansion describes 
the quantum-mechanical resonance data even below the line of convergence of the infinite 
curvature series, see |6.1| . On the other hand, the boundary line of the convergence regime 
of the Gutzwiller-Voros expansion is solely governed by those C m (k) which have nothing 
to do with the quantum analog Q m (k), i.e. solely by terms of character ( |6.23| ). The reason 
is, of course, that the convergence property of an infinite sum is governed by the infinite 
tail and not by the first few terms. Whether the Gutzwiller-Voros expansion converges or 
not is therefore not related to whether the quantum-mechanical data are described well 
or not. 

The convergence property of a semiclassical zeta function on the one hand and the approx- 
imate description of quantum mechanics by these zeta functions are two different things. 
It could happen that a zeta function is convergent, but not a good description of quantum 
mechanics (see especially the failure of the entire quasiclassical zeta function to approx- 



imate the exact cluster phase shift in Sec. 5.2). On the other hand it may not converge 



in general, but its finite truncations nevertheless approximate - at least to some order - 
quantum mechanics, as it is the case for the Gutzwiller-Voros zeta function. 

These findings hold for any re-writing of the Gutzwiller-Voros zeta function, as Z G y(k)\ n 
was already shown to be asymptotic in a regime where the curvature sum is still absolutely 
convergent and the limit lim^oo Z GV (fc)| n exists. Therefore, any re-writing of Z GV (/c), es- 
pecially the one of ref . [ 43| as the ratio of four quasiclassical zeta functions ( |6.12|) will at best 



work as an asymptotic expansion to the exact quantum-mechanical cumulant expansion. 
Note that, for finite curvature order n 
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(6.24) 



If the ratio is evaluated according to the r.h.s. of ( |6.24|) , one obtains exactly the same result 
as for the original Gutzwiller-Voros expansion using formula (|2.1| ). If, however, the ratio is 
evaluated according to l.h.s. of ( |6.24| ), the relation to the quantum-mechanical cumulant 
expansion is lost: the matching of the semiclassical coefficients of z m with the quantum- 
mechanical ones is spoiled, as the asymptotic terms, resulting from various curvature 
orders of the Z GV (fc)| n calculation, mix. If n is large enough, also the l.h.s. of ( |6.24|) 
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will deviate strongly from the quantum mechanics as the original formulation of the 
Gutzwiller- Voros expansion does — the difference is that this new expression approximates 
quantum mechanics at slower rate than the original formula, as the asymptotic terms of 
higher and lower curvature order are mixed. However, at high enough curvature order n 
also the new l.h.s. of ( 6.24 ) will encounter terms of class (|6.23| ) and will therefore — for 
large negative imaginary wave numbers - deviate strongly from the quantum-mechanical 
resonance data. 



What is the reason for the truncation at ka ~ 2 mcrit_1 L/a? This boundary follows from a 
combination of the uncertainty principle with ray optics and the exponentially increasing 
number of periodic orbits of the 3-disk repeller. For fixed wave number k, quantum me- 
chanics can only resolve the classical repelling set of the periodic orbits up to a critical 
topological order m cr i t . The quantum wave-packet which explores the repelling set, has 
to disentangle 2™ different sections of size d ~ a/2 n on the "visible" part of the disk 
surface between two successive collisions with the disk. Since these collisions are spatially 
separated by the mean length L, the flux spreads by a factor L/a. In other words, the 
non- vanishing value of the topological entropy for the 3-disk system, h ~ In 2, is the 
reason. For comparison, the uncertainty bound on the wave number in the hyperbolic, 
but non-chaotic two-disk system is independent of the curvature order (in case diffractive 
creeping is negligible), as there is only one geometrical periodic orbit and therefore the 
repelling set is trivial with zero topological entropy. 



The result that the semiclassical curvature expansion has to be truncated at finite or- 
der for a fixed wave number k, is different from the fact that the (in principle infinite) 
multiscattering kernel A mjm / = M mjm / — 5 m>TO / can be truncated to a finite matrix. The 
truncation in the curvature order is related to the resolution of the repelling set of pe- 
riodic orbits of the 3-disk system. The truncation in the size of the matrix is related to 
the semiclassical resolution of the single disks of the 3-disk system. The point particle 
classically only scatters from the disk, if its impact parameter is of the size or smaller 
than the disk radius a. Note that in the fundamental domain of the Ai disk system, one 
considers only one half-disk. Mathematically, this follows from the asymptotic behaviour 
of the ratio J m (ka) / H^\ka) which governs the scaling of the kernel A m m / and which is 
valid for m larger than m max , defined in Eq. (|6.2|), see App.|C|. In order to visualize this, 
we have plotted in Fig. OA the moduli of the eigenvalues (on a logarithmic scale and in 
descending order) of the multiscattering kernel Am^ik) of the Ai 3-disk repeller with 
R = 6a for the cases k = 100/a (on the real wave-number axis) and k = (100 — 1.25i)/a 
as function of the eigenvalue index j. The imaginary part of the latter wave number is 
characteristic for a domain where the subleading resonance bands emerge. The one-disk 
resolution is clearly visible in the exponentially decreasing tails of both curves above 
Refc ~ 140/a. In order to exhibit this feature, the matrix itself was truncated here at a 
large value of m = 220. Furthermore, from the curves, one can read off that only the first 
few eigenvalues (six for the upper curve corresponding to the case k = (100 — 1.25i)/a and 
two to four eigenvalues for the lower curves for the case k = 100/a) are "essential", i.e., 
are of the order unity or bigger. These numbers match very well the minimal topological 
order needed in the semiclassical calculation to approximate the relevant resonances at 
the specified k values. Whereas inside the negative complex /c-plane one has to go to 
higher curvature orders in the truncation of the semiclassical zeta function in order to 
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find all the subleading resonances (namely to order six for the specified k value), on or 
close to the real axis only the leading resonances are "visible" , in agreement with the data 
for the cluster phase shifts which, for the specified k value, can be well approximate by a 
semiclassical calculation of order three to four. One can also extract from the figure what 
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Fig. 6.4. The moduli of the eigenvalues of the multiscattering kernel -A m m /(/c) of the Ai 3-disk 
repeller with R = 6a for the cases k = 100/a and k = (100 — 1.25i)/a on a logarithmic scale. 
The eigenvalues are displayed as function of their index j in descending order. 

happens if the negative imaginary part of k is increased: the curve is basically parallelly 
shifted upwards. Thus the number of eigenvalues and the minimal curvature order for the 
semiclassical description of quantum mechanics increases the deeper one "dives" into the 
lower complex wave-number plane. Although only a few eigenvalues are essential for the 
computation of the resonances and phase shifts, the size of the matrix is determined by 
the much bigger number m max « (e/2)ka, such that many more eigenvalues are produced 
by a matrix diagonalization code. Unfortunately, one cannot escape this mismatch, as the 
model space for the matrix M has to be that large in order to guarantee stable numerical 
results for the leading eigenvalues. 

In summary, the minimal size of the matrix is determined by the resolution of the single 
disks, whereas the maximal topological order up to which the semiclassical curvature 
expansion makes sense, follows from the Heisenberg uncertainty limit on the quantum 
resolution of the repelling set. The topological exponential rise of the number of periodic 
orbits, with increasing curvature expansion order n, is the physical reason for the eventual 
breakdown of the curvature expansion of the semiclassical zeta function ( |2.1| ) as compared 
with the exact quantum-mechanical cumulant expansion which defines the determinant 
of the multiscattering matrix in an infinite-dimensional Hilbert space. 
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7 Conclusions 



Starting from the exact quantum-mechanical S-matrix we have tried to find a direct 
derivation of the semiclassical spectral function for a rather special class of classically 
hyperbolic scattering systems, namely the non-overlapping disconnected finite n-disk re- 
pellers in two dimensions. We have confined our investigation to these systems as they 
are on the one hand "realistic" enough to capture the essence of classically hyperbolic 
scattering problems (or, for certain geometries, even chaotic scattering problems) and on 
the other hand simple enough to allow for a "top-down" approach from exact quantum 
mechanics to semiclassics without, and this is the important point, any "formal" step 



in between. We have reason for this "pedantry" : It is known from the work of Refs. j|I 



43] that the standard spectral function, the cycle-expanded semiclassical zeta function of 
Gutzwiller and Voros, is not entire for the 3-disk system and therefore fails to describe 
subleading scattering resonances in the complex wavenumber plane below its boundary of 
convergence. The question is whether these failures are induced by unjustified formal steps 
in the semiclassical reduction or whether they are inherently a property of the semiclas- 
sics itself. Since we expected that the semiclassical spectral function must follow from the 
semiclassical reduction of the cumulant expansion of the corresponding multiscattering 



kernel [44,45 , we had to avoid any bias or unjustified assumption on the exact kernel as 



well as on the semiclassical spectral function, e.g., on the existence and regularization of 



the quantum mechanical expression , on the structure of the period orbits, especially 
on the structure of their stabilization, on the organization of the spectral function in terms 
of cycles or curvatures, on the classification of these curvatures by the topological lengths 
of the orbits, etc. 

Our first task therefore has been to ensure that the quantum mechanical starting point 
for the semiclassical reduction is well-defined. The T-matrix of the n-disk scattering sys- 
tems, derived by the methods of stationary scattering theory |T3" |, was shown to ex- 



ist on the real /c-axis and to be trace-class. Therefore the actual starting point for the 
spectral classification of the scattering system, the determinant of the S-matrix, exists 
also and can be manipulated by cyclic permutations, unitary transformations, splitting 
into sub-determinants, in other words, operations which are non-trivial for matrices of 
infinite-dimensional Hilbert spaces. With the help of these (then justified) operations, we 
succeeded in transforming the determinant of the S-matrix into a form that is well suited 
for the semiclassical reduction step, see Eq. (|3.9j). It separates into the incoherent superpo- 
sition of n one-disk scattering determinants and into the ratio of the determinant and the 
complex conjugated determinant of the genuine multiscattering matrix, M. Furthermore, 
the determinant of the multiscattering matrix can be decomposed into sub-determinants, 
if the n-disk system has additional symmetries. All of the above mentioned determinants 
are shown to exist separately. This is one of the key points for the semiclassical reduc- 
tion, since the existence of the S-matrix alone would not guarantee that the one-disk 
aspects can be separated from the multi-disk aspects in a well-defined manner. Note that 
the standard geometrical periodic orbits (without creeping) can only "know" about the 



# 3 In fact, even the quantum- mechanical expression is not entire in the whole complex k- 
plane, since it has a branch cut on the negative real axis and poles which cancel the one-disk 
singularities. 
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multi-disk aspects, and not about the single disk aspects. 

As the determinants are taken over infinite dimensional matrices, one has to worry about 
their very definition. The von-Koch criterion (the existence of the determinant in one 
orthonormal basis, see App. [A.3[ ) is not sufficient for this task, since implicitly in the 
derivation of the S-matrix and explicitly under symmetry-reductions unitary transfor- 
mations are mandatory. The multiscattering matrix must be reducible to a form "unit 
matrix plus a trace-class matrix" in order for its determinant to exist. Fortunately, we 
could prove that the multiscattering kernel A = M — 1 is trace-class for any n-disk 
geometry as long as the disks do not overlap nor touch each other. Furthermore, the 
determinant over the infinite matrix M is defined as a cumulant expansion which - as 
shown by us - semiclassically reduces to the curvature expansion. Thus already quantum 
mechanically the cumulant/curvature "regularization" emerges. Moreover, by working in 
the full domain of the n-disk system, we could show in Sec.[S] that the cumulants split via 
the quantum traces into quantum itineraries which can be classified by the very symbol 
dynamics of their semiclassical reductions - the semiclassical periodic orbits. Thus the cu- 
mulant/curvature ordering in terms of the topological lengths of the quantum itineraries 
and hence of the topological lengths of the semiclassical periodic orbits is already present 
on the quantum-mechanical level. One does not need to impose it from the outside (as 
it would be the case for the semiclassical reduction of the Krein-Friedel-Lloyd sums of 
two bounded reference systems) ; but it follows naturally from the properties of quantum 
mechanics; namely, from the defining cumulant expansion of the determinant of the exact 
multi-scattering matrix. 

Thus the classification by the quantum itineraries is a virtue, but it is unfortunately 
also a vice, as quantum traces and the Plemelj-Smithies recursion formulas are involved. 
The latter introduce (unnecessarily) large terms which finally cancel in the cumulants 
themselves. As the cumulant sum of a trace-class operator converges absolutely, a direct 
semiclassical reduction of a complete cumulant and not of the potentially large quan- 
tum itineraries or quantum traces would be highly desirable. Unfortunately, the direct 
semiclassical reduction of a complete cumulant is not known. It might correspond to an 
integration of the small differences between the direct motion and the shadowed motion of 
the quantum wave packet. Instead, the standard calculation for the complete curvatures 
proceeds through the shadowing of all full periodic orbits of a the pertinent topologi- 
cal length by all pseudo-orbits (=products of shorter periodic orbits) of the same total 
topological length. Because of these large cancellations, the semiclassical reduction on 
the level of the itineraries is potentially dangerous for the semiclassical equivalent of the 
quantum-mechanically absolutely converging cumulant sum, the curvature sum. There is 
no guarantee that it converges as well. 

As mentioned, in Sec.[| we have managed to construct the semiclassical equivalent for each 
specified quantum itinerary. By working in the full domain and utilizing the pertinent 
simple symbol dynamics, which is valid under the condition that the number of disks is 
finite and that the disks do not overlap nor touch, our semiclassical reduction applies 
for all n-disk geometries, with one exception: we have to veto geometries which allow for 
grazing periodic orbits. In this way, we could guarantee for any specified quantum itinerary 
that the sequence of disk labels transforms uniquely to a sequence of non-overlapping 
semiclassical saddles in the complex angular-momentum plane which corresponds to the 
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standard semiclassical periodic orbit, specified by the same symbolic sequence. The weight 
of the periodic orbits was shown to be identical to the one derived by Gutzwiller ||. 
Furthermore, we have shown that to each itinerary that generates a (non-creeping or 
creeping) periodic orbit with a "ghost tunneling" section straight through a disk there 
belong "parent" itineraries, such that the ghost and corresponding parent periodic orbits 
cancel exactly in the semiclassical curvature sum. This establishes how pruning emerges 
from quantum mechanics in the semiclassical reduction. We have also shown that, to each 
quantum itinerary of topological length n, there belong 3 n — 1 different periodic orbits 
which contain creeping sections and we have specified a generalization of the symbolic 
labelling. By the Watson contour method of Ref. |29fl we have derived their structure 



which agrees with the result of the semiclassical construction of Refs. ]3"2|,|3"1|] which in 



turn is based on Keller's semiclassical theory of diffraction ||30|| . The direct link of the 
determinant of the exact S-matrix (via the determinant of the multiscattering matrix, via 
its cumulants and quantum itineraries) to the periodic orbits is therefore established. If 



the operations are inverted, the right hand side of Eq. (|4.12|) emerges, modulo the caveat 



that the semiclassical curvature sum might not converge, in general. 

What is known about the convergence properties of the curvature expansion of the 
Gutzwiller- Voros zeta function from the literature? The Gutzwiller trace as well as the 
zeta-function for n-disk repellers is known to converge (even absolutely) in the complex 
wave-number plane above a line specified by the resonance with largest imaginary k- value, 
see e.g. Ref. [|T7| . As all resonances belong to the lower half of the complex wave- number 
plane, the zeta function converges at least on the real /c-axis. From Refs. PT^[43f it is known 
that the cycle or curvature expansion of the Gutzwiller- Voros sum converges even inside 
the resonance region above an accumulation line defined by the poles of the dynamical 
zeta function. Thus above this accumulation line (and away from the branch cut and sin- 
gularities of the exact quantum-mechanical side) our semiclassical limit ( |4 . X 2| ) (or ( |4.13| ) 



for symmetry-reducible problems) is established for the full, untruncated Gutzwiller- Voros 
zeta function; see, however, below for the discussion of the asymptotic behaviour of the 
curvature sum. 

The relations are compatible with Berry's expression for the integrated spectral density 



in Sinai's billiard (a bounded n — ► oo disk system, see Eq.(6.11) of Ref. [49|) and - in 
general - with the Krein-Friedel-Lloyd sums fl2.6| ). They justify the formal manipulations 
of Refs. [|l7||ll^|7ll . Furthermore, for these scattering systems, unitarity is automatically 
preserved semiclassically (without reference of any re-summation techniques a la Berry 
and Keating|21| which are needed in bounded problems). Quantummechanically, unitarity 



follows from the relation 



detS (n) (A;) t 



detS( ft )(Jfe*) 



(7.1) 



which is manifestly the case (see the first lines of Eq. Q4.12|) and Semiclassically, 
this follows from the second lines of ( |4.12| ) and ( |4 . 1 3| ) , with the caveat that curvature sums 
on the right hand sides must exist, i.e., they either converge or are suitably truncated. This 
is of course a very pleasant property. But, on the other hand, unitarity can therefore not 
be used in scattering problems to gain any constraints on the structure of Z G y, as it could 
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in bounded problems, see pi| . Why are bounded problems special? In the semiclassical 
treatment of scattering problems the poles of the determinant of the S-matrix result from 
the zeros of Zcy(k) in the lower complex /c-plane (where in general - except at the zeros 
- ZQy(k) dominates ZQy(k*)* which is small, but nonzero), whereas the zeros of the 
determinant of the S-matrix are produced by the ones of ZQy(k*)* i n the upper complex 
£;-plane (where in turn Zqv(^) is the small, but nonzero zeta- function). For bounded 
problems k is real and both zeta- functions become equally important. (A sign of this is 
the fact that the Hankel functions of either first or second kind which appear in the spectral 
determinants are replaced by the corresponding Bessel functions.) This obviously calls for 
a fine-tuning, hence, the re-summation. Note also the symmetry-breaking ie prescription 
which had to be added to the l.h.s. of the Krein-Friedel-Lloyd sums, see Sec.||. 



As stated above, the incoherent single-disk scattering decouples from the genuine multi- 
disk scattering. The 1-disk poles do not influence the position of the genuine multi-disk 
poles. However, DetM(/c) does not only possess zeros, but also poles. The latter exactly 
cancel the poles of the product over the 1-disk determinants, ITj=i detS^ (kaj), since both 
involve the same "number" and "power" of H^(kaj) Hankel functions in the denominator. 
The same is true for the poles of DetM(fc*)^ and the zeros of IT/=i detS^ (kaj) , as in this 
case the "number" of H^(kaj) Hankel functions in the denominator of the former and the 
numerator of the latter is the same — see also Berry's discussion of the same cancelation 
in the integrated spectral density of Sinai's billiard (Eq.(6.10) of Ref . f49|| ) . Semiclassi- 
cally, this cancelation corresponds to a removal of the additional creeping contributions 
of topological length zero, i.e., 1/(1 — exp(i2-7n^)), from Zqv via the 1-disk diffractive 
zeta functions, Zdisk(ji) and -Zdisk(jr)- The orbits of topological length zero result from the 
geometrical sums over additional creepings around the single disks, X^=o( exp(i27ri^) ) n ' w 
(see [0]), and multiply the ordinary creeping paths which are classified by their topo- 
logical length. Their cancellation is very important for situations where the disks nearly 
touch, as in such cases the full circulations by creeping orbits of any of the touching disks 
should clearly be suppressed, as it now is. Therefore, it is important to keep a consistent 
count of the diffractive contributions in the semiclassical reduction. 



What happens to the resonances, when the spacing between the disks becomes vanishing 
small such that bounded regions are formed in the limit of n > 2 touching disks? Because 
of the ratio DetM^^/DetM^), to each (quantum-mechanical or semiclassical pole) of 
detS^^/c) in the lower complex fc-plane there belongs a zero of detS(fc) in the upper com- 
plex fc-plane with the same Re k value, but opposite Im k. When the bounded regions are 
formed some of these opposite zeros/poles move onto the real axis (such that their con- 
tributions cancel out of ( |3.9|) ). We have convinced ourselves that for the 3-disk scattering 
system with e > separation these resonances approach infinitesimally the bound-state 
eigenvalues of the complementary calculation of the spectrum inside the bounded region, 
see, e.g., Ref. [[P7H for the billiard bounded by three touching disks. Semiclassically, this 
would be a non-trivial calculation as the eigen-energies have to be real which — without 
resummation a la Berry and Keating ||2lJ - they are not. In this situation, one really 
has to think about further resummation techniques. Most of the resonances, however, do 
not move onto the real axis at all, as n-disk repellers, even with bounded sub-domains, 
are still scattering systems. The would-be bound states, however, drop out of the exact 
formula for detS^(fc), as they should. 
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Let us come back to the numerical data of Sec.|6] and the existence of the curvature 
expansion. In this section we have reported on numerical results for the exact quantum- 
mechanical Ax resonances of the three-disk system with R = 6a in the complex /c-plane 
in the region: < Re A; < 250/a and > Imk > — 1.6/ a. The first observation is that 
the quantum-mechanical resonances in this window can be grouped into (leading and 
subleading) bands. In addition to the data presented in Sec.|6.1|, where we have related the 
band structure to the semiclassical curvature expansion, it has been numerically checked 
that the emergence of a new band is in fact linked to a new cumulant order. The data of this 
window up to Re k = 200/a can be fitted very well with a quantum- mechanical cumulant 
expansion which is truncated at order seven. This knowledge, together with the fact that 
any periodic orbit results from the semiclassical reduction of a quantum itinerary with 
the same symbol sequence (in the full domain), tells us that periodic orbits of topological 
length eight and higher are completely irrelevant for the description of the presented 
quantum-mechanical data, for regions below Refc ~ 200/a and above Imk = — 1.6 /a. 
Thus any deviation of semiclassical predictions from the exact data cannot be cured 
by the inclusion of higher topological orbits. At best, they should leave the resonances 
untouched. 



This finding seems to be at variance with the result of Sec. 6.1 where the quasiclassical zeta 



function of Ref. |43| approximates most of the exact resonances at curvature order twelve. 



However, this truncated zeta functions finds also six erroneous resonance bands which do 
not have quantum-mechanical counter parts. This means that its topological expansion 
does not match the cumulant expansion, as we know of course from the analytical results 
of Sec||. The semiclassical reduction of a cumulant sum is the Gutzwiller-Voros curvature 
sum, and not the cycle expansion of the quasiclassical zeta function nor of the dynamical 
one. In fact, as shown by the comparison of exact to semiclassical coherent phase shifts, 
see Sec. |6.2| , the latter two zeta functions are very inferior to the Gutzwiller-Voros zeta 
function which describes the exact phase shifts up to the resolution of the plot. Any 
competitor zeta function should do at least as well as the Gutzwiller-Voros one in order 
to be taken seriously. The question whether it converges or not is not a criterion for how 
successfully it approximates the quantum-mechanical data. 



In Sec.|6T3| we have finally executed what was already advocated by us in Refs. . For 

the 3-disk example we have numerically compared term by term and order by order the 
quantum-mechanical cumulants with their semiclassical counterparts, the curvatures of 
the Gutzwiller-Voros zeta function. The numerical data show that, for a fixed value of the 
wave number, the cumulants and curvatures agree in magnitude and also in phase only 
up to a finite cumulant order which is determined by the wave number and then deviate 
strongly. We have interpreted this result from the uncertainty bound on the quantum- 
mechanical resolution of the details of the classically repelling set which exponentially grow 
with the topological entropy of the system under consideration. Close to the real axis these 
deviations are hidden by the smallness in absolute value of the higher-order cumulants and 
curvatures. Therefore, the semiclassical phase shifts agree very well with the quantum- 
mechanical ones, even for very high curvature orders. However, with increasing value for 
— Im k, inside the lower half of the complex plane, the deviations are enhanced by the 
scaling laws discussed in Sec.EO such that they eventually become noticeable. 



This observation matches extremely well the results of the resonance comparison in 
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Sec. |6.1| . The resonances which are located above and away from the boundary of con- 
vergence are approximated by the Gutzwiller-Voros curvature expansion as soon as the 
curvature order is sufficiently high. The resonances at or below the boundary of con- 
vergence, however, are approximated only up to the curvature order which respects the 
uncertainty bound. The curvature expansion works there only as an asymptotic series. 

Our interpretation is that eventually quantum mechanics and (semi-) classics have to part 
ways, as the quantum-mechanical spectral data only need power-law complexity, i.e. N 3 
operations if the multiscattering matrix can be truncated as an iV x N matrix, whereas the 
resolution of the classically repelling set needs exponential complexity if the topological 
entropy is non-zero. In other words, whether the curvature expansion converges or not with 
respect to quantum mechanics it should be truncated at the cumulant order specified by 
the uncertainty bound. All curvature terms exceeding this order are - from the quantum- 
mechanical point of view - irrelevant. From this perspective, the semi classical side of 
the relation (|4.12|) (and ( |4.13| )) should be interpreted to be valid just for the truncated 
Gutzwiller-Voros curvature sums, where the order of the truncation increases with increas- 
ing value of Re k (or, since k = p/h, with decreasing h). The semiclassical limit Re k —>■ oo 
and the cumulant limit m — > oo do not commute, in general, if the topological entropy 
is non-zero. These facts should be kept separated from the ^-effects of Refs. [[^-|41J which 
investigate the 0(h) corrections to the periodic orbits. We discuss here the ^-corrections 
to the curvatures which result from the periodic orbits via large cancellations against the 
pseudo-orbits. Part of the ^-corrections of Refs. [j38| -|4"C| cancel out as well, as can be shown 



from the comparison of the difference between the m-th order exact and semiclassical trace 
which exceeds by far in magnitude the difference between the corresponding m-th order 
cumulant and curvature. In fact, from the discussion of the subleading Debye corrections 
in App. |F.3| one can deduce that each term H\ 2 \ka) / H\ l \ka) introduces a correction fac- 
tor of order (1 + ih/4pa), such that the quantum itinerary of topological length m has 
at least an 0(1 + mih/Apa) factor relative to the corresponding periodic orbit (assuming 
that all disks have the same radius a for simplicity). However, the pseudo-itineraries of 
order m (which are the quantum mechanical analog of the pseudo-orbits) have the same 
correction factor, such that it cancels in the corresponding cumulant. Thus, the 0(mh) 
terms cancel. But what about the 0{mh ) terms which might be of the same importance 
as the 0(h) terms, as the limits m — > oo and h — > do not commute? Only if the cu- 
mulant sum is truncated at a finite order, the 0(mh 2 ) terms become negligible relative 
to the 0(h) terms. In principle, the uncertainty boundary should be derivable from the 
semiclassical reductions of Sec.|5] and App.[F] to the quantum itineraries, if ^-corrections 
are taken into account systematically. In practice, however, there is a very long way from 
the /i-corrections extracted from the quantum itineraries to the surviving ^-corrections 
on the cumulant level because of the very large cancellations of the quantum itineraries 
with the pseudo-itineraries. 
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A Traces and determinants of infinite dimensional matrices 



This Appendix summarizes the definitions and properties for trace-class and Hilbert- 
Schmidt matrices and operators, the determinants over infinite dimensional matrices and 
possible regularization schemes for matrices or operators which are not of trace-class. 



A.l Trace class and Hilbert- Schmidt class 



This section is based on Ref. 721 and also Refs. 60,73-751 which should be consulted for 



further details and proofs. The trace class and Hilbert-Schmidt property will be defined 
here for linear, in general nonhermitean operators A G C{7i): 7i — > 7i (where 7i is a 
separable Hilbert space). The transcription to matrix elements (used in the prior chapters) 
is simply aij = ((pi, A<pj) where {0 n } is an orthonormal basis of 7i and ( , ) is the inner 
product in 7i (see Ref. ||74j| where the theory of von Koch matrices of Ref. [76] is discussed.). 
Thus the trace is the generalization of the usual notion of the sum of the diagonal elements 
of a matrix; but because infinite sums are involved, not all operators will have a trace 
and, if the trace exists in one basis, it is nontrivial that it exists also in any other basis: 

(A) An operator A is called trace-class, A G J7i, if and only if, for every orthonormal 
basis, {4> n }: 



l(0n, A0 n )| < oo 



(A.l) 



The family of all trace-class operators is denoted by J\. 
(B) An operator A is called Hilbert-Schmidt, A G 

orthonormal basis, {<f) n }'- 



EII A ^| 



< oo 



J7~2, if and only if, for every 



(A.2) 



The family of all Hilbert-Schmidt operators is denoted by J 2 - 
(C) Bounded operators B are dual to trace-class operators. They satisfy the the 
following condition: |(0,B</>)| < C||^||||0|| with C < oo and if),<f> G H. If they have 
eigenvalues, these are bounded as well. The family of bounded operators is denoted by 



B(H) with the norm ||B|| = sup < ^ '"j^fp' for </> G 7i. Examples for bounded operators 
are unitary operators and especially the unit matrix. In fact, every bounded operator 
can be written as a linear combination of four unitary operators |6(| . 
(D) An operator A is called positive, A > 0, if (A0, cp) > V0 G H. Notice that 
At A > 0. We define |A| = v / AtA. 



The most important properties of the trace and Hilbert-Schmidt classes can be summa- 
rized as (see Refs.[ |60| , |72| ): 



(a) Ji and J 2 are *ideals., i.e., they are vector spaces closed under scalar multiplication, 
sums, adjoints, and multiplication with bounded operators. 

(b) A G Ji if and only if A = BC with B, C G J 2 . 
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(c) Ji C J 2 . 

(d) For any operator A, we have A G J 2 if ||A0 n || 2 < oo for a single basis. 
For any operator A > 0, we have A G J7i if | (<f) n , A0 n ) | < oo for a single basis. 

(e) If A G Ji, Tr(A) = 2(0 n , A0 n ) is independent of the basis used. 

(f) Tr is linear and obeys Tr(A f ) = Tr(A); Tr(AB) = Tr(BA) if either A G Ji and B 
bounded, A bounded and B G J\ or both A, B G J 2 - 

Note that the most important property for proving that an operator is trace-class is the 
decomposition (b) into two Hilbert-Schmidt ones, as the Hilbert-Schmidt property can be 
easily verified in one single orthonormal basis (see (d)). Property (e) ensures then that 
the trace is the same in any basis. Properties (a) and (f) show that trace-class operators 
behave in complete analogy to finite-rank operators. The proof whether a matrix is trace- 
class (or Hilbert-Schmidt) or not simplifies enormously for diagonal matrices, as then the 
second part of property (d) is directly applicable: just the moduli of the eigenvalues (or - 
in case of Hilbert-Schmidt - the absolute squares) have to be summed in order to answer 
that question. A good strategy for checking the trace-class character of a general matrix 
A is therefore the decomposition into two matrices B and C where one, say C, should be 
chosen to be diagonal and either just barely of Hilbert-Schmidt character leaving enough 
freedom for its partner B or of trace-class character such that one only has to show the 
boundedness for B. 



A. 2 Determinants det(l+A) of trace- class operators A 



This section is mainly based on Refs. |6T| , |73| which should be consulted for further details 
and proofs. See also Refs. |71| , [7o] | . 

Pre-definitions (Alternating algebra and Fock spaces): 

Given a Hilbert space 7i, ® n 7i is defined as the vector space of multilinear functionals on 
H with 0i <g> ■ ■ ■ <g> n G ® n H if 0i, . . . , 0„ G H. A n (H) is defined as the subspace of ® n U 
spanned by the wedge-product 



il A • • • A 4> n = —j= e ( 7r )K(l) ® • • • ® &r(n)] > ( A - 3 ) 



where V n is the group of all permutations of n letters and e(n) = ±1 depending on 
whether n is an even or odd permutation. The inner product in f\ n (7i) is given by 



(01 A • • • A <j> n , r]i A • • • A T] n ) = det {(0,, t^)} , (A.4) 

where det{a^} = J2irev n e ( 7T ) a i7r(i) ' " •°7wr(n)- A n (A) is defined as functor ( a functor sat- 
isfies A n (AB) = A n (A) A n (B) ) on f\ n {H) with 



/\"(A) (0i A • • • A 0„) = A0i A • • • A A0 n . (A.5) 
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Properties: If A trace-class, i.e., A G Ji, then for any positive integer k, /\ k (A) is 
trace-class, and for any orthonormal basis {</> n } the cumulant 

Tr(A\A))= £ ((^A-.-A^MA^A-.-AA^J) (A.6) 

V 7 h<-<i k 

is finite and independent of the basis. TrA°(A) = 1. 
Definition: Let A G Ji, then det(l + A) is defined as 

det(l + A) = Xl Tr (A (A)) • (A.7) 

k=0 



Properties: Let A be a linear operator on a separable Hilbert space H and {<fij}f an 
orthonormal basis. 



( a ) Efclo Tr (A fe (A)j converges for each A G Ji. 

(b) |det(l + A) I < rijli (1 + A*j(A)) where l^j(A) are the singular values of A, i.e., the 
eigenvalues of |A| = VaTK, and |det(l + A)| < exp(Tr|A|). 

(c) For any Ai, . . . , A n G Ji, (zi, . . . , z n ) 1— > det (1 + E"=i ^Aj) is an entire analytic 
function. 

(d) If A,B G Ji, then 

det(l + A)det(l + B) = det (1 + A + B + AB) 

= det((l + A)(l + B)) = det((l + B)(l + A)) . (A.8) 

If A G Ji and U unitary, then 

det (U f (l + A)U) = det (l + U f AU) = det(l + A) . (A.9) 

(e) If A G Ji, then (1 + A) is invertible if and only if det(l + A) ^ 0. 

(f ) If A 7^ is an n-times degenerate eigenvalue of A G Ji, then det(l + zA) has a zero 
of order n at z — — 1/A. 

(g) For any A G Ji, 

N(A) 

det(l + A)= [] (1 + ^(A)) , (A.10) 
i=i 

where here and in the following {Aj(A)}^i ^ are the eigenvalues of A counted with 
algebraic multiplicity (N(A) can of course be infinite). 

(h) If A G Ji, then 

iv(A fc (A)) 

TtU\a)) = £ A,(aW)= £ A J1 (A)...A Jfc (A)<oo. 

i =1 i<ji<-<j k <N(A) 

(AH) 
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(i) If A G Ji, then 



det (1 + zA) = Y,z k £ K ( A) • • • X Jk (A) < oo . 

fc=0 l<h<-<j h <N{A) 

(j) If A G Ji, then for \z\ small (i.e., \z\ max|Aj(A)| < 1), the series J2kLi z fc Tr 
converges and 



(oo k \ 
-E T Tr ((- A ) fc )j 



= exp(Trln(l + ^A)) . 
(k) The Plemelj- Smithies formula: Define a m (A) for A G J\ by 



det(l + 2A) = ^ 



m=0 



m! 



Then a m (A) is given by the m x m determinant 



a m (A) 



Tr(A) m — 1 
Tr(A 2 ) Tr(A) 
Tr(A 3 ) Tr(A 2 ) 





m — 2 
Tr(A) 









Tr(A m ) Tr(A (m_1) ) Tr(A( m " 2 )) ••• Tr(A) 



(A.12) 
-A) k ) /k 



(A.13) 



(A. 14) 



(A.15) 



with the understanding that a (A) = 1 and «i(A) = Tr(A). Thus the cumulants 
Q m (A) = a m (A)/m\ (with Q (A) = 1) satisfy the recursion relation 



Q m (A) = -^(-l) fe+1 g m _ fc (A)Tr(A fe ) for m > 1 



(A.16) 



fc=i 



Note that formula (|A.14j ) is the quantum analog of the curvature expansion of the 
Gutzwiller-Voros zeta function with Tr(A m ) corresponding to the sum of all periodic 
orbits (primitive and also repeated ones) of total topological length m, see Eq. ( |4.8| ). In 
fact, in the cumulant expansion (|A.14|) (as well as in the curvature expansion there are 
large cancellations involved: Let us order - without loss of generality - the eigenvalues of 
the operator A G J\ as: 



I Ai| > |A 2 | > • • • > |Ai_i| > |Ai| > |A i+1 | > • • • . 

This is always possible because of Ya=i ^ < 00 • Then, in the standard (Plemelj- 
Smithies) cumulant evaluation of the determinant, Eq.( |A.14] ), there are enormous can- 
cellations of large numbers, e.g., at the k th cumulant order (k > 3), all the intrinsically 
large "numbers" , X^ 1 \2, ■ ■ ., A^~ 2 A2A3, . . . and many more have to cancel out exactly 
until the r.h.s. of ( |A.11| ) is finally left over. Algebraically, the fact that there are these 
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large cancellations is of course of no importance. However, if the determinant is calculated 
numerically, the large cancellations might spoil the result or even the convergence. 



A. 3 Von Koch matrices 



Implicitly, many of the above properties are based on the theory of von Koch matri- 
ces [ |71 , [7B ,0 : An infinite matrix 1 — A = \\Sjk — ajk\\f, consisting of complex numbers, is 
called a matrix with an absolutely convergent determinant, if the series \ a nk^nk 2 ' ' ' a j n ,k, 
with n = 1,2, •• • converges, where the sum extends over all pairs of systems of indices 
(ji, j'2, ■ • • , j n ) and (hi, &2, - • • , k n ) which differ from each other only by a permutation, 
and where j\ < 32 < • • • jn- Then the limit 



lim det||&;t — a,-i.||? = det(l — A) 

n— >oo J 

exists and is called the determinant of the matrix 1 — A. The matrix 1 — A is called von 
Koch matrix, if both conditions 



J2\ajj\< 00, (A. 17) 

i=i 

00 

J2 M 2 <oo (A.18) 
j,k=i 



are fulfilled. Then the following holds (see Ref . [ 74) , TJ\ ) : 



(a) Every von Koch matrix has an absolutely convergent determinant. If the elements 
of a von Koch matrix are functions of some parameter \x (a^ = a 3 fe(/i), j, k = 1, 2, • • •) 
and both series in the defining conditions, ( A. 17 ) and ( A.18 ), converge uniformly in the 



domain of the parameter /1, then as n — » 00 the determinant det||5jfc — ctjki^Wi t en ds 
to the determinant det(l + A(/i)) uniformly with respect to fi, over the domain of \l. 
(b) If the matrices 1 — A and 1 — B are von Koch matrices, then their product 1 — C = 
(1 — A)(l — B) is a von Koch matrix, and det(l — C) = det(l — A) det(l — B). 

Note that every trace-class matrix A e J\ is also a von Koch matrix (and that any matrix 
satisfying condition (|A.18| ) is Hilbert-Schmidt and vice versa). The inverse implication, 
however, is not true: von Koch matrices are not automatically trace-class. The caveat 
is that the definition of von Koch matrices is basis-dependent, whereas the trace-class 
property is basis- independent. As the traces involve infinite sums, the basis-independence 
is not at all trivial. An example for an infinite matrix which is von Koch, but not trace- 
class is the following: 



2/j for i — j = — 1 and j even , 
Ay = ^ 2/i for i— j = +1 and i even , (A. 19) 

else , 
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i.e., 



A 



/Ol 

10 

1/2 

1/2 

1/3 

1/3 



(A.20) 



Obviously, condition (|A.17|) is fulfilled by definition. Secondly, condition (|A.18|) is satisfied 
as Y£?=i 2/n 2 < oo. However, the sum over the moduli of the eigenvalues is just twice the 
harmonic series J2™=i l/ ra which does not converge. The matrix (|A.20| ) violates the trace- 



class definition (|A.l ), as in its eigenbasis the sum over the moduli of its diagonal elements 
is infinite. Thus the absolute convergence is traded for a conditional convergence, since the 
sum over the eigenvalues themselves can be arranged to still be zero, if the eigenvalues 
with the same modulus are summed first. Absolute convergence is of course essential, if 
sums have to be rearranged or exchanged. Thus, the trace-class property is indispensable 
for any controlled unitary transformation of an infinite determinant, as then there will be 
necessarily a change of basis and in general also a re-ordering of the corresponding traces. 



Nevertheless, the von-Koch-criteria ( A. 17 ) and ( A.18|) are useful, as any trace-class matrix 
has at least to meet these simple tests which can be easily performed in any specified basis. 



A. 4 Regularization 



Many interesting operators are not of trace-class (although they might be in some J p 
with p > 1: an operator A is in J v iff Tr|^4| p < oo in any orthonormal basis). In order 
to compute determinants of such operators, an extension of the cumulant expansion is 
needed which, in fact, corresponds to a regularization procedure |H],[F3|]: 
E.g., let A e J v with p < n. Define 



R n (zA) = (1 + zA) exp (^2 t-jf- A") - 1 



(A.21) 



as the regulated version of the operator zA. Then the regulated operator R n (zA) is trace- 
class, i.e., R n (zA) e J\. Define now det n (l + ;zA) = det(l + R n (zA)). Then the regulated 
determinant 



JV(zA) 

det n (l + zA) = [] 



<, wa^ AfU-^-(a))' 

v l + zA i (A))exp JJi 

\k=l 



k 



< oo 



(A.22) 



81 



exists and is finite. The corresponding Plemelj-Smithies formula for det n (l + A) results 
from the standard Plemelj-Smithies formula (|A.14j) by simply setting Tr(A), Tr(A 2 ), . . ., 
T^A 71 " 1 ) to zero J7J . 



See also Ref. 78 where the Fredholm determinant 



A(A) = fi ( 1 - y) 
k=o V A */ 

is regulated — in the case fi = d/m > 1 — as a Weierstrass product 



(A.23) 



a(a) = n 



k=0 



X_ 



cxp 



_*1 

2A1 



+ ••• + 



AM ' 

Ma? 1 . 



(A.24) 



Here {A.,} are the eigenvalues of an elliptic (pseudo)-differential operator H of order m on 
a compact or bounded manifold of dimension d (with < A < A x < ■ ■ ■ and A^ f +oo) 
and \p] denotes the integer part of ji. Eq.( |A.24|) is a unique entire function of order \i with 
zeros at {A&} and subject to the normalization conditions 



lnA(O) 



d 

dX 



lnA(O) 



d^ 



lnA(O) = 



(A.25) 



Clearly Eq.( [A.24[) is the same as (|A.22|) ; one just has to identify z = —A, A = 1/H 
and n — 1 = [/i]. An example is the regularization of the spectral determinant A(E) = 
det [(E — H)] which - as it stands - would only make sense for a finite dimensional basis 
(or finite dimensional matrices). In Ref.|79[ the regulated spectral determinant for the 
example of the hyperbola billiard in two dimensions (thus d = 2, m = 2 and hence = 1) 
is given as 



A(E) = det [(E - H)Q(E, H)] 



(A.26) 



where Q(E, H) = — H 1 e sH 1 . Thus the spectral determinant in the eigenbasis of H (with 
eigenvalues E n ^ 0) reads 



A ^=n(i-|) 



e E/En < oo 



(A.27) 



Note that H 1 is for this example of Hilbert-Schmidt character. 
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B Exact quantization of the n-disk scattering problem 



In this Appendix (which is based on M. Henseler's diploma thesis |4l| where also the 
corresponding formulas for the three dimensional n-ball scattering problem can be found, 
see also |5^,i8|) we will construct the scattering matrix for the scattering of a point particle 
from n circular hard disks which are fixed in the two-dimensional plane. The basic ideas go 
back to Lloyd's multiple-scattering method ]53], an application of the KKR- method |51 



to three-dimensional band structure calculations as the limiting case of n disjunct non- 
overlapping muffin-tin potentials (see also Ref. [|9] for the translation of these methods to 
the infinite two-dimensional Sinai-billiard) and to the work of Gaspard and Rice|l3]|, who 
introduced the techniques reported below to the scattering problem of a point particle 
from three equal disks in the two-dimensional plane. Here we will present a generalization 
of these methods to the scattering from n non-overlapping disks of - in general - different 
sizes. 



B.l The stationary scattering problem 



As stated in Sec.^|, the quantum-mechanical description of the scattering from n hard disks 
will be performed in the framework of the stationary scattering theory. Let ip k (r) be a 
solution of the scattering problem (for a fixed incident wave vector k). The decomposition 
of i/j into a sum over complex exponential (angular) functions 



im(f-$fc) 



(B.l) 



($£: and $ r are the angles of k and f, respectively, in the global coordinate system) leads 
to (V^+A; 2 )^(r*)=0. The corresponding separation of a plane wave in two dimensions 
into angular eigenf unctions reads: 



e ik-r = e ifcrcos4v = £ ^m$ r e im(|-* fc ) _ (B.2) 
m=— oo 

The ordinary Bessel and Hankel functions (Ji(z) = \(H[ (z) + h\ 2 \z))) of integer order 
satisfy the expressions (for \z\ ^> 1): 



H\ 2 \z)~J?-e-'< z -¥-T;) incoming, (B.3) 
h\ 1] {z) ^ ^-e +i ^ z -¥-D out-going. (B.4) 

The to-be-constructed solution can be written as a superposition of incoming and out- 
going spherical waves (kr 3> 1) 
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1 oo 



(B.5) 



where the matrix S is the scattering matrix of the two-dimensional scattering problem. 



B.2 Calculation of the S-matrix 



In order to describe a generic configuration of n disks we use the following notation (see 
Fig. |3.1|) : The index j G {!,-■■ ,n} labels the j th disk whose radius is aj. The distance 



between the centers of the disks j and j' is called Rjji = Rj>j- To specify the n disks we 
introduce n+1 different coordinate systems. First of all, a global coordinate system (x, y) 
is chosen with its origin in the neighborhood of the n disks. In case of symmetrical systems, 
as, e.g., three equal disks at the corners of an equilateral triangle, the origin is best placed 
in the center of symmetry. In order to fully use the symmetry of such configurations n 
local coordinate systems {x^\ yW) are introduced whose origins are placed at the centers 
of the n disks, respectively. The axes of these coordinate systems are chosen in such a 
way that they fully respect the symmetry of the configuration. The spatial vector to the 
center of the disk j, as measured in the global system, is called Rj, Rj is its length and 
its angle. Vectors called Sj or S are surface vectors. The unit vector Rj), = Rj),/Rjji 
is pointing from the center of disk j to the center of disk j', as measured in the (j)-system, 
ctjij is its corresponding angle. In general, vectors with an upper index (j) are measured 
in the (j)-system, vectors without upper index are measured in the global system. 

The Green's functions satisfy the differential equation (V^ + k 2 )G(f,f') = 5 2 (f— f'). In 
two dimensions the free Green's function reads 111 31: 



G(r,f') = - l -H^\k\f-f'\) . (B.6) 
For the following, we will apply the Green's formula: 



J d 2 r (0(r-)V r 2 ^(f) - ^(f)V r 2 0(r")) = J dS- (<f>0)V s iP(S) - s <f>{§)) 

V dV 

where V is the integration volume and dV denotes its boundary. After inserting the ex- 
pansion coefficients ip^ir) from ( B.l| ) and the (free) Green's function in the last equation, 
one finds: 



dv 



dS-[r m (S)V s G(S,f')-G(S,f')Vs^ m (S)] = [\k m{fl) p|y (B.7) 



The integration volume is chosen as a big disk whose center is in the origin of the global 
coordinate system and whose radius is large enough that the asymptotic equations ( |B.3| ) 
and (|B.4|) hold for the points far away from the origin but inside the integration volume. 
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From the large disk the small n disks (as given in the concrete disk configuration) are 
excluded; however, the radii of these subtracted disks have been increased by a small 
increment e > in comparison to the original disks. In the end, the case e — > is 
considered. In order to construct the S-matrix, one has to work out ( p.7|) for two different 
cases |T3]. In the first case the point f' is on the surface of the (original) scattering disk 
j, such that it is now outside the integration volume V. In the second case r' is in the 
integration volume; however, so far away from all n disks that the asymptotic equations 
( |B.3| ) and ( |B.4j ) are then valid. The boundary of V splits into n+1 disjunct regions: Into 
the outer layer of the large disk, d^V, and into the boundaries djV of the n subtracted 
disks which contain and cover the scattering disks. 



B.2.1 First case: r' = Xj £ boundary of disk] 

Because of the Dirichlet boundary conditions, the wave function vanishes on the boundary 
of the scattering disks; however, its gradient does not vanish there: 

+oo 

^ k m(X 1 ) = ; nrVitixf)^ £ B^e^V (B.8) 



m '= — oo 



Here the unit vector fij is chosen to point perpendicularly to djV into the complementary 
region of V. Note that \Xj | = aj. Furthermore, 9j labels the direction of as measured 
in the local coordinate system of the disk j. The coefficients B m ^ ; are unknown so far. 
Eq.( [B.7| ) now reads: 



o = 4 + E 4 ■ (B-9) 

j'=l 



The occurring integrals are: 



Il= J dS ■^i(S)V s G(S,X,)-G(S,X J )Vs^ k m (S)} (B.10) 



dooV 

4 = ~ J ds f ■ G(sfjP)V s .^ k m ( S f) . (B.ll) 
d sl v 

In the following we will repeatedly apply the addition theorems for Bessel and Hankel 



functions [80 



C n (w)e ±m ' 3 = ]T Cn^J^e^ , (B.12) 



l=— oo 



where w = \Ju 2 + v 2 — 2uv cos a, w cos B = u — v cos a, \v e 1Q | < \u\, w sin (3 = v sin a and 
C n {z) G {J n (z),Y n (z),H^(z),H^(z)}. 
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Calculation of /£ 



The calculation is performed in the global coordinate system. The addition theorem ( |B.12| ) 
is used to rewrite the free Green's function (|B.6|) . In addition, because of the large value 
for R s , the Hankel function Hi(kRs) is approximated according to ( p.4j) . The resulting 
expression and the asymptotic expression (p.5| ) for ip^S) are inserted into 1^. The terms 
proportional to S m ; cancel out, such that 



IL = e ^ Rj J2 Jm-mikR^Jmika^e™' 6 * . (B.13) 



m =— oo 



Only the coordinates of the center of the disk j are still expressed in the global system, 
whereas the coordinates on the disk surface have been transferred to the local coordinate 
system of disk j. 

Calculation of /j 

Here we work relative to the local coordinate system of disk j. Using the addition theorem 
flB.12 ) for the free Green's function ( ]B.6|) and performing the angular integration under 
the boundary condition ( p.8| ) we obtain 

OO 

Ij = E BjHPikaMikaj) , (B.14) 

l=— oo 

where all quantities are expressed in the local coordinate system of disk j. 
Calculation of I 3 -, , j ^ j' 

Working relative to the local coordinate system of disk j', we have in this case: 

1 oo 

G(s J ,,X { / ) ) = ± £ H} 1 \kXP)J l {ka jl )e- n ^'-^'^' )) . (B.15) 

l=— oo 

In writing down the last equation the addition theorem for Hankel functions has been 
used again. Here 4>^ ' = ctjj> + 7^/, where cpf ' is the angle of X.f \ ajj> is the angle of the 
ray from the center of disk j' to the center of disk j and 7^ is the difference angle. All 
three angles are measured relative to the local coordinate system of disk j'. After insertion 
into I?, and the angular integration we apply once more the addition theorems for Hankel 
functions. Then Ij, reads: 

00 

4 = ~ 1 Y L £ BjVKto j OJr(K0^k^)e i ^/-^)(-l)''e i ^, (B.16) 
1,1'=— 00 

where the entries of Ii do not depend on the global coordinate system. The j' dependent 
quantities are expressed in the local coordinate system of disk j', the j-dependent ones in 
that of disk j. 



86 



The computed integrals are now inserted into the formula ( |B.9| ), written as /£, = — J2f=i If, 
which leads to 

£Cje^=££BjM#e^ (B.17) 

l l j'M 



with C m j abbreviating the terms in Eq. flB.l3D , whereas M^,/ stands for the terms in 



( B.14|) and (|B.16| ). Equation (|B.17|) holds for all points Xj on the boundary of the disk j. 



Then, the coefficients C m j and M|,/ are normalized in such a way, that in the 1-disk case 
the new M-matrix is just the unit matrix. This corresponds to a division of the l.h.s. and 
r.h.s. of (|B.17|) by the diagonal matrix {Hi 1 \kdj)Ji(kdj)rcdj /2i}. Asymptotically (i.e., for 
|/| ^> \kdj\) the modulus of its matrix elements behaves as {H^^kdj) Ji(kdj)\ ~ l/(7r|/|). 
Therefore, this division does not affect the "trace-character" of the matrices C J , 
and B- 5 (see App. y). Thus one gets the matrix equation 

C j = B r ■ M j ' j (B.18) 

with 



\i=^ J ^ (LRj) — , (B.19) 
Hm(ka,j) Traj 

Ml£=8"'8 mm , + (1-^) aj ^ j \ H^ m/ (kR jf ) eK^-^)(-l)-' , (B.20) 



dj'H^i (kdj> 

where Rj and <pR. are the magnitude and the angle of the ray from the origin of the global 
coordinate system to the center of disk j, as measured in the global coordinate system. 
The angle atfj is the angle of the ray from disk j to disk j' as measured in the local 
coordinate system of disk j, Rjf = Rj/j is the distance of the centers of disk j and f, dj, 
dji are their radii, respectively. 



B.2.2 Second case: point r' = rEV , r large 
For this case we obtain from ( |B.7| ): 

n 

^(n = ^ + E J j- (B.21) 

3=1 

In analogy to the first case, the following abbreviations have been introduced: 

C= J dS-(^(S)V s G(S,r)-G(S,f)V s ^ k m (S)), (B.22) 

dooV 

~ I d*i • Gft, r W) f„,< ft) • (B.23) 



9,-V 
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Calculation of 7^ 

7^ can be calculated in close analogy to 7^. A single application of the addition theorems 
( B.12|) yields 7^ = J m (fcr) e im * r , where $ r is the angle of f in the global coordinate 



system. 



Calculation of 7J 



We have 

1 oo 

G(sj, ) = ± HPWJt-tikfyMkaj) e ±z '^ e^*^*^ , (B.24) 

41 {,{'=-00 

where the addition theorem for cylindrical functions has been applied twice. The angle 
<p Sj of Sj is measured relative to the local coordinate system of disk j. After integration 



over this angle, we get: 



_ oo 

£ H} 1 \kr)J l MkR j )Mka j )Bj e u <*--***> , (B.25) 

Z1 {,{'=-00 

where all quantities, except a,-, are now defined with respect to the global coordinate 
system. 



Both integrals are now inserted into Eq.( |B.21] ). Taking into account ( |B.3[ ) and ( |B.4| ), one 
gets Eq.( p.5| ) for fcr> 1. The S-matrix can now be written as 



S (n) = 1 - iB J • B j , (B.26) 

where we have introduced the superscript (n) in order to indicate that the S-matrix refers 
to the n-disk scattering problem. The matrix D J in the last equation is given by 



D w = -7rajJ m /_ m (A;i^)J m (fcOj-)e im * fl ? . (B.27) 



Using QB.18 ), we finally get the (formal) expression for the S-matrix which will be justified 
in App. 0: 

S (n) = 1 - iC j • (M" 1 )^" • D J " . (B.28) 
The S-matrix S*- 1 -* of the scattering of a point particle from a single hard disk is given by 

= -il^^ > (B-29) 



H^(ka 



3 / 



as can be seen by comparison of the general asymptotic expression ( |B.5| ) for the wave- 
function with the exact solution for the 1-disk problem. 



s8 



B. 3 The determination of the product D ■ C 



In order to rewrite the determinant of the S-matrix (see Sec.^) we have to determine the 
product DC (see Eqs. ( |Bl9|) and (|R27D ). We apply once again the addition theorem 
for Bessel functions ( |B.12|) using Rfj = Rj>—Rj, where Ry and Rj are the magnitudes 
of these vectors and , and the corresponding angles, as measured in the global 
coordinate system. We find the following expressions 



l=— oo 

oo 



E Dj« = -2i 



lO 



Hp\ka 3 
Su> ■ 



l=— oo 



Using the expression (P-20| ) for M^^, we finally get for X = M — iD ■ C: 



tf/, 2) (£;a 



^5 jj 'S lv - (1-5"')- 



OjJlikOj) (2) 



H\}\kaji) ' ajiH\}\kay 



H£>,(kR jf )(-iy e 



The r.h.s. of this equation can be reformulated in terms of the scattering matrix of the 
single-disk problem, S^{kaj>), and by the complex conjugate of M, namely Xfjy (k) = 
S\j)(ka,j') (iVF^ ^(k*)) . The matrix X can therefore be expressed as the product, 

x ^ = EE( M f)* Y l' > 



where the second factor is given by Y~J = 5 J - J '5 / - / ,S_/, _ i ,. Thus we get the formal expression 
for the determinant of X: 



Det L X(fc) = [ J] detiS (1) (fceij) ) Det L M(fc*^ 



;b.3o) 



The last step in the formal evaluation of the determinant of (as function of the wave 
number k) is the insertion of Eq. ( [B.30| ) into Eq.( |3.6| ) giving the final result (see (|3.9)): 



det , S W W = gg£W)| nd e t , S W(fa 



Det L M(A;) 



(B.31) 



coherent 



incoherent 
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C Existence of the S— matrix and its determinant in n-disk systems 



This appendix is based on M. Henseler's diploma thesis [it] and on Ref. [Ic . 



The derivations of the expression for S-matrix ( P-28|) in App. [B] and of its determinant 
(see Sec.fD are of purely formal character as all the matrices involved are of infinite size. 
Here, we will show that the operations are all well-defined. For this purpose, the trace- 
class (J7i) and Hilbert-Schmidt(j72) operators will play a central role. The definitions and 
most important properties of these operator-classes can be found in App. [A]. As shown in 
App. || the S^-matrix can be written in the following form (see (|B.26| )): 

S(») = l-iT , T = B j B j . (C.l) 

The T-matrix is trace-class on the positive real /c-axis (k > 0), since it is the product of 
the matrices D J and B J which will turn out to be trace-class or, in turn, are bounded there 
(see App. 0.1 for the definitions). Again formally, we derived in App. [B|that C J = B J 1VP J 
implies the relation B J = C^M -1 )^ . Thus, the existence of M _1 (/c) has to be shown, 
as well - except at isolated poles in the lower complex /c-plane below the real /c-axis and 
on the branch cut on the negative real /c-axis which results from the branch cut of the 
defining Hankel functions. As we will prove later, M(/c) — 1 is trace-class, except of course 
at the above mentioned points in the /c-plane. Therefore, using property (e) of App. |A|.2 
we only have to show that DetM(/c) ^ in order to guarantee the existence of M _1 (/c). 
At the same time, M _1 (/c) will be proven to be bounded as all its eigenvalues and the 
product of its eigenvalues are then finite. The existence of these eigenvalues follows from 
the trace-class property of M(/c) which, together with DetM(A;) ^ , guarantees the 
finiteness of the eigenvalues and their product. 

We have normalized M in such a way that for the scattering from a single disk we simply 
have B = C. Thus the structure of the matrix C J does not dependent on whether the 
point particle scatters only from a single disk or from n disks. Hence the properties of 
this matrix can be determined from the single disk scattering alone. The functional form 
( P.19|) shows that C cannot have poles on the real positive /c-axis (k > 0) in agreement 



with the structure of the S^-matrix discussed in App.fB| If the origin of the coordinate 
system is placed at the origin of the disk, the matrix is diagonal. In the same basis C 
becomes diagonal. Thus one can easily see that C has no zero eigenvalue on the positive 
real /c-axis and that it will be trace-class there. Thus neither C nor the 1-disk (or for 
that purpose the n-disk ) S-matrix can possess poles or zeros on the real positive k- 
axis. The statement about follows simply from the unitarity of the S-matrix which 
can be checked easily. Since, for real positive k, S^/c) = S _1 (/c), we have |S(/c)| = 1 
on the real axis, such that poles (and also zeros) of S are excluded there. Actually, for 
the exclusion of poles and zeros on the real positive /c-axis, only the weaker condition 
that |detS (n) (fc)| = 1, k > 0, is needed. That this is fulfilled for all non-overlapping n 
disk systems is obvious from the final expression (|3.9| ) for detS^ in Sec. |j. This formula 
even holds for DetM(/c) — > if k approaches the real positive axis, since then DetM(fc*)' 
approaches zero as well, such that both terms cancel in formula (|3.9|) . Thus the fact that 
|detS^(/c)| = 1 on the positive real /c-axis cannot be used to disprove that DetM(/c) 
could be zero there. However, if DetM(/c) were zero there, the "would-be" pole must 
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cancel out of S^ n ^(A;). Looking at formula flB.28|) , this pole has to cancel a zero from 
C or D where both matrices are already fixed on the 1-disk level. Now, property (g) 
and (f) of App. [A].2 leave for M(fc) (provided that M — 1 has been proven trace-class) 
only one chance to cause trouble on the positive real fc-axis, namely, if at least one of its 
eigenvalues (whose existence is guaranteed) becomes zero. On the other hand M still has to 
satisfy C J = B J M J J with C completely determined by the 1-disk scattering alone, where 
DetM = 1 everywhere in the fc-plane. The fact that C^(k) cannot have zero eigenvalue 
for k > can be used to show that the following inequality holds for the modulus of the 
diagonal matrix element, |C m £,(A;)| > 0, for the state \m) of any orthonormal basis. Now 
choose as the basis the eigenbasis of 'M J mJ ^,(k) and \m) as the state there M(fc) has a 
candidate for a zero eigenvalue in the \m,j) space. Comparing the left and the right-hand 
side |C m £j(A;)| = |B m / M^| one finds a contradiction if the corresponding eigenvalue of 
M(fc) were zero, i.e., the l.h.s. would be greater than zero for k > 0, whereas the r.h.s. 
would be zero. Hence, such a zero eigenvalue cannot exist for k > 0, hence DetM(/c) ^ 
for k > 0, hence M(&) is invertible on the real positive A;-axis, provided M(/c) — 1 is 
trace-class. From the existence of the inverse relation B J = , '(M -1 ) J ' J ' , the trace-class 
property of G 7 to be shown and the boundedness of (M" 1 )-" follows the boundedness of 
B J and therefore the trace-class property of the n-disk T-matrix, T^(k), except at the 
above excluded fc-values. 

What is left to prove is 

(a) M(&) -16 Ji for all k, except at the poles of H^(kaj) and for k < 0, 

(b) C^(k),Ty(k) G J\ with the exception of the /c-values mentioned in (a), 

(c) T^(kaj) G J\ (with the same exceptions as in (a) and (b)) 

(d) M _1 (A;) does not only exist, but is bounded. 

Under these conditions all the manipulations of Sec.f|, are justified and S^ n \ as in (|3.2|), 
and detS 1 -™), as in (13. 91), are shown to exist. 



C.l Proof that (ka,j) is trace-class 



The S-Matrix for the j th disk is given by ( ggg ). Thus V= - iT W {kaf) = S^{kaj)-1 is 
diagonal: 



-2Ji(ka,j) 
H ( ^\ka j ) 



Vu-fc^P 2 . (C2) 



Hence, we can write V = U|V| where U is diagonal and unitary, and therefore bounded. 
What is left to show (see property (a) of 0.1) is that |V| G J\. This is very simple since 
we can now use the second part of property (d) of App. [SJ.1: we just have to show in a 
special orthonormal basis (here the eigenbasis) that 
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+oo 



E iv|« = E 2 

l=— oo l=— oo 



Ji(ka,j 



Hp-\k aj ) 

as |V| > by definition. The ordinary Bessel and Hankel functions of integer order satisfy 



< oo 



(C.3) 



J-n(z) 



(2), 



(C.4) 



v — > oo 



Thus: 



1 fez 
V2^ V2z/ 



v real : 



7TZ/ \ZV 




ez 



(C.5) 



tr(|V|) S 4gi(5!|i! 



2/ 



2E(^ 



,2/ 



(C.6) 



These a; satisfy: a^ < a; < 1 for I > l and l > _ The series X^o( a 'o) 2/ converges, 
and hence also the sum X^o( a z) 2 ' as ^ * s bounded from above by the previous sum. That 
means that |V| G Ji and (because of property (a) of App. 0.1) — 1 G J{, as well. 
This, in turn, means that detS^ (kcij) exists (see property (i) of App. [A[2) and also that 
the product U]=i detS (1) (fcaj) < oo in the case where n is finite (see property (d) of the 
same Appendix). The limit lim^oo does not exist, in general, as the individual terms 
detS^^/ca.,) can become large, of course. 



C.2 Proof that A(k) = M(k) - 1 is trace-class 



The determinant of the characteristic matrix lSA(k) is defined, if A(k) G J\. In order to 
show this, we split A into the product of two operators which - as will be shown - are both 
Hilbert-Schmidt. Then according to property (b) of App. |A].l the product is trace-class. 

Let therefore A = E • F, where A follows from ( |B.20| ). In order to simplify the decompo- 
sition of A, we choose one of the factors, namely, F, as a diagonal matrix: 

A#' = E#'if, 1$; = Fist? 8 W (C.7) 

and 



H^ikaj 



Already this form leads to the exclusion of the zeros of the the Hankel functions H^lkaj) 
and also the negative real fc-axis (the branch cut of the Hankel functions for k < 0) from 
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our final proof of A(k) G J\. First, we have to show that ||F|| 2 = X^/F+F);/ < oo. We 
start with 



i=i «=o |#, v (HOI 2 i=i «=o 



This form restricts the proof to n-disk configurations with n finite. Using the asymptotic 
expressions ( |C.5| ) for the Bessel and Hankel functions of large orders, it is easy to prove 
the absolute convergence of J2i a z i n the case a > 2. Therefore ||F|| 2 < oo and because of 
property (d) of App. |A].l we get F 6 J 2 . 

Using the decomposition ( |C.7| )) and the definition of F ( |C8| ), the second factor E, is 
constructed. We then have to show the absolute convergence of the expression 



],]'=1 \ J / 1,1 — — OO \ 1J -' 

m 1 



m/(\kaa f \)\ 

in order to prove that also E 6 J 2 . This is fulfilled, if J2i v a «', < oo, where 



\H£){\kaa r \)\ 



Necessary conditions for the convergence of the double sum over aw are: J2i> < °o 
as well as J2i a w < oo. For the case I — > oo , V fixed, we obtain with the help of the 
asymptotic formulas (|C.5|) the expression: 

\ -2/' 
e\kR jjt \ X -' 

7T'" 



/ — ► oo : a,^l 7/ ) Ul^Yil + ir-^^-Y . (C.12) 



For any e > this yields the estimate: 



MO < (2/f fil±^) , />/ with7<e. (C.13) 

For x = (1 + e)a,j/Rjj> < 1, the series S^o^ 2 ' converges absolutely. As X^(2Z) 2/ 'x 2 ' — 
(xd/dx) 21 J2i% 21 < oo, the series MO converges absolutely, as well. Therefore we 
have the absolute convergence of J2i &w for % < i?^/ with fixed /' in the limit e — > 0. In 
the opposite case, I' — > oo , Z fixed, the absolute convergence of J]/' a zz' f° r f %' < -Rjj' 
can be proven analogously. 

We must of course show the convergence of a w f° r the case 1,1' — > oo. Using again the 
asymptotic behavior of the Bessel and Hankel functions of large order we get the following 
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proportionality for /, V — > oo: 

The double sum Y^i°i'=o a «' * s convergent, if J2i°i'=o hi' converges. In order to show this, we 
introduce two new summation indices (M, m) as I + 1' = 2M and I — V = m. Hence, we 
have 



oo oo 2M 

E h lV = E E C Mrn (C.15) 
l,l'=0 M=0 m=-2M 

with 



(2MT I a, V IM+f V a a V' M - ¥) , C16 > 

Mm (M+f) 2 l" + f)(M-f) 2 (*'-?> \2%./ ' 1 ' 

For sufficiently large M, the powers occurring in the last expression can be approximately 
estimated with the help of the Stirling formula, n n ~ n\ e n / In this way, we get for 
M -> oo: 



. \ M- 



m \ 2 
2 



Hence, the total sum reads 



-2M oo / / . , a n \ 2M ~~ 



2 



E E c Mm <2nj: ((^5^ ] , (Clf 



M=0m=-2M M=0 \ \ u 3f 



where the sum over m has been performed with the help of the binomial formula. The 
remaining series in ( |C.18|) converges for aj + ^ay < Rjf. Therefore, under the stated 



conditions J2ii> a w converges absolutely, as well. We finally get the desired result: 
The operator E belongs to the class of Hilbert-Schmidt operators (j^), if the conditions 
-a,-' + cij < Rjj>, (1 + e)aj/ < (1 + e)^aj> < Rjf and (1 + e)aj < Rjj> are met in the 
limit e — ► 0. In summary, this means: E(k) ■ F(k) = A(k) G J\ for such finite n disk 
configurations for which the disks neither overlap nor touch and for those values of k 
which lie neither on the zeros of the Hankel functions H^\kaj) nor on the negative 
real /c-axis (k < 0). The zeros of the Hankel functions H^(k*aj) are then automatically 
excluded, too. The zeros of the Hankel functions H^{kaaj) in the definition of E are 
canceled by the corresponding zeros of the same Hankel functions in the definition of F 
and can therefore be removed. A slight change in a readjusts the positions of the zeros 
in the complex fc-plane such that they can always be moved to non-dangerous places. For 
these ("true") scattering systems the determinants DetM(fc) and DetM(fc*)''' are defined 
and can be calculated with the help of one of the cumulant formulas given in App. |A].2, 
e.g., by the Plemelj-Smithies formula ( |A.14j ) (with Det = e Trlog , see ( |A.13| ), for small 
arguments) or by Eqs.( |A.i"0|) or ( |A.12j ) if M or A can be diagonalized. 
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C.3 Proof that C J and D J are trace- class 



The expressions for D J and C J can be found in ( [B.27|) and ([B.19| ). Both matrices contain 
- for a fixed value of j - only the information of the single-disk scattering. As in the proof 
of T^ 1 ) G Ji, we will go to the eigenbasis of In that basis both matrices D J and C- 7 
become diagonal: 



mm' 



~7ldj J m (^kdj 



-im<J>f 



p 3 _ p i™*fl, 



2i 



Jim 



(C.19) 
(C.20) 



Using the same techniques as in the proof of T« G Ji, we can show that C J and D J 
are trace-class. In summary, we have D J G J7i for all k since the Bessel functions which 
define that matrix possess neither poles nor branch cuts. The matrix C J G J\ for almost 
every k, except at the zeros of the Hankel functions H^(kaj) and the branch cut of these 
Hankel functions on the negative real fc-axis (k < 0). Note that the values of trD J or trC- 2 , 
are finite and the same whether one uses the non-diagonal expressions ( |B.27| ) / ([B.19Q or 
the diagonal ones (|C.19|) /( p.20|) . This is, of course, in agreement with property (e) of 
App.0.1. 



C.4 Existence and boundedness of M l {k) 



As M(fe) — 1 G J\ except at the zeros of H^(kaj) and on the negative real fc-axis (k < 0), 
M _1 (fc) exists everywhere, except at the points mentioned above and except at A;- values 
where DetM(/c) = 0. In other words, except at the poles of the S^ n \k) matrix, see Eq.( |3.9|) . 
With the exception of the negative real axis and the isolated zeros of H^(kaj), M(A;) is 
analytic. Hence, the points of the complex /c-plane with DetM(/c) = are isolated. Hence, 
DetM(A;) ^ almost everywhere. Thus, almost everywhere, lSA(k) can be diagonalized 
and the product of the eigenvalues weighted with their degeneracies is finite, see App. |A|.2 
for both properties. Thus M~ 1 (fc) exists and can be diagonalized as well. Hence, all 
the eigenvalues of M^ 1 (/c) (and their product) are finite in the complex fc-plane, where 
DetM(/c) is defined and nonzero. Thus M _1 (/c) is bounded (and DetM _1 (/c) exists) almost 
everywhere in the complex fc-plane. 

In summary, the formal steps in the calculation of the n-disk S-matrix (see App. 0) and 
its determinant (see Sec.^) are all allowed and well-defined, if the disk configurations are 
such that the disks neither touch nor overlap. 
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D Comparison to Lloyd's T-matrix 



As mentioned in Sec.[3], Lloyd has constructed a formal expression for the T-matrix of a 
finite cluster of muffin-tin potentials in three dimensions, see Eq.(98) of Ref. f53|. Tran- 



scribed to the case of a cluster of n disk-scatterers fixed in the two-dimensional plane, 
Lloyd's T-matrix reads as T(fc) = C(fc) (m(/c)) _1 D(fc) with 



CJ = Jm-iikRj) ( ] , (D.l) 



( —2iJi(ka j y 

Bi: m , = J m ^(kR f )e- im '^' , (D.2) 
M# = S*5 W + (1 - 8*) ^ Hi%(kR fj )r jf (l, I') , (D.3) 



where the tilde is discriminating the matrices in the Lloyd representation from the cor- 
responding matrices in the Gaspard-Rice representation, defined in Q3.3| ), ( ft.4 | ) and fl3.5| ). 
The Lloyd representation allows for a very simple interpretation. The matrix C- 7 describes 
the regular propagation (in terms of the homogenous part of the free propagator) from 
the origin to the point Rj and a one-disk scattering from a disk centered at this point, 
as given by the one-disk T^-matrix. The matrix D J describes the (regular) propagation 
back from the disk f to the origin. The matrix (M JJ ) _1 parametrizes the multiscattering 
chain. If it is expanded around 5^ , it describes the sum of no propagation and no scat- 
tering plus the propagation from disk j to disk f (in terms of the full propagator) and a 
scattering from disk f and so on. The disadvantage of the Lloyd representation is that 
the trace-class character of A = M — 1 is lost, as the terms J m (kaj) and (H^^kaj))" 1 
"stabilize" only the asymptotic behavior of the index I, but not of the index I' any longer, 
as the asymmetric Gaspard-Rice form did. The infinite determinant detM is therefore 
no longer absolutely convergent, but only conditionally. Any manipulation in the Lloyd 
representation of the matrix M and the corresponding S-matrix has therefore to be taken 
with great care. Note, however, that the (formal) cumulant expansions of M and M are 
the same as the corresponding traces satisfy Tr(A n ) = Tr(A n ). In other words, if the 
cumulants of M are summed up according to the Plemelj-Smithies form of M, the result 
of detM(fc) and detM(fc) is the same. In fact, one can derive the Lloyd representations 
C j , M jj ' and D J " from the expressions for C j , M j ' j and B j of App. @.l (see Eq. flETTp ) 



by the following formal manipulations: First, C m j and M|,/ are divided by the diagonal 
matrix {H^ (ka,j)/(— 21)}. This produces already G 7 . Second, B m /^ in the (now changed) 
relation ( |B.17|) and in flB.261 ) is rescaled as 



B <j' = Bj ' ^n T (kn \ ' (D ' 4) 

such that D J and M" emerge. Both manipulations are only of formal nature as they 
change the "trace-character" of the corresponding matrices. 
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E 1-disk determinant in the semiclassical approximation 



In App. |B.2| we have constructed the scattering matrix for the 1-disk system (see Eq.( |3.8| )): 



(sW(fa 



J-L A 

r (l) . . . u mm' 



(E.l) 



Instead of calculating the semiclassical approximation to its determinant, we instead do 
so for 



1 d 



d{k) = — — lndetS (1) (fc( 



2ni dk 



(E2) 



the so-called time delay. Recall that the corresponding T' 1 ' -matrix is trace-class. Thus, 
according to properties (j) and (c) of App. [A].2 the following operations are justified: 



1 d 



d(k) = -£-tr f lndetS (1) (A;t 

27ri dk v 



-tr 



(HW(ka 3 ) d Hi\ka 5 ) 



2vri \H%\kaj) dk H&\k aj ) 



2vri I Hm (kdj) H^{kaj 



(E.3) 



Here the prime denotes the derivative with respect to the argument of the Hankel func- 
tions. Let us introduce the abbreviation 



Xv 



Hi 2 \kaj) Hl l) {ka 



(E.4) 



Following Ref.[]29[], we apply the Watson contour method [B3J] to ( |E.3| ) (see also Sec.|5] and 
App© 



a ■ + °° n 1 r e~ iu7V 



Xv 



(E.5) 



Here the contour C encircles in a counter-clock-wise manner a small semi-infinite strip 



D which completely covers the real z/-axis Ea, but which only has a small finite extend 



into the positive and negative imaginary v direction. As in Ref.Kj, the contour C will 
be split up in the path above and below the real z/-axis such that 



a ( 1 +0 °r +le 1 + T" e"^ 

d{k) = -^{-- / dz/-— + - / dz/-— -x„} ■ (E.6) 
2tti 2i J sm(z/7r) 2i J sm{h>ir) 

\ — oo+ie — oo— ie 



# 4 In App.[F| symmetrized expressions have been Watson transformed. Thus, the corresponding 
D + only has to cover the the real positive z^-axis . 
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Then, we perform the substitution v 



— v in the second integral so as to get 



d ( fc ) = TT^r < — — / dv— / dv— -x- 

Z7ri 2i J sm(z/7rj 2i J sm(z/7r) 

v — oo+ie — oo+ie 



+oo+ie 



,2iv7T 



a M 2 / d "rb^ /'"^> 

— oo+ie 



where we used the fact that X-v = X^- The contour in the last integral could be deformed 
to pass over the real z/-axis since its integrand has no Watson denominator any longer. 
We will now approximate the last expression semiclassically, i.e., under the assumption 
kcij ^> 1. As the two contributions in the last line of QE.7Q differ by the presence of the 
Watson denominator, they will have to be handled semiclassically in different ways: the 
first will be closed in the upper complex plane and evaluated at the poles of Xu, the second 
integral will be performed under the Debye approximation for Hankel functions. We will 
now work out the first term. The poles of Xu i n the upper complex plane are given by the 
zeros of Hjp(kaj) which will be denoted by u^kaj) and by the zeros of H^ 2 '(kaj) which 
we will denote by —p£(kaj), i = 1,2,3, In the Airy approximation to the Hankel 
functions, see H30| , they are given by Eqs. (|4.2|) and (|4.3|) : 



v t {ka) = ka + e +iT/3 (ka/Q) 1/3 q e + ■ ■ ■ = ka + ia e (k) + ■■■ , (E.8) 
-P e (ka) = -ka - e- i7T/3 (ka/6) 1/3 qe + --- = -ka + i{a e (k*a))* + ■■■ 

= -(v t (k*a))* , (E.9) 

where ae(kdj) is defined in and qi labels the zeros of the Airy integral ( |F.27| ), for 
details see p9| , [44l . In order to keep the notation simple, we will abbreviate U£ = Uf (kaj) 
and T>£ = Pe(ka,j). Thus the first term of ( |E.7| ) becomes finally 



fi i +oo+ie g2ii/7r 1 00/ g2ii^7r p — 2iz^7r 



^{2 
2vri 



/ =2 % E (3-^ + 3-^1. (E.10) 



— oo+ie 



In the second term of ( |E.7|) we will insert the Debye approximations for the Hankel 
functions [1501: 



H^ 2 \x)~\l — exp ( ±i\/ x 2 — v 2 =F arccos — =F i— ) for \x\ > v , (E.ll) 

V vrVa; 2 -v 2 \ x A. 



Hj?- /2) (x)~Ti\ — exp(-Vi/ 2 -x 2 + t/ArcCosh-) for Id < v . (E.12) 

V Tiyv 2 - x 2 V x) 

Note that for v > kaj the contributions in Xv cancel. Thus the second integral of ( |E.7| ) 
becomes 
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a + °° +fcQj (—2i) d ( v \ 

S/ dz/ ^ = Si / du -^ik u2 ~ u 8x0008 



-oo — feaj 



a? 



- — / du^k 2 a 2 - v 2 + • • • = -- jf fc + ■ ■ • , (E.13) 



/u7r 7 v j 2 



-fca-i 



where • ■ ■ takes care of the polynomial corrections in the Debye approximation and the 
boundary correction terms in the v integration. 

In summary, our semiclassical approximation to d(k) reads 



d(fc) = 2a 7 -V — - + — - —k H . (E.14) 

Using the definition of the time delay (|E.2p , we get the following expression for detS^ (kaj): 

lndetS (1) (A;a,) - lim lndetS (1) (fc a,-) 

-i2vr + 2(i2vra,) V + + ■ • • 

2 i^x 1 1 - e i2?r ^( fca i) 1 - e - i27rp HK) / / 

oo * _ J 

-2vriiV(A;) + 2^ / dA;^ {-ln(l-e i27r ^( fca ^) +ln(l-e- i2 ^ (fca ^)} + ■ ••, (E.15) 



o 



o 



where in the last expression it has been used that semiclassically 4rU(.{kaj) ~ 4rT>i(kaj) ~ 
a,j and that the Weyl term for a single disk of radius aj goes like N(k) = na 2 k 2 / (4tt) + ■ ■ ■ 
(the next terms come from the boundary terms in the ^-integration in (E.13|) ). Note that 
for the lower limit, k Q — > 0, we have two simplifications: First, 

lim 8^,(^0%) = Km H ™ ^^- 5 mm , = 1 x 5 mm , Vm, m 
~» lim detS (1) (A;oa,) = l . 

fco-»0 

Secondly, for fco ~~ ^ 0, the two terms in the curly bracket of ( |E.15|) cancel. Hence, we finally 
obtain the semiclassical result for the determinant of S^(kaj) 



2 

poo I -| p — 2'mVi (ka,j) N 



rr°° ( i _ e - 

det S«(fca J ) ^ g^AWiifziV 7 (E16) 

n£i (l - e 2i7r ^( fea ^)J 

which should be compared with expression ( |4.1|) of Sec.^. For more details we refer to 
App.0 
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F Semiclassical approximation of two convoluted A-matrices 



In this appendix we introduce the necessary apparatus for the semiclassical reduction of 
Tr[ A m (k) ] for the n-disk system where 



A*f = (1 - 8 ir ) A -lY ^~ Va ^ H&ikRjf) ■ (F.l) 

a r H K v '{ka f/ ) 



As usual, dj, dji are the radii of disk j and f, 1 < j, f < n, Rjji is the distance between the 
centers of these disks, and aj>j is the angle of the ray from the origin of disk j to the one 
of disk j', as measured in the local coordinate system of disk j. The angular-momentum 
quantum numbers I and I' can be interpreted geometrically in terms of the positive- or 
negative- valued distances (impact parameters) l/k and I'/k from the center of disk j and 
disk j', respectively, see Figs. |K.lHF.3| . 

The semiclassical approximation of the convolution of two kernels J2v ^\<l« contains 
all (but one) essential steps necessary for the semiclassical reduction of the quantum 
cycles and traces themselves. What is missing is the mutual interaction between successive 
saddles of the quantum itinerary, including the final saddle which "closes" the semiclassical 
open itinerary to a period orbit. This is studied in Sec.[5| 

The idea here is to construct the convolution of the two kernels J2v Ajjj Aj,/» and then to 
compare it - in the case j ^ j" - with the single kernel Aff,, (see ( [F. 1|) ) in the semiclassical 
limit, where the Hankel function H^ v ,(kRjjir) is evaluated in the Debye approximation 
(|K.11|) to leading order f8(J. Let us start with 



V A"' A j ' j " 
7 j w A m" 



ajJi(ka,j 



dfiHi}\ka 



l" p i(iay j—l" ay ,■/;) 



■lfe 



x V (-!)'' Jv{ka ^ trW ao...Mr(i) 



Z'=oo 
+oo 

=w$" £ ( 

l'=oo 



-H}_\, (kRjjf)Hfli„ (kR 



. Jfj „) e u ' (0, J"i'- a «') 



\V Jl'( ka j') „(1) /, D srr(l) 



H&ikRjfWlAkRfj") e u 'K'v-,y) (F.2) 



H^kR^H^kR^e 



(i) 



■jl' (aijii ji — ctjji) 



l'=0 



+ H^^kRjj^H^i^^kRfj") e 



(i) 



-«'(a 



(F.3) 



where we have introduced the abbreviations Wj'J, for the /'-independent pieces and the 
weight factor d(V) = 1 for I' ^ and d(0) = 1/2. We have symmetrized this expression 
with respect to I' for simplicity using that J_ii(kaj/) = (—l) l 'Ji'(kaj>) and H^,(kaj') = 
(— Vf' Hp (kaji)^ valid for V integer. We will furthermore abbreviate Aatjt = aj<<j> — ajjt 
where < Aay < 2n. However, in order to be able to get three domains for this angle 
(which we will later identify with the three different cases: specular reflection from disk 
j' to the right (see Fig. |b'.l| ), to the left (see Fig. |b'.2| ) and the ghost "tunneling" case (see 
Fig. F.3|) ) we define Ao^ = aifiji — a^i — an where a = 0, 2, 1, respectively, and balance 
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this by multiplying accordingly the right hand side of Eq. (|F.3| ) with the phase factor 
(— \) la which is only nontrivial for a — 1. We denote this nontrivial phase by (— l)~ i<J 
where a' = a for a = 1 and zero otherwise. The three choices for the value of a are still 
equivalent at this stage. 




Fig. F.l. The geometry belonging to a trajectory, j — > f — > j", specularly reflected to the right. 
Shown are the geometrical path (full line) and the shortest allowed right handed (dashed line) 
and left handed (full line) creeping paths. All paths originate from an "impact parameter" circle 
of radius \l/k\ centered at disk j, then contact the surface of disk f (of radius a') and end on 
an "impact parameter" circle of radius \l" jk\ centered at disk j" . Note that the impact radii 
do not have to be equal to the disk radii, a and a" . 




F.l The Watson contour resummation 



It will be shown that ([F.3|) contains, in the semiclassical limit kay 3> 1 - depending on 
the pertinent angles and distances - a classical path or a possible ghost path and two 
creeping paths, all starting under the impact parameter l/k with respect to the origin 
of disk j and ending at an impact parameter I" /k relative to the origin of disk j". This 
calculation will be performed with the help of the Watson contour method [|64] ,f29[ 



i.e., 
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Fig. F.3. The geometry belonging to a ghost trajectory j — > f —* j" which passes straight through 
the disk j' (of radius a'). Shown are the geometrical ghost path (full line and short-dashed line), 
and the shortest allowed right handed (dashed line) and left handed (full line) creeping paths. 
All paths originate from an "impact parameter" circle of radius \l/k\ at disk j and end on an 
"impact circle" of radius \l" /k\ centered at disk j". Note that the impact radii do not have to 
be equal to the disk radii, a and a" . 

under the replacement 

£ {-l) l '^d{l')X v = ± I di/ — L — e--'-' X v , . (F.4) 
,T^n 2i J sin(z/'7r) 



Here 



Jy (kdji) 



H\%{kRjj>)Hfl v ,{kR jl: jn) e 
+ HillikRj^Hty^ikRff,) e 



(kdji) 



-Yv 



(F.5) 



stands for the integrand and Yy for the symmetrized square bracket in Eq. ( |F.3| ). The 
contour C + is the boundary of a narrow semi-infinite strip D + which completely covers 
the positive real z/'-axis. 



® 



Fig. F.4. The Contour C+ in the complex z^-plane. 

C + has been chosen in such a way that it encircles in a positive sense all poles of the 
Watson denominator sin(z/7r) at v' = 1, 2, 3, • • • exactly once (see Fig. |F.4j ). At v' = the 
weight factor d(0) = 1/2 is taken into account by a principle value description, i.e., by 
the average of the two contour integrals whose paths cross the real i/-axis symmetrically 
just to the right and left of the point v' = 0. A precondition on the validity of the Watson 
replacement is the analyticity of X u i in this strip D + . This is the case if D + has been 
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chosen narrow enough in the imaginary v direction that the poles of X u >, the zeros of 
the Hankel function H^\kdj>) and H^l^kdj/) lie either above or below the strip in the 
complex z/'-plane [for k real and positive]. The contour can now be split into four parts: 



0+ie 



/p LIS ttU r. p LIS li U r. p lis /I (J 

dv' . . - A = + / dis' , — r Xs + P / dis' . — r l, 
sin(i/7r) j sin(i/7r) J sm(z/7r) 

C+ +oo+ie 0+ie 



+oo+ie 
0-ie 



0+ie 
oo — ie 



dv 1 . . - r X„+ / du' X v , , 

V ' 0-ie V ' 



(F.6) 



where P / • • • denotes the principal value integration. The next step in the evaluation is a 
shift of the contour paths below the real z/'-axis to paths above this axis by the substitution 
is' — > —is' in the corresponding integrals: 



e -ii/W °°f ie g-ivW °+ ie ii/vrcr' 

d*/--— X„, = - / dis'——X u ,- / dz/'--— 
sm(z/'7r) J sm(z/7r) J sm(z/7r) 

C+ V ' 0+ie V 7 -oo+ie V ' 



+ P / dl/ 



0+ie 



sin(z/7r) 



(F.7) 



We insert JT„/ = (J v '(kaf)/ H$ (kaji))Y v i and use that for general complex-valued angular 
momenta z/, the transformation laws for the Hankel and Bessel functions read 



H {2) u ,(ka f ) = e- il/7r H$\ka.jf) , 
such that 



(F.8) 
(F.9) 



= - i e-"* sin^V) ^l^ 03 ' , (F - 10) 

H_}l(ka r ) Hj}\ka r ) H^\ka r ) 

H W u ,(k ajl ) H^{ka r ) ' 

Furthermore, by definition, we have Y_ u = Y +v . Thus 



oo+ie 



- i dv' 
2i 7 

c+ 



sin Z/'7T 



+ 



2i y 

— oo+ie 

+oo-ie 

jP /d,/ 

— oo+ie 



J u i(kcij> 



sm(i/n)Hj}\kaj,) 



-ii/'ir(l-cr') 



^(W) 



J v i{kcij 



Hj}\kaj 



(F.12) 
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where ( |K.11| ) and the symmetry of Y v > has been used in order to reflect the resulting 
sin(z/7r)-independent integrals at v' — such that they combine to the symmetric integral: 



0+ie 



0-ie 



+00 — ie 



+P J dv' - ■ ■ = J dv' ■ ■ ■ + P J dv' - ■ ■ + P J dv' ■ ■ ■ + J dv' ■ ■ 



(F.13) 



-0+ie 



0-ie 



Furthermore, in the case a' = 1, the identity e ll/7r = e _u/7r +2i sin(z/7r) has been employed 
in order to group the terms resulting from the paths below the real i/-axis into the terms 
belonging to the paths above this axis. Altogether we have so far that 



-oo+ie 



£Ag'A$' = W#M i "/'d^' 6 "H\ J m kaf) 
^ ' ] 2i J 8m{j/ir) H $\ka r 



x 



oo+ie 

(!) 



+ H^, (kRjj' ) H^l,_ v , {kRj'j" ) e-^'^i'- 



du > e -i^(i- CT ') 



— oo+ie 
X 



H$\kaj. 



H^\ka r 



J v i(k<lj 



+ H^kR jy )HV_ v ,{kR nil )e 



(F.14) 



Note that both integrals on the right hand side exist separately. The one with the Watson 
"sin" -denominator is finite, because the zeros of the sin(V7r) function in the denomina- 
tor are avoided by the +ie prescription and because the rapid convergence of the ratio 
J u i(kdji) I 'H^' \kdji) counterbalances the diverging Rjj> and -Rj/j/z-dependent Hankel func- 
tions, as long as the disks do not touch. This is basically the same argument by which 
one can show the existence of the sum on the left hand side. However, we do not have 
to prove this separately, because we already know from App.|C| that A is trace-class. The 
existence of the principal value integral follows from the symmetric nature of the path 
and of the integrand (see below for more details). 

It will be shown that the term with the Watson "sin" -denominator, — l/[2i sm(u' 71)] = 
J2^=o e 1 ( 2n+1 ^ u n , will lead in the semiclassical reduction to paths with left handed and right 
handed creeping sections around the middle disk f [where the index n counts further 
complete turns around this disk]. On the other hand the term without this denominator 
will give either a semiclassical path specularly reflected from the disk f (to the left or 
right) or a ghost path passing undisturbed through disk f. 



F.2 The integration paths 



Thus the third step is to close the path of the "sin" -dependent integral in the upper 
complex //-plane. 
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Fig. F.5. The path for the "sin" -dependent integral. The lines denoting the zeros of Hi 1] (ka) m 

(2) 

the upper and of Hy (ka) in the lower complex z^-plane are shown as well. 

For given values of oijj/, aj»j>, l/k and l"/k, i.e., for a given geometry, this selects which 
value of a has to be inserted into Eq. ([F.14|) . The reason is that the closing of contour 
will be performed under the condition that the corresponding semicircular integration 
arc vanishes, such that the integral is solely given by its residua which are here the 
zeros v\ (P. = 1, 2, 3, . . .) of the Hankel function H^\kcij>) in the upper complex z/'-plane. 



At "optical boundaries" this clear separation is not possible |3T|j3^|. This is the realm 
of "penumbra" scattering. In order not to be plagued by these difficulties, we exclude 
geometries which allow for grazing classical paths from the further consideration. 

In the Airy approximation to leading order, the zeros of these Hankel functions are given by 
Eq.( f4.2|) , modulo O ([^ctj']" 1 ^ 3 ) corrections. A necessary condition for the vanishing of the 
semicircular arc, which, in turn, determines the choice among the three values for a, is that 
the total angle j3 of the integrand's "creeping exponential" exp{iz//3(z/)} (including the 
terms resulting from the Hankel functions) must be positive [and large enough to exclude 
the penumbra region in the "optical shadow" and "optically lit" region] for v' given by 
Eq.( [4.2|) , i.e., v' ~ kaj>. A violation of this condition would correspond semiclassically to a 
negative creeping path which has to be excluded for physics reasons: during the creeping 



the modulus of the wave has to decrease and not to increase |29|| , as tangential rays are 
continuously ejected, while the path creeps around a convex bending. The positivity of 
the creeping exponential actually only guarantees the vanishing of the integrand on the 
arc to the left of the line of zeros v\ of the Hankel function H^\kaji) and to the right 
of the line of zeros v'^ of the Hankel function H^\kaji) in the upper complex z/'-plane. 
The vanishing of the remainder of the arc is a consequence of the strongly decreasing 
J v /(kaji) / H^) {kaji) term which dominates the behavior of the integrand to the right of 
the z/^'s and to the left of the z/^ 2 ^'s. Whereas the v'^ line does not cause any problems, 
the v'z line is potentially dangerous as the Hankel function in the denominator is vanishing. 



The remedy is to put the path right in between two adjacent zeros |29 



As already mentioned, the "sin" -independent integral is symmetric in path and integrand. 
Because of this the path can be symmetrically deformed as follows (the preserved sym- 
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Fig. F.6. The original and the deform ed co ntour of the "sin" -independent integral for the case 
a' = 0. The lines of zeros are as in Fig, [F,5| . 

metry takes care of the original principal value description): It is replaced by an arc 
from — oo + ie to +ioo(l + i5), a straight line from +ioo(l + i5) to — ioo(l + i5) and fi- 
nally a symmetric arc (to the first one) from — ioo(l + i<5) to +00 — ie, where in the case 
a' = 0, the parameter 5 is chosen positive and small enough such that |Rez/| -C kay. 
[This allows later to use the Debye approximation of the Hankel functions, H^\kaj/) 
and H^\kaj/).] See Fig. [K.6| . The deformation of the path is justified as the sum of the 
new path and the (negatively traversed) original one do not encircle any singularities of 
the integrand. Since the integrand is symmetric under the exchange v' — > —u', the in- 
tegrals over the two symmetric arcs completely cancel, such that only the straight line 
segment from +ioo(l + i<5) to — ioo(l + i5) gives a contribution. This expression is fi- 
nite since it is symmetric under v' — > —v' and since the integrand vanishes rapidly for 
\v\ — > 00, as long as the slope of the straight line section is negative. In the case a' = 1 
the parameter S has to be chosen negative since the integrand only vanishes rapidly for 
a straight line section with positive slope (see Fig. [F.7| ). The reason for this difference 
is the presence of the ratio H^\kaj>) / ' H^,\kaji) in the first case which is replaced by 
(H^\kaji) / H^) (kaji)) — 2J l/ /(kaj') / H^\kaj') = —1 in the case a' = 1. (See also below 
the discussion of the pertinent Fresnel integrals in the semiclassical saddle-point approx- 
imation). 

As mentioned, the actual result depends on the concrete geometry and on the impact 
parameters l/k and l"/k, i.e., on the value of Acty = atjiiji — ctjj', on the value of a and 
the angles 4>±{l,i'') = arccos((z/ ± l)/kRjjn) and 4>'±{l" , v') = arccos((z/ ± l")/kRff) 
resulting from the asymptotic Debye approximation of the Hankel Hi± u ,(kkRjj>) and 
H±li_ v ,(kkRjijn), respectively. Since a can take three values there exist three mutually 
exclusive alternatives: 

The first one corresponds to a = and < 27r =f Aaji — <f)±(l, kaj>) — (f)±(l", ka^) < 2tt 
(this geometry allows only a classical path from disk j [under the impact parameter / / k] 
to disk j" [under the impact parameter I" jk\ which is specularly reflected to the right at 
disk see Fig. [F.l|) : 
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Mv = w if { Z 

„=o^ Hy(kaj> 



HP,{kR jf )H$>_ v ,(kR j , jH )e 



(i) 



(7r+ Aa •/ +2n7r) 



u dC+ 



+ H&(kRj j ,)H?>,_ v ,(kR jfj »)e 



(i) 



3 if'(7r— Act ./+2n7r) 



-ioo(l+i<5) ( 2 ) 



+ 



dz/ 



ioo(l+i<5) 



Hl}\ka f 



(kRjf ) H$L V , (kRfju ) e 



(i) 



— 7T+A« ■/) 



+ HilUkR jf )H^,_ l ,,(kR rf ,)e 



(i) 



— 7T — Ady ) 



(F.15) 



The second case is cr = 2 and < Act.,' — </>_(/, fca,/) — <p'!_(l" , kaf) < 2tt and < 
47r — A«j/ — (f) + (l,kajr) — (f)'+(l" , kcij') < 27r (this geometry allows only a classical path 
from disk j to disk j" which is specularly reflected to the left at disk f, see Fig. [F.2| ): 



HtUkRjfWLAkRff^e 



(i) 



_,ij/(A« •/+2rwr— 7r) 



^(fci2 ii 0Hi 1 i_ lw (A ; ^ i »)eK(^^+*-) 
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ioo(l+i<5) ( 2 ) 



dz/ 



+ioo(l+i5) 



H$\ka f 



H^UkR^H^ikR^e 



(i) 



3 ii/(-37r+Aa,-/) 



HlllikRj^H^ikR^e 



(i) 



3 ii/'(7r— Aa.() 



(F.16) 



The third alternative is cr = a' = 1 and < A«j/ — /ca.,/) — <p'L{l", kajr) < 2tt and 
< 27r — A«j/ — (f>+(l, kcij') — <f>+(l", kajr) < 2tt (this geometry allows only a "classical" 
path from disk j [under the impact parameter l/k] to disk j" [under the impact parameter 
I" /k] which goes directly through disk j', see Fig. F.3|) : 
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()..,/. Hy(ka 



K 1 ) (lcTt...\Tl { - 1 ) .JhTl \ p i^(Aa,,+2ri7r-7r) 



-ioo(l-i|<5|) 
ioo(l-i|<5|) 

+ H^UkRj^y^kR^) <yi*-*°v)] } . (F.17) 



Note that Eg. (|F.17|) can also be derived from the "+Aa,-/ part" of (F.15|) plus the "— Actf 
part" of (|F.16|) , by a rearrangement of the corresponding creeping and geometrical terms, 
i.e., by the addition of an extra term of smaller creeping length than the smallest one 
before and the subtraction of the very same piece from the geometrical terms. 

The contour integrals of these three alternatives are evaluated at the zeros u' e of the Hankel 
functions Hj)\ka,j>), such that 

di/ J ^ haf) Y v e^+iX- = 2vri £ - Y Vt e^+D^ . (F. 18) 

H u , (ka r ) i=1 g^H u , (ka r )\ v , =u > t 



Up to this point all expressions are still exact. The steps introduced so far just served 
the purpose of generating the three distinct "classically" allowed angular domains and 
of transforming the original expression ( F.3| ) into a form ready for the semiclassical ap- 
proximation. This will be taken next under the condition Refcaj' ^> 1. Note that this 
inequality automatically induces RekRjji ^> 1 and KekRj'j" 1. 

The contour integral (which, in fact, is now a sum over the residua) and the straight line 
integral are now treated semiclassically in different ways. 



F.3 Semiclassical approximation of the straight-line integrals 



The straight line integrals will be evaluated in the saddle-point approximation at a saddle 
v' s where the path crosses the real axis. For evaluating the saddle-point integral, the 
Debye approximation QE.ll ) will be inserted for the given Hankel functions. For the first 



and second alternative, an internal consistency check on the validity, is the condition 
\y' s jk\ < ay which in physics terms means that the impact parameter at disk j' has 
to be smaller than the disk radius a'j. For the third alternative the weaker conditions 
|(z/ — l)/k\ < Rjji and \{y' s — l")/k\ < Rj'j>> are sufficient: The difference in the impact 
parameters at successive disks should be smaller than the distance between the disks. Its 
validity is guaranteed by the triangular condition. 

The saddle-point integral is evaluated by expanding the exponents of the Debye approxi- 
mates to second order and a successive integration. The reflection angle is determined by 
the saddle-point condition itself, the geometrical length of the path can be read off from 



108 



the total exponent at zeroth saddle-point order, i.e. from the sum of square root terms of 
the Debye exponents divided by k. Under the Gauss' integration the second-order fluc- 
tuations about the saddle determine the stability factor \j\J i? eff and, together with the 
already present phases, the overall phase. 

The straight line integral of the first two alternatives corresponds then to the standard 
geometrical path from disk j [under impact parameter l/k] to disk j" [under impact 
parameter I" jk\ where there is a specular reflection from the boundary of the disk ay 
either to the right for the first alternative (Fig. [F.l| ) or to the left for the second alternative 
(Fig. |F.2| ). The slope of the path of this straight line integral, which asymptotically is 
i(l+i5), has to join smoothly the slope of the saddle-point path. This condition determines 
the sign of the slope. The saddle-point integral, which is of the Fresnel-type, results when 
the pertinent exponents of the Debye-approximated Hankel functions are expanded to 
second order around the saddle point v s : 



J dtfz/'e l2{&V) We-) 2 -' 2 V^W 2 -''-"'* 52 V (fc ^v )2 - ( ' A - i " )2 v' . (F.19) 



-e-W-^oo 



Here Su 1 = v — v s is the integration variable. By the substitution 8v' = e m ^x this 
integral becomes a standard Gauss' integral 



c 



i7r/4 / dx e- xH/2 = e~' mli J2^Jb (F.20) 



with b = 2{(ka) 2 - v'\Y x l 2 - {{kR jf ) 2 -(I- u's) 2 }- 1 ' 2 - {(kR fj ,,) 2 - (u' s - I") 2 }- 1 / 2 
positive as ay < Rjjr — aj,Rj/j" — af. The right hand side of Eq.( |b'.2(J| ) together with 
the prefactors and phases of the Debye-approximated Hankel functions determine the 
remaining terms (see below). 

Perturbative higher-order ^.-corrections (see Refs.[38-f40||) result here from higher-order 



terms in the Debye approximation through expansion terms proportional to (l/kr) n = 
(h/pr) n (with r = a,ji, Rjj> or Rj>j») and from the integration of polynomial second- and 
higher-order (fiu 1 /pr) 2n terms under the Gauss-type saddle-point integral. The polynomi- 
als are generated by a consistent expansion of all prefactors and exponential terms of the 
Debye series up to a given order. The Debye series reads 



H^ 2 h 



z) ~ 




1 - ^ 



i 

2\ 4 



3 ±iV z 2 — y 2 q=ii/arccos(-)=Fi j 



x < 1 Ti 
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(F.21) 

where the upper signs apply for the Hankel function of first kind and the lower ones for 
the Hankel function of second kind. Note that this expansion is of asymptotic nature and 
therefore induces the asymptotic nature of the ^-expansion itself. Here we will limit our 
discussion just to the leading term, such that no ^-corrections arise. 



F.4 Semiclassical approximation of the residua sum 



In the contour integral (or residua sum) the Debye approximation is not justified for the 
ratio 



JAkaj 



(F.22) 



since Rez/^ ~ kay. It is still valid, however, for the Rjj> and Rj' 3 .//-dependent Hankel 
functions, since Rjf > aj + aji and Rff > aj'+aj". Instead, the ^/-dependent cylindrical 
functions are evaluated under the Airy approximation. 

The latter step is justified as we evaluate the ratio ( |F.22| ) at the zeros z/£ of the Hankel 
function H^(ka). In the Sommerfeld representation the contour of a Hankel function 
H^(z) has normally two saddles [25]. For Re v -C z or Re v 3> z one of these two saddles 
dominates such that the Hankel function can be approximated by a single exponential 
times polynomial corrections (the Debye approximation). In such a case the Hankel func- 
tion can obviously not vanish. In order for it to vanish the contributions of the two saddles 
have to be of the same magnitude. In other words, we have to be in a region of competing 
saddles, where the standard saddle-point approximation (which has been used for the 
purely geometrical calculation) is replaced by the Airy approximation. This is the case 
when Rev ~ Re z. There we have, to leading order {y = v' and a = ay), 



£rW(te)~-e- 1 * (AV A(qV) + 0((ka)- 1 ), (F.23) 



7r \ka 

(ka) ~ - e +i f ( 3 A(qW) + O Uka)- 1 ) . (F.24) 
7r \ka I v ' 



Here 



q V)= e -i%(JL\ (u - ka) + O ((ka)- 1 ) , (F.25) 



= e +i f ( [y _ ka) + O ((ka)- 1 ) (F.26) 



are the zeros of the Airy integral [29 
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A(q) = J dr cos(gr — r 3 ) 
o 

= 3- 1/3 7rAi(-3- 1/3 g) , (F.27) 



with Ai(z) being the standard Airy function p0| ; approximately, qi m 6 1 / 3 [37r(£— 1/4)] 2//3 /2 



Thus Eq.( f4.2| ) arises. Note that this is the first term in an asymptotic expansion where 
the corrections are of relative order O (( ka)~%J ~ O (^ 3 ) as = p/fi. The first correction 
to the Airy approximation is therefore more important than the first correction term to 
the Debye approximation as the latter term scales as O (K). Up to order O ((A;a)~§), the 
zeros Ui read as 



i , 
ka\ 3 ; tt / 6 \ 3 qj ] 



v e ~ ka + es I I % e 3 f j ^ . . M ^ 



6\§ 1 (2% t 2%\q\ 



The //-derivative of the Hankel function H^{ka) at v — ui has the form [^J 



9 rrm/ , 2e _i Tr / 6 \§ ,,. . f e l f ( 6 \t q £ / 6 \s 37 /gA 2 



6 V / 37 563 
W V 36 • 5 • 630 ~ 5-630-9 V"6 



(F.29) 



where A'(qe) is the derivative of the Airy integral A(q) at the position qi |[29|| . The Airy 
approximation to the Bessel function J Ue (ka) = \H < ^{ka) reads as 



~t (=)' I 1 " * (e) * (I) " rf n (e) * (I 



fca/ \ 45 2 -7 6 V45/ 



(F.30) 



Applying the Wronsky-relation A(z)A'(ze ±27Ti / 3 ) - A'(z)A(ze ±27Ti / 3 ) = -f e^/ 6 one gets 
for 2 = qf ] (with = 0) 



(2)> 7T e '6 



m>) = 7r—7T)7- (F-31) 
6 A'{qf ] ) 



Thus, under the Airy approximation, each of the residua in Eq.( |K.18 ) becomes 



2vri J 4^r) y e Vn+iy e , = _ e -we Ct-K-VHka;,) 1 /* -Y u , e i( - 2n+1 >l« (F.32) 

ib R v '{ka f )\ u , =u > * 
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with the coefficient Cg = Cg(ka) 




(F.33) 



The values of the first zeros qe and the corresponding coefficients Ce, truncated at order 
0({ka}- 2 ' 3 ) = 0(h 2/3 ), can be found in Ref.@ and are listed in Table Q 

Table F.l 

The first zeros qe of the Airy Integral A{q) and the corresponding coefficients Ce of the creeping 
wave under Dirichlet boundary conditions in the leading Airy approximation. 



i 


qe 


c t 


1 


3.372134 


0.91072 


2 


5.895843 


0.69427 


3 


7.962025 


0.59820 


4 


9.788127 


0.53974 



We will limit our discussion to the Airy expansion of this leading order, i.e., v\ as in 



Eq^O) and Cg as given by the first term of Eq.( F.33 ), since all the higher terms vanish 



at least as fast as h 1 ^ 3 and h 2 ^ 3 , respectively, in the limit h — > 0. 

Finally, the Debye approximation ( |h].ll| ) is inserted in Y v > for the Rjj> and Rj'j" dependent 
Hankel functions. The two square root terms in the exponential of the Debye approximate, 



e.g. ^J(kRjji) 2 — (I — u' e ) 2 } etc., under the approximation v' t ~ kaf, give the length of the 
two straight sections of the path times k. All exponential terms proportional to e.g., 
i^arccos(- ■ •), rmv'^ correspond to the creeping sections (of mode number £) of the path. 
The latter include, of course, the creeping tunneling suppression factor linked to the 
imaginary part of the vi. The product of the two Debye prefactors is just the stability of 
the path times —i2/ir. The latter factor cancels the exposed factor in ( |F.32| ). 



In summary, the residua of the contour integrals in the Airy approximation correspond to 
those paths from disk j [under impact parameter l/k] to disk j" [under impact parameter 
I" /k] that have straight sections and circular creeping sections of mode number I which 
join tangentially at the surface of disk f . For the first term of Y Vv the creeping is in the 
left hand sense and for the second term in the right hand sense around disk j. The sum 
over n counts n further complete creeping turns around this disk. Note that the smallest 
creeping angle is less than 2ir, but larger than zero (see Figs.F.l-F.3). 



F.5 Resulting Convolutions 

The first alternative (Fig. |F.l| ) reads now 



V A jj ' A j ' j " 
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-^"EE«- W6 -- 1/2 cv(W) 1/3 



t=\ n=0 



X 



exp^(kR jf y-(l-v' e y exp^(kR ff ,) 2 -(v>-l») 2 

x e -i(Z-^)arccos[(«-^)/fcR J y]-i(^-Z'0arccos[(^-r0/fcfl J , j ,,]+i^(Aa J ,+(2n+l)7r) 



exp i<J (kRjj') 2 —(l+u' e ) 2 exp i^J(kRj'j") 2 — {y'f+l 



i(kR Jf y-(i+v>n/* [(kR fj „) 2 -M+i»n/* 

x e -i(H^)arccos[(^^)/fcK J y]+i(^+rOarccos[(-^-rO/fci?yy/]+i^(-Aa J / + (2ri+l)7r)} 

+ W j l" — e~ i7r//4 e"" arccos [('-^)/*:-Rjj']+««" arccos[(i^-H ARjV'] 
H V 7T 



x exp {i + iy/ikRjrf-ii/.-l") 2 ~ 2^{ka r ) 2 -v> 

x[(A^) 2 -^]V4{[2(^ 

- [(H,) 2 -^ 2 ] 1/2 {[(fc%0 2 -(^K) 2 ] 1/2 + [(^,) 2 -(^-r)T /2 }" 1/2 • 

(F.34) 

Here is given as in Eq. (|4. 2|) and as in Eq.( [F.33| ). The value of 1/ follows from the 
saddle-point condition 

Aaj' + 2 arccos[z///caj/] — arccos[(i/ — l)/kRjji] — arccos[(z^ — l")/kRjiju] = .(F.35) 

which fixes the scattering angle Qy = arcsin[v£/A;ajv] as 

^j/ = Aaj' + arcsin[(i/ — l)/kRjf] + arcsin[(i/ — l")/kRj'j"] . (F.36) 

One might wonder why there do not appear two different geometrical segments corre- 
sponding to the two terms of the straight line integral in Eq.( [F.15|) . The answer is that 
the second term of this integral gives the same contribution as the first one, since the 
values of the pertinent saddles just differ by a minus sign. [In fact, it is easy to show with 
the help of the transformation laws (F.8|) and (F.9|) that the second term of the straight 
line integrals is identical to the first one.] The effective radius R' eS belonging to Eq. ( |F.34|) 
results from the prefactors of the Debye-approximated Rjj'- and -RjY'-dependent Hankel 
functions, combined with the r.h.s. of ( |F.20p , and reads 

R/ _ 2 djj' d 3'j" - Pj'( d jf + dj'j") 3 y\ 
eS Pi' 

with 



d jf = jR 2 -({l-u' s }/k) 2 (F.38) 



d ff/ = ^R ffl -({K-l"}/k) 2 (F.39) 
Pf = ^4 - (u'Jk) 2 . (F.40) 
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This should be compared with effective radius generated by the standard evolution of the 
curvatures in the corresponding classical problem (see Eqs.( |5.56| ) and (|5.57 )) 



Rett — + 



(F.41) 



i=l 



Here L^+i is the length of the leg between the i th and the (i + l) th reflection. The quantity 
Ki is the curvature just after the i th reflection, i.e., 



K.j 



^_ 1 + L i _i i i r^cosi 



(F.42) 



where, in turn, rj and 0j are the local radius of curvature and the deflection angle at the 
2 th reflection. [Note that Kq 1 = 0.] By identifying L ,i = djji — Pj', L 12 = dj'j" — Pf, 
Ti = aji and (pi = 6j> (such that py = ay cos 9 f) one can easily show that R' ef{ and R e s 
give the same result. 

The result of the second alternative (Fig. F.2| ) is as in Eqs. QF.34|) and (F.35|) with Aoy 
replaced by Actjr — 2n. The third alternative (Fig. F.3|) reads as 



V \ J ' J " 



t=\ n=0 



expiJ{kR jf ) 2 -(l-v' e ) 2 expiJ{kR rr ) 2 -(iy' f -l") 2 



[(kR jf ) 2 -(l-u' e ) 2 ]^ [{kRffl ) 2 -(u'-l") 2 }^ 

x e -i(l-v' e )a.Tccos[(l-v^)/kR jj ,]-i(v' t -l")8xccoB[^^^^ 



expi v /(fc%/) 2 -(/+z/;) 2 expi^(kR fj ») 2 -(is' e +l") 2 
+ [(kR jf ) 2 -(M e ) 2 ]^ [(kR fj „) 2 -(v> e +l'W 

x e -i(l+v' e )a,rccos{(l+v' e )/kR ]j ,]+i(v' e +l'')a,rccosl(^ 

-\-Wjij," J— e~ i7r//4 e - i/arccos [('-^0/ fefi jj']+ i ''' arccos [(^- i '0/ fc R J 'j»] 



x exp {i ^(kR n/ ) 2 -(l-u' s ) 2 + ^(kR jf ) 2 -(v'-l») 2 } 

^J[(kR Jf ) 2 -(i-K) 2 } 1/2 + [(kR ffl ) 2 -(w-i") 2 }y 2 ' 

Here v' s has to satisfy the saddle-point condition 

Aa r - arccos[(i/ - l)/kR jr ] - arccos[(z/ - l")/kR rj „] = . (F.44) 

Again, the two terms in the straight line integral of Eq. (F.17f) give the same contribution, 
as the saddle 1/ of the latter term is — v' of the first one. The minus sign in front of the 
straight line integral is cancelled by an additional minus sign [relative to alternative one 
or two] resulting from the positive slope of the straight-line section (see Fig. |F.7| ) and the 
corresponding changes in the Fresnel integral 
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-ioo(l-i|i|) 



e -2W4 / dJz/ / e 2 ■ VV^y) 2 -^-^' 2 v/(^v') 2 -K-'")V . (F.45) 

ioo(l-i|5|) 

The latter, by the substitution 5v' = e 17r//4 x, becomes a negatively transversed Gauss' 
integral 



-m/4 



y Ax e~ xHl2 = -e-' m/i Jh^b (F.46) 



where 6 = {(kRjff -{I- y' s fY 1/2 + {(kRj/jn) 2 - - /") 2 }" 1/2 . In fact, all dependence 
of the disk f is finally gone from this expression. If the third alternative exists, the 
pertinent straight line integral corresponds to a "ghost" segment starting at disk j [under 
the impact parameter l/k] and ending at disk j" [under the impact parameter I" jk\ which 
is equivalent to the corresponding geometrical segment of the direct term A"„ (j" ^ j). 
Because of the angular conditions, specified before Eq.( [F.17p , the ghost path has to cut 
disk j', i.e. the modulus of the impact parameter v' s jk has to be smaller than the disk 
radius ay (see Fig.|F.3|). 



F. 6 Ghost segment 



Let us now discuss the "ghost" segment, i.e., the non-creeping terms of Eq. (|F.43|) . The 
ghost cancellation presented here is, of course, related to Berry's work on the ghost cancel- 
lation for periodic orbits in the Sinai billiard, see Ref . . However, here the calculation 



is based on Watson's method which specifies the integration paths, the signs of the ghost 
contributions and encodes the geometries (the choice of the three alternatives for a) into 
the creeping orbits. 

After restoring W^,, it reads 



aj"Hi„(kaj 



exp {i ^ k R 3f y-{l-v' s y + i^kR ]r f-{v'-l"Y} 

j[(kR 3jl y-(i-v> s ) 2 } 1/2 + [(kR j > j »)*-te-n 2 ] 1/2 



x e 
x e" 



+U (a ■/ -Q ■// i -arccos[(i-i/j)/fc_R . ■/]) 



il"(a ■; jfi —0.AAI1 —a,rccos[(u' s —l")/kR ,•/,•//]) 



with 



ozj»ji — o>jji + arcsin[(z/ — l)/kRjf] + arcsin[(z/ — l")/kRj>f] = n 



(F.47) 



(F.48) 



which is equivalent to condition ( F.44 ). As this saddle-point condition implies that the 
impact parameter v' s jk at disk j' lies on the straight line joining the impact parameter 
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l/k at disk j, with the impact parameter l"/k at disk j", the following relation between 
the lengths of the segments on this line holds 



R]f - {(l ~ <)/k} 2 + - i« - l ")/k} 2 = - {(l - l")/k} 2 , (F.49) 



i.e., the length of the straight line from the impact parameter l/k to the impact parameter 
I" /k is the sum of the lengths from l/k to v'/k and from u'/k to I" jk (see Fig. [F.3|) . 



The "ghost" segment ( |K.47 ) should be compared with Eq.( |b'.l| ), in the semiclassical ap- 
proximation ( [E.ll| ), for the Hankel function H^^kRjj//) 

A gr ~ (i - s jfl ) a ^ a ^ (-if e i(^-^,.,,) 

aj"Hi,,(kaj'/) 



x \l - e ~ W4 n 211 /l e ~ i( '~ n arccos[( ^" )/fci V ] . (F.50) 



Condition (F.49) implies that the lengths and stabilities of the ghost segment ( |b'.47| ) and 



of the direct path ( |F.50|) are the same. The comparison of the phases implies the relations 

7r/2 = (Xj>j — ctj"j + arccos[(i/ — l)/kRjj'] + arcsin[(/" — l)/kRjfr) (F-51) 
7r/2 = aijj» — cxjijn + arccos[(i/ — l")/kRjijn] + arcsin[(/ — l")/kRjju] (F-52) 



which are valid under the condition ( |F.48|) . Thus, we finally have in the semiclassical 
approximation 



ghostJ^K) = (Af ost Agl)^ ; ~ AgT (F.53) 
under the condition, of course, that the saddle v' s satisfies Eq. (|F.48|) . 
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Figures of 3-disk resonances 
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Fig. G.l. The Ai resonances of the 3-disk system with R = 6a. The exact quantum- mechanical 
data are denoted by plusses. The semiclassical ones are calculated up to 1 st order in the curvature 
expansion and are denoted by crosses: (a) Gutzwiller-Voros zeta-function (|2.1| ), (b) dynamical 
zeta- function fl2.2|), (c) quasiclassical zeta function (|2.3|), 
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Fig. G.2. The Ai resonances of the 3-disk system with R = 6a. The exact quantum-mechanical 
data are denoted by plusses. The semiclassical ones are calculated up to 2 nd order in the cur- 



vature expansion and are denoted by crosses: (a) Gutzwiller-Voros zeta- function (2A), (b) 
dynamical zeta- function Q2.2|), (c) quasiclassical zeta function (|2.3[). 
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Fig. G.3. The Ai resonances of the 3-disk system with R = 6a. The exact quantum-mechanical 
data are denoted by plusses. The semiclassical ones are calculated up to 3 rd order in the curvature 
expansion and are denoted by crosses: (a) Gutzwiller-Voros zeta-function (|2.1| ), (b) dynamical 
zeta- function (p.2|), (c) quasiclassical zeta function (|2.3|), 
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Fig. G.4. The Ai resonances of the 3-disk system with R = 6a. The exact quantum-mechanical 
data are denoted by plusses. The semiclassical ones are calculated up to 4 th order in the curvature 
expansion and are denoted by crosses: (a) Gutzwiller-Voros zeta-function (|2.1| ), (b) dynamical 
zeta- function fl2.2|), (c) quasiclassical zeta function (|2.3|), 
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Fig. G.5. The Ai resonances of the 3-disk system with R = 6a. The exact quantum-mechanical 
data are denoted by plusses. The semiclassical ones are calculated up to 5 th order in the curvature 
expansion and are denoted by crosses: (a) Gutzwiller-Voros zeta-function (|2.1| ), (b) dynamical 
zeta- function fl2.2|), (c) quasiclassical zeta function (|2.3|), 
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Fig. G.6. The Ai resonances of the 3-disk system with R = 6a. The exact quantum-mechanical 
data are denoted by plusses. The semiclassical ones are calculated up to 6 th order in the curvature 
expansion and are denoted by crosses: (a) Gutzwiller-Voros zeta-function (|2.1| ), (b) dynamical 
zeta- function fl2.2|), (c) quasiclassical zeta function (|2.3|), 
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Fig. G.7. The Ai resonances of the 3-disk system with R = 6a. The exact quantum- mechanical 
data are denoted by plusses. The semiclassical ones are calculated up to 7 th order in the curvature 
expansion and are denoted by crosses: (a) Gutzwiller-Voros zeta-function (|2.1| ), (b) dynamical 
zeta- function fl2.2|), (c) quasiclassical zeta function (|2.3|), 
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Fig. G.8. The Ai resonances of the 3-disk system with R = 6a. The exact quantum-mechanical 
data are denoted by plusses. The semiclassical ones are calculated up to 8 th order in the curvature 
expansion and are denoted by crosses: (a) Gutzwiller-Voros zeta-function (|2.1| ), (b) dynamical 
zeta- function fl2.2|), (c) quasiclassical zeta function (|2.3|), 
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Fig. G.9. The Ai resonances of the 3-disk system with R = 6a. The exact quantum-mechanical 
data are denoted by plusses. The semiclassical resonances are calculated up to 12 th order in the 
curvature expansion and are denoted by crosses: (a) Gutzwiller-Voros zeta-function ( |2.1| ), (b) 
dynamical zeta-function Q2.2|), (c) quasiclassical zeta function (|2.3[). 
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